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It is shown that if the equivalent magnetic dipole describing the earth's magnetic field at large 
distances is placed at the magnetic center of the earth as calculated by Schmidt, the theory of 
the latitude and azimuthal effects developed by Lemaitre and Vallarta also accounts for the 
longitude effect found by Clay, van Alphen and Hooft and, independently, by Millikan and 
Neher. Further it is shown that if the distribution function of magnetically deflectable par- 
ticles capable of reaching the earth's tropical belt is substantially an exponentially decreasing 
function of their energy, the theory quantitatively agrees with experiment. 


_Fpuad van Alphen and Hooft' and, inde- 
5 pendently, Millikan and Neher® have found 
experimentally that the total intensity of cosmic 
radiation measured by an ionization chamber at 
sea level is different for the same geomagnetic 
latitude and different longitudes. They have 
suggested that this longitude effect is to be 
interpreted as due to an asymmetry of the earth's 
magnetic field about the earth’s center. The 
purpose of this note is to examine this point of 
view in the light of the theory of the latitude and 
azimuthal effects developed by Lemaitre and 
Vallarta,’ and to find out what conclusions can 
be drawn from the experimental data concern- 
ing the energy distribution function of magneti- 
cally deflectable particles capable of reaching the 
tropical belt. 

The magnetic field at the surface of the earth 
(and elsewhere) can be represented by an ex- 

1]. Clay, P. M. van Alphen and C. G.’t Hooft, Physica 
1, 829 (1934). 

?R. A. Millikan and H. V. Neher, Phys. Rev. 47, 205 
(1935). 

* A discussion of the present status of the theory is given 
in a paper by Lemaitre, Vallarta and Bouckaert, Phys. 


Rev. 47, 434 (1935), where full references to previous 
literature are given. 


pansion of Legendre’ polynomials about any 
arbitrary point chosen as origin, in particular 
about the center of the geoid, but as Schmidt‘ 
has pointed out the expansion takes a particu- 
larly simple form, only one term of the first 
order and one of the second remaining if the 
origin is chosen at the magnetic center of the 
earth, a point about 300 km from the center of 
the geoid at geographic coordinates lat. 10 N, 
long. 168 E. If the equivalent magnetic dipole 
is placed not at the center of the geoid, but at the 
magnetic center of the earth with its axis point- 
ing towards the geomagnetic poles, the im- 
mediate consequence is that the distance from 
the geomagnetic equator, or from any geo- 
magnetic parallel, to the dipole is a periodic func- 
tion of longitude. Since changing the distance 
to the dipole also changes the opening of the 
allowed cone’ for any given energy, one must 

‘A. Schmidt, Zeits. f. Geophys. 2, 38 (1926). The in- 
troduction of the actual magnetic field of the earth, using 
Schmidt's results, as done in the present paper, was already 


suggested by Lemaitre and Vallarta, reference 6, p. 88, 
footnote (7). 

5 The theory of the determination of the cone is outlined 
in reference 3 and more fully discussed in Lemaitre’s and 
in Bouckaert’s papers there quoted. J 
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Fic. 1. The solid angle of opening of the cone at three 


geomagnetic latitudes. 


expect different intensities at different longi- 
tudes and consequently a longitude effect. We 
now show not only that this asymmetry of the 
earth's magnetic field is fully competent to 
account for the observations of Clay, van Alphen 
and Hooft, and of Millikan and Neher, but also 
that their intensity measurements along the 
geomagnetic equator throw valuable light on the 
energy distribution function of magnetically 
deflectable particles which give rise to both the 
latitude and the azimuthal effects. 

In Fig. 1 the solid angle cut off from the unit 
sphere (i.e., the sky) by the allowed cone is 
plotted as a function of energy, measured in 
terms of our energy parameter® x», for geo- 
magnetic latitudes 0°, 10° and 20°. The solid 
angle in each instance has been obtained by 
graphical integration from the calculations of 


TaBLe I. Calculated distances from the earth's magnetic 








APPROXIMATE 





Pornt GeomaG. Geocr. GEocrR. ’ 
No. LAT. LAT. LONG. (KM) Ax/xo LOCA PION 

1 0 128 76W 6519 0.0234 Lima, Peru 

2 0 0 160 W 6120 —0.0393 SSroute Honolulu 
to Sydney or 
Melbourne 

3 0 10 N 75 E 6390 0.0031 SS route Bombay 
to Colombo 

4 0 6S 1S W 6650 0.0440 SS route South- 
ampton to Cape 
Town 

5 0 12N 110E 6220 ~-0.0236 SS route Hong 
Kong to Singa- 
pore 

6 0 10S 100 W 6370 0 1000 miles SW of 
Galapagos 

7 —10 2N 104 E 6240 —0.0204 Singapore 

& 10 258 80 W 6520 0.0236 Guayaquil 

9 —20 6S 106 E 6220 —0.0236 Batavia 

10 20 9N 79 W 6460 0.0141 Panama 





*G. Lemaitre and M. S. Vallarta, Phys. Rev. 43, 89 
(1933), Eq. (12). This amounts to choosing as our energy 
unit 1/10.60 ergs, or 5.96 X 10" electron volts, for electrons 
or positrons. 
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Fic. 2. The solid angle of opening of cone at several points 
on geomagnetic equator. 


Bouckaert.’ In Table I we collect the results of a 
calculation of the distance from the earth's 
magnetic center to several points on its surface, 
together with the fractional change from the 
standard value of the energy parameter due to the 
fractional change in distance from the standard 
earth's radius (6370 km). 

From the value of Ax/x» given in this table 
and from Fig. 1 one easily finds the solid angle 
of the cone (the full light region in the language 
of previous papers) for the different points in the 
table. The solid angle is plotted as a function of 
energy, expressed in terms of our energy param- 
eter, for these points, in Figs. 2, 4 and 5. 

Referring to Fig. 2 it is seen that for any given 
energy the full light region at the geomagnetic 
equator is greatest along curve 4 (for the point 
on the earth where the steamer lane from South- 
ampton to Cape Town intersects the geomag- 
netic equator, approximately), is less for the 
same energy along curve 1 (Lima, Peru, ap- 
proximately), then less along curve 3 (where the 
steamer lane from Bombay to Colombo inter- 
sects the geomagnetic equator), curve 5 
(steamer lane from Hong Kong to Singapore) 
and least along curve 2 (steamer lane from 


7 L. Bouckaert, Ann. de la Soc. Sci. de Bruxelles, A54, 
174 (1934). Along the equator the angle of opening of the 
cone in the east-west vertical plane may be calculated 
from Stérmer’s formula sin 6= —2y,/x cos \+cos A/x* 
with \=0 and y, =1, because there the vertex of the cone 
touches Stérmer’s circular cone as defined by the above 
expression. The angle of opening of the cone in the east- 
west vertical plane may not be calculated by this formula 
for any other latitude, since for no other latitude does this 
condition hold, nor may the solid angle of opening of the 
cone be calculated from this formula for any latitude 
including the equator. 
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Fic. 3. The function ¢~*0’w(x9) for different points on 
geomagnetic equator. 


Honolulu to Melbourne). Consequently the total 
intensity must decrease in the same order. 

We now calculate the percentage change in 
total intensity due to the fractional change 
Ax/x» in the energy parameter xo, and show that 
the experimental results are consistent with the 
assumption that the energy distribution func- 
tion, defined as the number of particles outside 
the atmosphere having a given energy, is an 
exponentially decreasing function of the energy. 
Let f(x) be the distribution function and w(x») 
the function giving the solid angle of the cone 
in terms of xo. Then the number of particles 
coming through the cone corresponding to any 
particular value of the energy is, by Liouville’s 
theorem,® f(xo)w(x9). We now define the total 
intensity as the total number of particles which 
come through all the cones, beginning with those 
particles which have the least possible energy 
to reach the top of the atmosphere at the given 
point. This total intensity is given by 


co 
T= f(Xq) w(xXo)dxo. 
“ 2 
* A discussion of the application of Liouville’s theorem to 
the present theory, with full literature references, is given 
in reference 3. 


Let us now assume that the distribution 
function is an exponentially decreasing func- 
tion of energy. Remembering that for energies 
of the order of magnitude required to reach the 
tropical belt the energy is proportional to the 
square of the parameter xo, we have for the total 


intensity, 
1={ e~ 8 w(x9)d xo. 
r 


In Fig. 3 we have plotted as a function of x» 
the integrand of the above expression for the 
points in the geomagnetic equator given in 
our Table I. From this we have calculated by 
graphical and analytical integration the inten- 
sity J» at those points on the geomagnetic 
equator the distance of which from the magnetic 
center is equal to the earth's radius, i.e., the 
value of the above integral along curve 6 (or 3). 
This we have taken’ as our standard intensity 
with respect to which the percentage change of 
intensity due to the longitude effect is calcu- 
lated. To find the change of intensity at the 
different points in our Table I we merely have to 
find by graphical integration the areas of the 
strips between corresponding curves in Fig. 3. 
For points 1, 2, 3 or 6, 4 and 5, respectively, we 
thus obtain AJo/J>=1.9 percent, —3.7 percent, 
0, 4.2 percent and —2.5 percent, respectively. 
The value of AJ,/Jo is different for the same 
Ax/xo| because the areas of the corresponding 
strips are different. 

Similarly for points 7 and 8 at A= +10° the 
theory demands that the total intensity at the 
point on the earth described by curve 8 be 
greater than that described by curve 7 (Fig. 4) 
and the same calculation outlined above for 
equatorial points, with the same distribution 
function gives 5.6 percent as the percentage 
intensity change A/,/J,; between these points 
Likewise for points 9 and 10 (Fig. 5) at A= +20 
the theory requires that the intensity at point 
9 be less than that at 10, and a similar calcula- 
tion gives for the percentage intensity difference 
Al,/1,=4.0 percent. 

We now proceed to compare our theoretical 
results with experiment. It is very fortunate 
that the steamer lane from Bombay to Colombo 
intersects the geomagnetic equator at a point, 
the distance of which from the magnetic center 
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Fic. 4. The solid angle of opening of cone for two points 
at A=10°, 


is very nearly the standard earth’s radius. 
Consequently we have taken 1.73 as our standard 
ionization from which to reckon the intensity 
changes due to the longitude effect along the 
geomagnetic equator in the case of Clay, van 
Alphen and Hooft’s measurements. Conditions 
are not so favorable in the case of Millikan and 
Neher’s measurements because unfortunately 
they have not reported separately their measure- 
ments on the run from Hong Kong to Singapore 
and from Colombo to Bombay. However, the 
other point (point 6 in our Table I) which enjoys 
the same property of being at a distance from 
the magnetic center equal to the earth's radius 
is between Lima and the point where the steamer 
lane from Honolulu to Sydney intersects the 
geomagnetic equator. By interpolation we have 
taken as standard ionization in the case of 
Millikan and Neher's measurements the value 
2.60 and have referred the percentage change of 
intensity due to the longitude effect to this 
value. 

To compare the experimental value of the 
intensities at Guayaquil and Singapore we 
multiply the intensities given as ordinates in 
Millikan and Neher's® Fig. 2 by the factor 
0.1191 so as to bring their high latitude values 
in their Fig. 2 into agreement with their high 
latitude values in their Fig. 5. This gives at 
Singapore the ionization 2.51 which is 5.2 
percent less than the ionization at Guayaquil 
(2.64). 

If one now desires to compare the intensity at 
Batavia as measured by Clay, van Alphen and 
Hooft with the intensity at Panama as measured 
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Fic. 5. The solid angle of opening of cone for two points 
at A=20°. 


by Millikan and Neher one meets with a grave 
difficulty. Clearly the ionization measured by the 
former is not directly comparable with that 
measured by the latter, presumably because of 
differences of construction in their respective 
apparatuses. If one attempts to reduce the 
former's values to the latter's by multiplying by 
the factor (1.374) which is required to make the 
high latitude value of the former agree with that 
of the latter one obtains the value 2.45 at 
Batavia which gives about 8 percent less in- 
tensity there than at Panama, where according 
to the latter the ionization is 2.64. If one at- 
tempts the reduction by multiplying by the 
factor (1.395) which is needed to make the 
former's measurements at Singapore agree with 
the latter’s at the same point one obtains 
2.47 as the ionization at Batavia, which is about 
7 percent less than at Panama. It would appear 
that the conversion factor, besides depending 
on the apparatus, also depends on latitude. 
Further clarification of this point from the 
experimental side is earnestly to be hoped for. 
The calculated and experimental results are 
compared in Table II, and the calculated longi- 


TABLE II. Calculated and observed values of AJ,/ I; 








and Al, q;. 
Pornt Alo/le ‘Alo/le Ali/li Ahjli 
No. CALc. OBS. CALC. OBs. OBSERVERS 
i 1.9 1.2 Millikan and Neher 
2 —3.7 —O0.8 ‘ ” 
3 0 0 Clay, van Alphen and Hooft 
4 4.2 4.2 = . , , 
5 —2.5 —3.8(? Millikan and Neher 
7 —§.4 —5§.2 . = 
8 ; ; - 
9 - —~4.0 —7.0(?) | Clay, van Alphenand Hooft 


10 | Millikan and Neher 
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Fic. 6. The longitude effect along the geomagnetic equator. 


tude effect at the geomagnetic equator is 
plotted in Fig. 6 together with the experimental 
points. The curve is symmetrical about the 
longitude 12 W approximately, but asymmetrical 
with respect to the axis of the abscissae. 

With the exception of point 2 (SS route 
Honolulu to Melbourne at geomagnetic equator) 
it is seen that the agreement is very satisfactory. 

The influence of atmospheric absorption will 
be considered in a future paper, but even without 
taking this factor into account it is seen from 
the above discussion that the theory is sufh- 
ciently in agreement with experiment. An addi- 
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tional factor which may affect the longitude 
effect is the bulging of the earth's atmosphere 
along the geographic equator due to the earth's 
rotation. Estimates of the magnitude of this 
effect are difficult to give at the present time, 
but in general it is clear that absorption would 
be greatest along the geographic equator and 
least at the poles. It is possible that the in- 
tensity at Lima may be affected by one percent 
due to this factor, as compared with the intensity 
at the west coast of Africa (SS lane South- 
ampton to Cape Town at geomagnetic equator, 
point 4) but there seems to be no hope of ac- 
counting in this way for the high intensity 
measured by Millikan and Neher on the run 
from Honolulu to Melbourne, for there magnetic 
and absorption effects cooperate to give the 
least intensity. 

We conclude that the theory of the latitude 
and azimuthal effects developed by Lemaitre and 
Vallarta is competent to account quantitatively 
for the longitude effect, and further that if the 
energy distribution function of the cosmic radia- 
tion which is responsible for these effects is 
substantially an exponentially decreasing func- 
tion of the energy, the theory agrees with the 
measurements of the longitude effect. 

My thanks are due to Dr. L. Bouckaert and 
Mr. E. J. Schremp who have helped me with the 
graphical integrations and with the discussion 
of the distribution function suggested in this 


paper. 
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Difference Bands in the Spectra of the Major Planets 


ArTHUR ApeEL AND V. M. SiipHer, Unt 


rsil f Michigan and Lowell Observator 


Received March 11, 1935 


Chis brief paper contains the first identifications of difference bands in 


pl inetary spectra 


hese absorption bands lie at 816uu, 802uu, 683uu, 673uu and 584up~. They are due to methane 


' 


~ 


T is now well established that the atmospheres 
of the giant planets (Jupiter, Saturn, Uranus, 
Neptune) contain vast amounts of methane 
CH,) in a large excess of hydrogen, and that 


s and are found in the spectra of Jupiter, Saturn, Uranus and Neptune. 


in the two former ones there is also present a 


relatively small amount of ammonia (NH 


A. Adel an \ MSI pher Phys. Rev. 46, 92 1934): 
H N Russell scrence 81 l 1935 
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The temperatures' which prevail in these 
atmospheres are all well below — 100°C, so that 
the percentage of molecules occupying excited 
vibrational states is indeed small. The over- 
whelming molecular population residing in the 
ground state can be appreciated from the fol- 
lowing simple example. The ratio of the number 
of methane molecules excited thermally to the 
mode of vibration v»»= 1520 cm™', to the number 
of methane molecules in the ground state is 
very nearly given by 2e~*"/*", the factor two 
being present in virtue of the isotropic nature of 
ve in two dimensions. In the Jovian atmosphere, 
above the clouds, where 7 =163, this ratio 
becomes ~2e-"=1/200,000; that is, about five 
molecules in every million are thermally ex- 
cited to the vz fundamental vibration. 

Despite the small percentage population 
resident in excited states, the actual number of 
thermally excited molecules appears to be 
sufficient to render certain of these levels ap- 
parent in absorption. Consider, for example, the 
absorption band at 683." which appears in the 
major planet spectrum. We have observed it also 
in the laboratory, leaving no doubt that it is due 
to the methane molecule. In structure it is 
very narrow, contrasting markedly with the 
broad regions occupied by most of the thirty- 
six planetary methane bands already identified 
as having arisen in transitions from the ground 
state.' The spectrum given below (Fig. 1) shows 
the 683un band in the presence of several 
others. It was obtained in the laboratory with 
2140 meter-atmospheres of CH,, at 47.6 at- 


2 V. M. Slipher, Lowell Observatory Bulletin No. 42. 
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mospheres and 7 =293°. The narrow nature of 
the 683uy band is what might be expected in a 
band arising in absorption from an _ excited 
state whose population is completely insufficient 
for the development of the P and R branches 
and only barely adequate to exhibit a trace of 
development in the Q branch.’ Further evidence 
in support of this view is obtained from the fact 
that no plausible combination of frequencies 
can be found which would trace this band to a 
transition from the ground state, whereas it is at 
once correlated as a difference band; namely, as 
due to the transition from v2 to 6v3. This places 
the band 6v; at 6187A in excellent agreement 
with observation.' Besides 6v3;—v, at 683uy, 
there appears also in the spectrum of Jupiter a 
structurally similar band at 816uu.‘ In the same 
way, this band is to be identified as 5¥;— v2. By 
way of description, we may say very briefly 
that v2 is a symmetrical, double vibration, 
inactive in the infrared; and the sequence of 
bands mvs constitutes the strongest set in the 
major planet spectrum. 

When we come to consider the probable ap- 
pearance in the planetary spectra of other differ- 
ence bands, it appears that just as mv3+ »; gives 
rise to the strongest combination bands,' so 
the difference band sequence nv;—v; is the most 
likely to appear, in virtue of the great intensity 
of the set mv3. m is any integer, and i, of course, 
may be 1, 2 or 4.5 We have already discussed 

$In virtue of the spherical symmetry of the methane 
molecule all of its absorption bands possess P, Q and R 
branches 

*V. M. Slipher, Pop. Astron. 37, 140 (1929 

5 The methane molecule possesses the four fundamental 


frequencies of vibration: »=2915, ».=1520. , 3014 
v= 1304 cm". 
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Fic. 1. Absorption spectrum of methane showing the 683uu band in the 
presence of several others. 
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the system mv3;—v2; only the members 5y3— v2 
and 6v3—v2 have been located. The system 
nv;—v, may be expected to lie at approximately 
921 uu, 754uy, 644uu and 566uy for n=5, 6, 7 
and 8, respectively. 7v3;—», and 8v3—», cannot 
be located in the spectra of the major planets. 
However, the spectra of Neptune and Uranus, 
where the bands 5v3;—¥v, and 6v3—»,; are most 
likely to be discovered, are unknown in this 
region. 

We come finally to the group of bands 5v3— v4, 
6v3—vs4, 7vs—v, and 8v3—v4, which may be ex- 
pected to lie approximately at 802uy, 673yuy, 
584uu and 519uy, respectively. The band 5y3— 4 
at 802yuu is very prominent in the spectrum of 
Jupiter. In Uranus and Neptune, 6v3—¥», is 
hidden by the mass of absorption 5v3+»;,? while 
7v3—v, is definitely present at 584yy.’ 

It is evident that nv3;—v», is the strongest 
sequence with bands observable out to n=7: 
nv3;—v_ is second, with m reaching to 6; while 


s The absorption bands 5, 6, 7 and 8», are centered at 
approximately 726yuy, 619 uy, 543uu and 486uu, respectively. 
7 V. M. Slipher, Lowell Observatory Bulletin No. 13. 
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nv3—v, is very weak and, thus far, completely 
unobserved in the planetary spectrum.* 

The allocation of observed difference bands 
relative to the several planets is given in Fig. 2. 
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Fic. 2. Allocation of observed difference bands relative to 
the several planets. 





* The correlation of these five difference bands brings 
the total of identified methane bands in the major planet 
spectrum to forty-one.! 
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Provisional Wavelength Standards for the Extreme Ultraviolet 


J. C. Boyce anv C. A. Rieke, George Eastman Research Laboratory of Physics, Massachusetts Institute of Technology 
(Received March 7, 1935) 


This paper discusses the criteria for suitable wavelength 
standards in the extreme ultraviolet, and presents provi- 
sional values (obtained with a two-meter focus normal inci- 
dence vacuum spectrograph) for a number of lines of C, N, 
O and A falling in the spectral range between \1850 and 
\800. These are compared with the previous results of 
Edlén and of Bowen and of Bowen and Ingram. Further 


HE report of the commission on standards 

of wavelength of the International Astro- 
nomical Union contains the following criteria for 
auxiliary standards: “It is recommended that the 
mean wavelengths of the lines in certain groups 
of lines which give high orders of interference, 
and are satisfactorily distributed in the spectrum, 
and have been referred to the primary standard 


independent determinations are desirable. The present 
values, as those obtained by previous investigators, may be 
subject to a small systematic error inherent in the method 
of overlapping orders. This source of error is believed to be 
small, but must be eliminated before permanent standards 
can be established. 


by three observers whose measurements are in 
satisfactory agreement, be adopted as working 
standards.’""' Tertiary standards are further de- 
fined as those obtained by interpolation between 
these secondary standards. It is well known that 
no standards which come anywhere near meeting 


! Trans. Int. Astron. Union 1, 35 (1922). 
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TABLE I. Hydrogen. 








BOYCE AND C. A, 





RIEKE 


TABLE II. Helium (He II). 














x a 
OBSERVED Wr. ComMPUTED* 
1215.666 23* 1215.664 
1025.725 9 1025.717 

972.538 6 972.532 
949.740 4 949.739 
937.792 1 











these requirements have ever been established 
for the vacuum region of the extreme ultraviolet. 
Neither have interferometric determinations 
been made nor have three investigators obtained 
satisfactory agreement by any other method. The 
purpose of the present paper is to present the 
best values the authors have been able to obtain 
by the method of overlapping orders, together 
with a statement of the limitations of that 
method. It is believed that these preliminary 
values are more dependable in the region they 
cover than any previously published. It is hoped 
that independent work in this laboratory and 
elsewhere will soon confirm or revise them, and 
that in the interim they will be of use. 

Comparison spectra in this region of the spec- 
trum are most conveniently introduced by some 
form of gas discharge or by the presence of gases 
driven out by spark electrodes. Working stand- 
ards should therefore be sought among the more 
common elements easily excited in a gas dis- 
charge, namely H, He, C, N, O, Ne and A. Lines 
of these light elements may not be so sharp as 
those of metals; nevertheless their widths are less 
than the available resolution of most vacuum 
spectrographs. For special purposes another set, 
such as the lines of Cu II, would be desirable. 
The present discussion wil be restricted to lines 
of these lighter elements. 

Helium I and the various stages of ionization 
of neon are omitted, as their principal lines lie at 
shorter wavelengths than are considered here. 
While not literally satisfying the criteria quoted, 
the wavelengths calculated by Paschen? for the 
Lyman series of hydrogen and for the two series 
of ionized helium rest on a sound theoretical 
foundation and appear to be by far the most 
reliable in this region of the spectrum. For con- 





?F. Paschen, Preuss. Akad. (1929), 662. See also W. G. 
Penney. Phil. Mag. 9, 661 (1930). 


d d 
ComMPUTED* COMPUTED" 
303.7788 1640.409 
256.3145 1215.129 
243.0244 1084.940 
237.3297 1025.270 
234.3452 992.361 

972.109 
958.696 


949.326 


venience these values (for a few of the strongest 
lines of each series) have been included in Tables 
I and II. The spectral range covered by these 
lines is limited, and their higher series members, 
particularly those of helium, are hard to excite. 

Comparisons of these wavelengths have been 
made with values published by Edlén* and by 
Bowen‘ and by Bowen and Ingram.* The average 
difference between the wavelengths of Bowen 
and of Bowen and Ingram and the present pre- 
liminary standards is not more than a few 
thousandths of an Angstrom unit, though the 
individual values differ considerably, espe- 
cially above 41500 (here their wavelengths are 
usually given to only 0.01A). There is a good 
comparison with Edlén’s standard lines only be- 
tween A800 and 41200. The difference is constant 
and around 0.005A from A800 to 1050; it then 
increases to about 0.017A at 41200. Below A800, 
where we have compared with additional lines 
given by Edlén, we find the systematic difference 
to be less than 0.005A. We have made some use 
of Edlén’s values, corrected for this systematic 
difference, below \800 in drawing our reduction 
curves for neon and argon. The present article 
will therefore be restricted to the region from 
4800 to 41850. Adequate standards over such a 
range are urgently needed. When there is op- 
portunity to obtain further independent measure- 
ments in the region of shorter wavelength they 
will be published, since it is important that 
Edlén’s results obtained at grazing incidence be 
checked at more nearly normal incidence. 

The standards of Bowen and Bowen and In- 


?B. Edlén, Nova Acta Reg. Soc. Sci. Upsaliensis, Ser. 
IV 9, No. 6 (1934), and Zeits. f. Physik 85, 85 (1933). 

‘1. S. Bowen, Phys. Rev. 29, 231 (1927). 

+I. S. Bowen and S. B. Ingram, Phys. Rev. 28, 444 
1926). 
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gram, those of Edlén, and the provisional values 
proposed in the present paper, were all obtained 
by a comparison, direct or indirect, of high orders 
of the lines in question with first orders of the 
secondary and tertiary standards obtained from 
an iron arc. They are all subject to the possible 
error inherent in the use of overlapping orders. 
This source of error has been investigated by 
Michelson*® and by Kayser.’ The amount of the 
error varies greatly with the imperfections of the 
grating used. It is also very sensitive to focal 
adjustment, and a grating in which this error is 
known to be small may suddenly develop a large 
error of this type if its adjustment is even slightly 
disturbed. The plates used in the present in- 
vestigation are believed to be reasonably free of 
this type of error, to perhaps 0.005 to 0.010A. 
This belief is based upon a set of intercomparisons 
between different series of orders and upon the 
good agreement obtained with some of the wave- 
lengths of H calculated by Paschen.? 

The instrument used in this investigation was 
the two-meter normal incidence broad range 
vacuum spectrograph of the Carnegie Institution 
of Washington, which is located in the spectro- 
scopy laboratory of the Massachusetts Institute 
of Technology. It has already been described in 
a paper,® which discusses the method of reduction 
of its plates. Results with neon have been pub- 
lished ;* these depend to some extent on Edlén’s 
wavelengths below A800. For the wavelengths 
given in the accompanying tables, twelve expo- 
sures on five plates have been measured against 
the iron arc in the region 42250 to A2540, and 
three long exposures (on neon, argon and krypton 
plates) have been measured over their entire 
range. For the argon standards two additional 
exposures on argon were reduced, using the pre- 
liminary standards for C, N and O, which oc- 
curred as impurities in these exposures. Each 
plate was measured independently in both direc- 
tions by two observers. The iron comparisons 
provide standards for those lines whose higher 
orders fall in this region of the plate. The reduc- 
tion curves for other regions were drawn using as 


*A. A. Michelson Astrophys J. 18, 278 (1903). 

7H. Kayser, Astrophys. J. 19, 157 (1904); 20, 327 (1904 

*K. T. Compton and J. C. Boyce, Rev. Sci. Inst. 5, 218 
(1934). 

* J. C. Boyce, Phys. Rev. 46, 378 (1934 
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standards these lines and their intermediate 
orders, and on the krypton plate, lines of Kr II 
computed using the combination principle.” In 
addition the various orders of a given line were 
intercompared, and a few corrections were ap- 
plied to those sections of the plate which had a 
relatively small number of standards. The three 
plates used have, of course, been reduced inde- 
pendently. 

The probable error of a single observation has 
been computed from the differences of the meas- 
ures on lines in their first order, according to 
Peters’ formula. The lines used are those in the 
accompanying tables—all of shorter wavelengths 
than the iron comparisons. Hence the probable 
error includes the error in drawing the reduction 
curves in the regions where iron was not directly 
available; it is slightly less than +0.007A for an 
observation of unit weight (one measurement in 
the first order). No considerable sources of 
systematic error suggest themselves with the ex- 
ception of an error in the mean value of the 
twelve iron comparisons, and the possible error 
arising from the method of overlapping orders. 
The first of these might arise if the light from the 
arc failed to fill the grating, or if slight tempera- 
ture changes occurred during the long gas ex- 
posure. Due precautions were taken to reduce 
both these sources of error, and they are believed 
to be less than the accidental error. The probable 
error of a single iron comparison has been found 
from the differences between the twelve plates 
measured. It equals +0.007A for a single com- 
parison and hence +0.002A for the mean. Those 
lines measured against iron on fewer than the 
twelve exposures have been corrected to the 
mean system. Obviously, the accidental error of 
a line measured in the second order is one-half 
that in the first, etc. Some close multiplets were 
suspected of being affected by photographic 
errors in their first orders; for these lines the 
adopted wavelength depends only on measures in 
higher orders. The weight assigned to a meas- 
ured wavelength is the summation of the prod- 
ucts, number of plates times order. 

In Tables III to VI the adopted wavelengths 
are given, followed by the results of Bowen or of 


‘© T. L. de Bruin, C. J. Humphreys and W. F. Meggers, 
Bur. Standards J. Research 11, 409 (1933). . 
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TaBLe III. Carbon. 

A BOWEN Ep.Lén Int. Spectrum Wr. | a BOWEN EpLtn INT Spectrum Wr 
1930.900 5 I 3 (1194.08) 1— I 
1760.818 85 4 Il 3 1193,243 3 I 3 
1760.412 44 3 II 3 1193.012 1 I 3 
1658.126 13 3 I 3 1176.370 .359 351 10 Ill 20* 
1657.916 .92 3 I 3 1175.986 .988 .973 9 Ill 14* 
1657.380 37 2 I 2 1175.716 711 .700 15 III 18* 
1657.005 01 6 I 3 1175.582 577 5 Ill 16* 
1656.280 .27 4 I 2 1175.259 .261 .248 9 Ill 18* 

(1561.42) 381 8 I 1174.934 .922 .916 10 Ill 18* 
1560.699 .660 5 I 3 1141.745 1 Il 2 
1560.316 .267 4 I 3 1141.623 .61 Il 8* 
1550.790 174 10 IV 2 1139.343 3 Il 4 
1548.214 .189 12 IV 2 1066.138 5 II 6 
1463.33 3 I 3 1065.895 7- II 9 
1335.71 .705 .684 40 II 3 1037.020 021 .017 10 II 8 
1334.54 .539 515 30 II 3 1036.332 .336 .330 9 II 5 

(1329.58) .583 6 I 1010.374 .382 369 7 II 6 
1329.101 .100 5 I 2 1010.092 .090 .074 5+ II 6 
1328.831 .839 4 I 3 1009.862 .870 854 5 II 6 
1323.94 .916 15 II 2 977.020 .031 .026 30 Ill 6 
1280.340 2 I 2 945.566 3 ] 3 
1277.550 4 I 3 945.336 2 I 3 
1277.274 4 I 3 945.193 1 I 3 
1261.560 1 I 3 904.482 «472 468 s Il 9 

(1261.12) 1 I 904.144 133 .134 8+ II 6 
1260.955 2 I 2 903.952 .960 .950 8— II 6 
1247.387 391 .368 7 Ill 3 903.614 .620 .609 8 II 9 
1194.491 3 I 3 858.559 561 561 6 II 6 

5 Il 5 


858.091 .088 .094 








Bowen and Ingram, and those of Edlén. The in- 
tensity was estimated on an arbitrary scale. The 
state of ionization is denoted by a roman num- 
eral, and the final column gives the weight as- 
signed to the determination. The probable error 
of each wavelength is then +0.007A divided by 
the square root of the weight. Lines which have 
been measured in their higher orders against iron 
are indicated by an asterisk; a few below A800 
have been included. Some members of multiplets 
are included for identification purposes, although 
on account of blending they are not considered 
suitable for standards; these are given to two 
decimals and placed in parentheses. A few other 
lines are given to two decimals because the 
measurements are somewhat uncertain because 
of their great intensity or close companions. The 
computed hydrogen values are reprinted here for 
convenience and with them are our observed 
values of those lines for comparison. The com- 
puted values seem to us to be the more funda- 
mental and should be used as standards, though 
in practice the lines are unfortunately often so 
broad as to be unsuitable. 

The provisional values given here should be 
checked by other investigators. A few preliminary 


results provided by Professor Harrison in the 
current program of the Massachusetts Institute of 
Technology 21-foot normal incidence vacuum 
spectrograph indicate good agreement. Both the 
21-foot spectrograph and the broad range two- 
meter instrument have gratings ruled at The 
Johns Hopkins University. It is important that 

TABLE IV. Nitrogen. 











» BOWEN INT SPECTRUM Wr 
1745.246 .260 2 I 3 
1742.734 .740 3 I 2 
1494.669 4 I 3 
1492.630 5 I 3 
1243.297 3 I 8* 
1243.170 4 I he 
1200.706 681 4 I 19* 
1200.220 .200 5 I 22° 
1199.547 .533 6 I 23* 
1134.980 .987 7 I 13* 
1134.419 420 5+ I 13* 
1134.171 .180 5 I ha 

(1085.74) 701 10 Il 
1085.546 .540 7 Il 4 
1084.579 566 9 II 9 
1083.991 .983 7 Il 9 
(991.58) 571 10 Ill 
989.804 .803 7 Ill 9 
916.708 .698 8 II 7 
(916.01) .018 6 II 
(915.96) .963 5 Il 
915.612 .603 5 II 9 
775.966 6 Il a 
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TABLE V. Oxygen. 
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TABLE VI. Argon. 











Spectrum Wr. 














» BOWEN Ep.Lén Int. A Int. SrEcTRUM Wr 
1306.038 9 I 3 1066.660 9 I 9 
1304.864 12 I 3 1048.218 8 I 9 
1302.192 12 I 2 932.046 10 Il 9 
1217.645 10 I 10* (919.78) iS II 
999.494 6 I 9 887.404 10 Ill 9 
990.797 4 I 6 883.179 2) III 9 
990.213 8 I 4 879.622 8 It! 9 
990.121 3 I 4 878.728 il Ii! 9 
988.775 8 I 9 875.534 9 Ill 9 
(988.64) 3 I 871.099 10 Ill 9 
935.183 + I 3 - ——— —— " 
898.956 2 Ill 3 
835.293 .288 .292 7 III 18* 

(835.10) 094 .096 5 Ill : . , P : a. 
834.467 462 “462 20 il 30° tion, Shenstone has examined the extreme ultra 
833.749 741 .742 11 iil 18* violet spectrum of copper as excited in a Schiiler 
? > . = 

Rng ony ane aoe . . 30° lamp. He finds that many of the lines of Cu II 
832.762 .756 754 s II 30* in this region are those with wavelengths which 
617,060 ‘Cea TT 4 - a can be predicted with great accuracy from 
616.304 309 9 Il g* terms calculated by visible or near ultraviolet 
600.585 .583 5 Il 8* ‘ : : ca ; . ake 
$09,594 400 508 10 HI 16° combinations. It is now proposed to make 

7 i exposures with the Carnegie Institution vacuum 


580.974 .975 


| 
| 
| 
| 
i} 
| 
| 
| 
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these results be checked by another investigator 
having a grating ruled by another engine. 

When half the criteria for standards of the 
International Astronomical Union are satisfied, 
namely, satisfactory agreement between the re- 
sults of three different observers, there still re- 
mains the question of a more certain method of 
comparison with the primary standard, or at 
least with the present secondary standards. Two 
lines of approach suggest themselves to eliminate 
the difficulty of the possibility of systematic 
error due to the use of the method of overlapping 
orders. The first is most direct and involves the 
use of some form of interferometer in a vacuum. 
Experiments along this line, under the direction 
of Professor Harrison, are in progress at this 
laboratory. The second depends on the combina- 
tion principle and was suggested by Professor 
Shenstone. In an as yet unpublished investiga- 


spectrograph under conditions which will give 
from the same source these copper lines and 
some of the provisional standards suggested in 
the present article. Direct comparisons of lines in 
the same orders will then eliminate the possible 
error inherent in the present results. 

The program of research with its two-meter 
vacuum spectrograph was made possible by a 
grant to Dr. K. T. Compton from the Carnegie 
Institution of Washington. Grateful acknowledg- 
ment is also made of a grant to one of the writers 
from the Permanent Science Fund of the Ameri- 
can Academy of Arts and Sciences. We wish to 
thank Mr. D. H. Clewell for his careful assistance 
in some of the measurements. The subject of this 
paper has been under discussion for several years 
with a number of spectroscopists, all of whom are 
thanked for many helpful suggestions. It is a 
pleasure to thank particularly Dr. Compton, 
Professor G. R. Harrison and Professor A. G. 
Shenstone. 
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The Pure Rotation Spectra of NH; and ND, 


R. BowLinG Barnes, Palmer Physical Laboratory, Princeton University 
(Received March 11, 1935) 


The pure rotation spectra of NH; and ND, have been measured in the spectral region from 
40u to 170u. Ten rotation lines of NH; corresponding to final J values from 3 to 12 have been 
observed. These fall upon a curve given by »=19.880/ —0.00176/*. Thirteen lines of ND, 
corresponding to final J's from 8 to 20 have been observed. These fall upon a curve given by 


vy = 10.267 —0.00045J*. 





HE pure rotation spectrum of NH; was 

first studied by Badger and Cartwright,' 
who located five lines which fitted the formula 
v= 19.957J —0.00508/* with an average deviation 
of +0.173 cm™. Using a spectrometer which 
had a very much larger resolution, Wright and 
Randall? measured the three lines corresponding 
to J=5, 6 and 7, respectively, and made possible 
an accurate determination of the doublet width. 
This they found to be 1.33+40.03 cm™. The 
centers of the three lines fell upon a curve given 
by the formula v= 19.880/—0.00176/* with an 
average deviation of +0.02 cm™. The accuracy 
of these results is further guaranteed by their 
extremely close agreement with combination 
differences from the accurately measured near 
infrared bands. 

It was the purpose of the present investigation 
to try to extend the measurements on NH; both 
to longer and shorter wavelengths, and to make 
similar measurements for the ND, molecule. This 
has been done and the preliminary results re- 
ported in a short note.* In this present paper the 
experimental procedure and the final results are 
given. 

Since it was desirable to be able to study as 
great a portion of the far infrared as possible a 
wire grating spectrometer was selected rather 
than one employing an echelette grating. The 
instrument used was recently described by the 
author‘ and proved to be quite satisfactory from 
30u to 180. Its main features were a Welsbach 
mantle source, f-2 mirrors of 6-inch diameter, a 


*R. M. Badger and C. H. Cartwright, Phys. Rev. 33, 
692 (1929). 

*N. Wright and H. M. Randall, Phys. Rev. 44, 391 
(1933). 

*R. Bowling Barnes, W. S. Benedict and C. M. Lewis, 
Phys. Rev. 45, 347 (1934). 

*R. Bowling Barnes, R. S. I. 5, 237 (1934). 


wire grating having 25 lines to the cm so spaced 
as to eliminate the even ordered spectra, and a 
radiometer as a receiving device. In order to 
obtain pure radiation in the various regions of 
the spectrum a set of filters, such as those de- 
scribed in the paper just mentioned, was used. 
As usual the entire spectrometer was enclosed in 
a metal case, the inside of which could be dried 
with P,O; in order to remove the water vapor 
from the air. In spite of this precaution, the 
energy spectrum contained very clear indications 
of the rotation spectrum of water vapor. The 
samples of the gases studied were contained in 
glass cells of 4 cm diameter which had windows 
either of 1 mm sheets of paraffin (melting point 
68-72°C) or 0.5 mm crystal quartz. These cells 
were filled by allowing the gases to flow through 
them until a sufficient quantity of the air had 
been displaced by the gas to give the desired 
absorption. They were always used at atmos- 
pheric pressure, and ranged in length from 4 cm 
in the case of NH; to 1 cm for NDs. Transmission 
measurements were made at each wavelength 
with the cell and a dummy cell alternately in the 
path of the radiation, using a 0.5 mm plate of 
NaCl as a shutter. The deflections finally ob- 
tained were of the order of several mm and so a 
great many readings were necessary for each 
point, at least 8 being taken at each setting. The 
individual values of the percent transmission at 
any given wavelength are believed to be accurate 
to about +2 percent. The many other experi- 
mental difficulties which are encountered in this 
region of the spectrum have been discussed else- 
where, and so need not be mentioned here. 

The degree of resolution, and therefore the 
accuracy of the wavelength determinations, is 
limited in most cases by the finite width of the 
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slits which must be used in order to obtain 
readable deflections. In these measurements, the 
slits actually used were as follows: from 30 to 
85u, 0.75 mm; from 85y to 122y, 1.5 mm; from 
122u to 180u, 2.0 mm. These represent “‘slit- 
widths” at 50, 100 and 150y of 4.42, 1.67 and 
1.14 cm, respectively. In order to obtain as 
much resolution as possible readings were taken 
at intervals of } of the slit-width. From the 
absorption lines of water vapor which were ob- 
tained in the background spectrum one sees that 
minima were located which are separated by 
about 1.0 cm". 

The spectrometer was not equipped with a 
divided circle. It would, however, have been 
possible to calculate the wavelengths from the 
geometrical arrangement of the instrument, but 
as there appeared to be very slight mechanical 
defects, manifested by small shifts in the meas- 
ured position of the central image, a semi- 
empirical method of determining the calibration 
curve seemed preferable. This curve must have 
the form uw=a sin (6-Al) where Al is the differ- 
ence between the scale reading at the central 
image and at the wavelength yu, and a and bd are 
constants. In the measured energy and absorp- 
tion curves, several lines appear at long wave- 
lengths which can be accounted for only as lines 
of stated wavelengths appearing in the third 
order. For these lines 3 sin (6-Al,)= sin (b- Alyy), 
thus determining }. The three lines accurately 
measured by Wright and Randall; the neighbor- 
ing NHs lines on either side, which should lie 
closely on their formula; and four strong lines of 
H:O, whose calculated position is well known 
were now chosen as calibration points, and the 
best value of the constant a found by fitting the 
curve to these points. 

In Fig. 1 the rotation lines of NH; are shown, 
the ordinates giving the actual percent trans- 
missions, no allowance having been made for 
the slight differences of absorption or selective 
reflection of the windows of the cell and dummy. 
These curves have been left as transmission 
curves purposely, in order to show clearly just 
exactly what has been measured. Particular at- 
tention is directed to the regions between the ab- 
sorption lines, since they give a very good indica- 
tion as to the spectral purity of the radiation 
used. Table I, which gives a summary of the re- 
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Final Obs. Obs. Calc. Obs.? Obs." 


J as ¥ y dl 


12 42.47 235.40 235.48 
11 46.15 216.61 216.31 
10 50.76 196.93 197.02 








9 56.26 177.69 177.622 176.10 

8 63.14 158.33 158.129 156.8 

7 72.18 138.50 138.549 138.54 

6 84.07 118.92 118.896 118.90 118.60 

5 100.78 99.20 99.178 99.180 99.06 

4 126.04 79.32 79.406 79.79 
10/3 152.35 65.62 65.67 

3 


168.42 59.36 59.592 





sults available on the far infrared of the NH; 
molecule, is self-explanatory. The calculated 
values of vy were obtained by using the formula of 
Wright and Randall given above. The observed 
values agree with the calculated v's to +0.11 
cm. As may be seen from the curves lines 3 to 6 
show definite signs of being double. Line 3 falls 
in a region where part of the energy comes from 
overlapping third orders, and so the exact deter- 
mination of its wavelength is difficult. The line 
at 152.35 is the third order of line number 10 
which occurs at 50.76,. 

Samples of heavy ammonia, NDs, were fur- 
nished us by Professor H. S. Taylor. These were 
prepared by reacting 99+ percent pure heavy 
water with magnesium nitride.’ The resulting 
ammonia contained 97 percent D, as shown by 
ultraviolet absorption measurements. Fig. 2 
shows again the actual transmission measure- 
ments. The simplicity of the spectrum found indi- 
cates that the absorbing gas in our cell was mostly 


TABLE II. The far infrared of ND3,. 


Final Obs. Obs. Calc. 
J yu v v 
20 49.46 202.10 201.60 
19 51.85 192.80 191.85 
1s 54.93 181.99 182.05 
17 58.09 172.10 172.20 
16 61.66 162.12 162.31 
15 65.76 152.02 152.38 
14 70.31 142.19 142.41 
13 75.13 133.07 132.42 
12 80.53 124.14 122.35 
11 88.40 113.09 112.26 
10 97.70 102.32 102.15 
9 108.89 91.81 92.01 


8 122.60 81.54 81.85 


>H.S. Taylor and J. C. Jungers, J. Am. Chem. Soc. 55, 
5057 (1933). . 
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Fic. 1. Pure rotation spectrum of NHs, showing absorption lines corresponding to final J values from 3 to 12. The line 
at 152y is a third order line. (The last numeral in each of the wavelength numbers is a decimal and should be preceded 


by a decimal point.) 
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Fic. 2. Pure rotation spectrum of ND,, showing absorption lines corresponding to final J values from 8 to 20. (The last 
numeral in each of the wavelength numbers is a decimal and should be preceded by a decimal point.) 


NDy. A general background absorption (absent at all surprising. Table II gives a summary of the 
in the NH; curves), however, does suggest the NDs results. The calculated frequencies are given 
presence of small percentages of unsymmetrical by the formula v= 10.26/—0.00045/J*. It is seen 
molecules such as ND,H and NDHz,:. This is that the agreement with observation, particu- 
confirmed by ultraviolet measurements on the larly in lines 11-13, is not good. This is undoubt- 
sample as used in our cell which indicated a_ edly due to the disturbing effect of the NHD, 
purity of 85 percent D. Considering the crude absorption, the background being particularly 
method of filling the cells made necessary by the _ strong in this region. However, the question must 
thin windows used, the presence of these impuri- be raised how accurately the simple formula, 
ties, due to exchange with atmospheric H,O, isnot »=2BJ—4DJ* derived for linear molecules, is 
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Fic. 3. Schematic representation of the pure rotation 
spectra of NH; and ND). Intensities as shown are only 
approximate. 
applicable to the symmetric rotator type. Dr. W. 
S. Benedict has shown that additional terms in 
the cubic formula are necessary, and are particu- 
larly important in NDs. Since the present experi- 
mental results are not accurate enough to test 
this point, details will be reserved for a later time. 
For the present, the simple formula is taken, with 
Dyp,= }Dwu,; the resulting value of Byp,= 5.13. 
From this and Byy,=9.94, the height of the 
pyramid may be calculated, since 1/B~IJ and 

3mMyMy 


Ixu, = (3/2)mya? +————¢? 
My +3my 


to an H or D atom, and q is the height. Solving, 
we find g=0.36A, ry_qg=a?+¢@=1.02A. The ac- 
curacy of this result is not high, but the agree- 
ment with the value calculated indirectly, 0.37A, 
is good.® 

Fig. 3 shows schematically the relative posi- 
tions of the rotation lines of the two symmetrical 
ammonias. 

The author wishes to express his appreciation 
to Drs. C. M. Lewis and W. S. Benedict who 
aided greatly in making these measurements, and 
to Professor Taylor and Dr. J. C. Jungers who 
prepared the heavy ammonia. 





* Manning, J. Chem. Phys. 3, 136 (1935). 
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On Isotopic Shifts in the Spectra of Diatomic Molecules 


G. H. Drexe, Johns Hopkins University 
(Received March 8, 1935) 


An inspection of the theory of diatomic molecules shows 
that the elementary theory of the dependence of the 
energy of diatomic molecules on the mass of the nuclei 
must be supplemented in two ways. In the first place the 
reaction of the nuclei to the precession of the electronic 
angular momentum about the internuclear axis gives a 
term 


[B,—a( V+ }) JL L(L+1) —A?]. 


Furthermore there is a second positive term of uncertain 


HE study of the spectra of molecules in 

which one or more of their atoms are sub- 
stituted by the corresponding isotope can often 
furnish information about the structure of the 
molecule which cannot at all, or only with 
difficulty, be obtained in any other way. This is 
especially true for those molecules in which a 
hydrogen atom is replaced by its heavy isotope 
deuterium, because for this substitution the 
relative change in the masses is so much greater 
than for all other isotopic substitutions. In order 
to make full use of the material which can be 
obtained in this way it is necessary to know 


magnitude. Both these terms are of importance only for 
very light molecules especially for the hydrogen molecule. 
The large electronic shifts observed for some band systems 
of H:, HD and D, are well accounted for by this theory. 
A very accurate comparison of the constants of very light 
isotopic molecules necessitates also taking into account 
small corrections which appear in the higher order ap- 
proximations of the energy values for a rotating anhar- 
monic oscillator. 


exactly the changes which must be expected in 
the structure of the corresponding energy levels 
and wave functions. The approximate theory 
which was hitherto quite adequate to give a 
satisfactory account of all the changes as long as 
only heavy atoms were concerned is not able to 
explain all the details for proton deuteron substi- 
tutions. Kronig' called attention to the fact that 
if certain terms in the Hamiltonian of a diatomic 
molecule which are usually omitted are taken 
into account, there is for many states a slight 
difference in the equilibrium distance even for the 


“e R. de L. Kronig, Physica 1, 617 (1934). . 
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nonvibrating and nonrotating molecule. It is 
very difficult to prove such a difference experi- 
mentally, and it is very doubtful whether the 
effect discussed by. Kronig is the cause for the 
discrepancies found by Holst and Hulthén® in 
the spectrum of AIH and AID. 

An investigation of the spectra of He, HD and 
D, must be expected to give the most detailed 
empirical information about the deviations from 
the elementary theory, as these deviations are all 
inversely proportional to the mass and therefore 
especially large for these light molecules. It is 
found that there are large electronic isotope shifts 
up to 25.6 cm™. The electronic isotope effect 
usually is understood to be the shift in the elec- 
tronic frequency of a band system when one iso- 
tope is substituted for another. The electronic 
frequency is the frequency which the nonrotat- 
ing and nonvibrating molecule would emit. It 
does not correspond to an actual line but can be 
obtained by extrapolation from the observed 
band lines of the system. In the following it is 
shown that shifts of this magnitude are de- 
manded by theory. Furthermore, it is shown 
which terms have to be considered in order to 
find out a more accurate expression for the 
dependence on the molecular mass of the vibra- 
tional frequency and other constants. 


$1. 
We shall restrict ourselves to the nonrotating 
molecule. The vibrational energy is then deter- 
mined by 





[ nr) — van) + We lRe) =0, (1) 
8r2ur* ar ” he 


in which R(r) is the vibrational eigenfunction 
and V,(r) the potential energy for fixed nuclei. 
If the electronic wave function for fixed nuclei 
is , the total wave function for the nonrotating 


molecule is 
Y=-R. (2) 


(1) and (2) are valid only approximately, as cer- 
tain terms of the complete Hamiltonian have 
been neglected.’ We can introduce these terms 


?W. Holst and E. Hulthén, Nature 133, 496 (1934); 


Zeits. f. Physik 90, 712 (1934). 

3J. H. Van Vieck, Phys. Rev. 33, 467 (1929). See also 
R. de L. Kronig, Band Spectra and Molecular Structure, 
Cambridge, 1930. 
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now-as perturbations with the perturbation func- 
tion H’ such that 
h? 0° 
H'y= (M?+M,’)¢—r— |i (3) 


8r2ur? or? 


M:, M,, M; are the components of the orbital 
angular momentum‘ with respect to axes fixed 
to the molecule. (The ¢ axis is the internuclear 
axis.) We shall consider in this paper only the 
first term in this expression. For an evaluation of 
the influence of the second term a knowledge of 
the ® would be necessary. The empirical data 
seem to show that its influence for the two levels 
for which sufficiently accurate data are available 
is small compared to that of the first term. 

The influence of (3) on the energy is in the first 
approximation 


W,'= f VII’ pdr = f RRadr f @Il'bdv’. (4) 


We have 
(M;?+M,")*=(M’?—M,;’) 
=[L(L+1)—A*?]}#, (5) 


in which M is the operator corresponding to the 
total orbital momentum. L is the quantum num- 
ber corresponding to this operator and A that 
corresponding to Mf. (5) is correct only if the 
angular momentum is constant. This will be the 
case especially for higher levels of light mole- 
cules, but (5) can be expected to be a fairly good 
approximation also for most other molecular 
levels. The contribution of the first term in (3) 
to the energy is therefore, if we take for r its equi- 
librium value, B..L(1+1)—A?*]. To obtain a 
better approximation, we remember that (4) 
means averaging over the electronic motion and 
the nuclear vibrations. The averaging over the 
electronic motion is taken care of by (5), whereas 
the averaging over the nuclear vibrations con- 
cerns only the quantity B=//8x*yr*. This same 
average occurs when the rotational energy of a 
rotating oscillator is calculated and its value is 


: We restrict ourselves then to singlet levels, as in general 
M is the éofal angular momentum. We can nevertheless 
ply the results obtained in this way to the triplet bands 
are drogen, as in this case the spin is so loosely coupled 
to the rest of the molecule that for all practical purposes 
the properties of the triplet levels are the same as those of 
singlet levels. 
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ISOTOPIC SHIFTS IN 


B.—a(V+4)+---. It is essential for this thet 
the deviations from harmonic oscillations are 
taken into account, for otherwise a would have 
the same order of magnitude but opposite sign. 
The contribution of the first term in (3) to the 
energy is thus 


W,' =(B.—a(V+4) ][.L(L+1)—-a?]. (6) 


The part independent of the vibrational quantum 
number V would in the usual empirical analysis 
of the energy levels of a diatomic molecule be 
classed with the electronic energy. It differs how- 
ever from it because it depends on the mass of the 
nuclei and does not occur in the energy of the 
molecule with fixed nuclei. Its physical signifi- 
cance is that it represents the reaction of the 
nuclei to the precession of the electronic angular 
momentum about the internuclear axis. (3) 
shows that it is proportional to the square of the 
component of the angular momentum perpen- 
dicular to the internuclear axis, and it is evident 
that this must be present also for the nonvibrat- 
ing and nonrotating molecule. The magnitude of 
this term can be easily calculated as all the con- 
stants are known from the rotational and 
vibrational analysis of the corresponding bands. 
The term is zero only for SZ levels. As this and 
all other similar terms are proportional to B it is 
evident that they are most significant for the 
hydrogen molecule because on account of the 
light masses of the hydrogen nuclei B is so much 
larger than for any other molecule. 


§2. 
COMPARISON WITH EXPERIMENTAL VALUES 


It is obvious that the constant part of (6) can- 
not be separated from the true electronic energy 
if only the H, levels were known, but if also the 
corresponding HD or D, levels are known the 
two parts can be determined separately. There 
are in the spectrum of the hydrogen molecule two 
levels which are analyzed with sufficient detail 
and accuracy to allow an experimental test of the 
above formula. These are the 3p*II level® for 
which L=1 and A=1 and the 3p*> level® for 
which L=1 and A=0. The observed origins of 


°G. H. Dieke and R. W. Blue, Phys. Rev. 47, 261 (1935). 
*G. H. Dieke, to appear shortly in Phys. Rev. 
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the band systems corresponding to transitions 
from these two levels to 2s*= are: 


H; HD D: 
3p*Il—-2s*= 16,793.87 16,788.86 16,783.92 
3p'E —2s8*Z 11,838.91 11,825.68 11,813.40 


The observed shifts and those caiculated from 
(6) are: 


H,—HD HD—D, 
obs. calc, obs. calc. 
3p'Tl—+2s*= 5.01 7.55 4.94 7.61 
3p'M—2s8E 13.23 13.44 12.28 12.62 


Lt 
It is difficult to judge the accuracy of the ob- 
served origins. They must be obtained by extra- 
polation from many observed lines and there is 
always the danger of systematic errors of uncer- 
tain magnitude if the formula used is not a suffi- 
ciently good approximation. The uncertainty 
due to errors in the measurements is only a few 
hundredths of a wave number and it is estimated 
that the total uncertainty is not more than one 
wave number, probably much less. The only ex- 
ception is perhaps the value for the origin of 
3p*2—2s*S of Hz for which Richardson's value 
was taken which was computed from less complete 
data, and may have therefore possibly a larger 
error. 

It is seen that the general agreement between 
the empirical and theoretical values is very satis- 
factory, which shows that very probably the 
largest part of the observed shifts is given by 
(6). In particular the isotopic shifts have the 
right sign and the shift for the p*= levels is 
about twice that of the p*II level as it should be. 

As a rule there is no difficulty in finding the A 
value for a given electronic level from the struc- 
ture of the bands in question. To find the L 
values is much more difficult and requires a much 
more complete knowledge of the whole spectrum. 
With the isotope effect discussed here it may be 
often possible to fix the L-values in cases of 
doubt. 

The remaining discrepancies between the ob- 
served and calculated values may be partly due 
to the neglected second term in (3). This term 
must also be positive, but there is no reason why 
it should be zero for s= states. It must be ex- 
pected that its contribution will be larger for the 
lower states, so that the difference between the 
contribution of the initial and final state which 
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enters into the band origins will be negative. The 
corrections due to the second term in (3) have 
therefore the right sign to reduce the discrepan- 
cies between the above theoretical and observed 
values. 

There are two other band systems of hydrogen 
for which at present there are sufficient data for a 
calculation of the electronic shifts. They are the 
ultraviolet band systems 2p'II—1s'E’ and 2p'D 
—1s'd.* The difference between the H, and HD 
origins of both systems has the opposite sign 
from that of the two visible systems. The values 
of the shifts are 8 cm™ for 2p'II—+1s'= and 12 
cm for 2p'II-+1s'E with a probable error of 
about 5 to 10 cm. If these values are correct it 
shows that the second term in (3) must be much 
larger for the normal state 1s'S than for any 
excited state, so that it can compensate for the 
positive shifts due to the first term in (3) which 
are of the same order of magnitude as those 
found for the visible bands. Such an assumption 
does not seem altogether improbable. 

For the lines of the hydrogen and deuterium 
atoms we find a shift due to the common motion 
of nucleus and electron about the center of grav- 
ity. It is 8.25 cm™ for Ha and 44.5 cm™ for La, 
the line of the heavier atom having the higher 
frequency. Such an effect, because of the fact 
that the center of gravity of the nuclei and elec- 
trons do not coincide, is not taken into account 
in Van Vleck’s treatment of diatomic molecules 
on which the foregoing considerations are based. 
It must be expected, however, to be much smaller 
for systems with more than one electron. The 
motion which the center of gravity of the nuclei 
executes because of the motion of one electron is 
at least partly counterbalanced by the motion 
of the second electron. Therefore the total effect 
will be much smaller. This is analogous to the 
fact that the energy level series of the helium 
atom are best represented with the Rydberg 
constant for an infinite mass and not with that 
for a mass four. 


§3. 
Tue EFFECT ON THE VIBRATIONAL FREQUENCY 


The term linear in +4 of (6) has the effect 
that the frequency w, calculated according to 


ICR. je Je n, Phys. Rev. 45, 480 (1934). 
*K. Mie, eits. f. Physik 91, 475 (1934). 
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(2) with V,(r) as potential must be diminished 
by af L(L+1)—A?] in order to obtain the total 
coefficient w,’ of V+ 4 in the energy formula. w,’ 
is the quantity which is determined empirically 
and usually designated as the vibrational fre- 
quency. 

In going from one isotope to another w, varies 
with the first power of the transformation coeffi- 
cient p, whereas @ varies with the third power.’ 
w,’ which is the quantity empirically determined 
from the experimental data varies, therefore, not 
exactly with p but should show slight variations. 
A rigorous experimental proof of this is at present 
not yet quite possible as the empirical values of 
the constants are not quite reliable enough. 





Table I shows the effect of this correction. This 
TABLE I. 

HD: H H, D,: H: D.: HD 

3p*l {without corr. 0.86634 0.70764 0.81682 

with a@ corr. 0.86620 0.70741 0.81668 

3p°= = { without corr. 0.86751 0.70881 0.81706 

(with a@ corr. 0.86722 0.70832 0.81676 

calc. 0.70744 0.81675 


0.86616 





table is given with the intention to show how 
much these corrections affect the w,’ ratios 
rather than to show the agreement with the calcu- 
lated values of p of the last row. The corrections 
of the next paragraph may be of the same order 
of magnitude, and it would be necessary to take 
them into account also if the comparison is in- 
tended to be of such a high degree of accuracy. 
It may be remarked, however, that the agree- 
ment between the observed and calculated values 
for 3p*II is of nearly the same order with and 
without the corrections. The 3p*> state offers a 
better case as the corrections are about twice as 
big as those for the 3p*II state. The agreement in 
the first two columns is poor in any case. This is 
very probably because of a too low w,’ value for 
Hy, which was calculated by Richardson and 
Das” from much scantier and less accurate em- 
pirical material. The last column shows, however, 
a much better agreement with the correction 
than without it. 


~ * This statement is subject to the corrections of § 4. 
OQ. W. Richardson and K. Das, Proc. Roy. Soc. A125, 
309 (1929). 
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§4. 


CORRECTIONS FOR THE ENERGY VALUES OF 
THE ANHARMONIC OSCILLATOR 


The constants in the familiar expression, 
we( V+4) —x(V+4)?+9(V+9)?—2(V+49)*+--- 


for the vibrational energy of an anharmonic 
oscillator are according to the elementary theory 
multiplied by the first, second, third and fourth 
power of p, respectively, in an isotopic substitu- 
tion. This is true, however, only approximately 
and if one goes to higher approximations, cor- 
rection terms have to be added which require a 
different power of p as isotopic transformation 
factors. These corrections do not show in the 
formulae which Fues" first derived with the 
wave mechanics because he did not go to high 
enough approximations, but they are fully given 
in a treatment by Dunham.” If the potential 
energy is written 


V(r) =a0d*{ 1+41:5+a26*+ --- ] 


in which 6= (r—79)/1; all the constants are sup- 
posed to be independent of the nuclesr masses. 
The energy of a rotating oscillator having this 
potential can be written in the form 


E ¥n(V+)) (I+) 
. @ 


and Dunham gives the first 15 of the coefficients 
Y,. A few examples are sufficient to show what 
kind of corrections appear. If 8= B,/w,, we have 


Yoo = (B./8)(3a2—7a;?/4), 
Vio = w{ 1+ (87/4) (25a,4—95a3a3/2 

— 67a_*+459a,7a2/8 — 1155a,*/64) ], 
Yo: =B[1+(8?/2)(15+14a,—9as 


+ 15a; —23a,a2+21(a,?+<a;')/2], etc. 





" F, Fues, Ann. d. Physik 80, 367 (1926). 
'? J. L. Dunham, Phys. Rev. 41, 721 (1932). 


It is not the purpose of the present paper to 
give an accurate evaluation of these corrections,"* 
but a few remarks may show of what order of 


magnitude they are. 

Let us consider first the term Yo which can be 
considered as a correction to the electronic 
energy. It can be evaluated with ease if we make 
use of the approximate expressions 


Yo = —x=(3B,/2)(a2—5a,7/4), 
Yi; = —a=6B8(1+a,). 


It is sufficient to give the constants for Ds. 
They are 


sg a; a: Yoo 
2s*= 110.4 —1.59 1.83 2.24 
3p" 110.4 — 1:67 2.03 2.34 
3pz 113.1 —1.70 1.95 1.34 


For HD Yq is 1.5 times and for Hy, twice as big 
as for Dz, which means that the shift between H, 
and D, is just equal to the value of Yo for D2. given 
in the table. The isotopic shift of the origin of 
the band systems which is the difference between 
the shifts of initial and final level is therefore 
+0.10 for 3p*II--2s*= and —0.90 for 3p*= 
—2s*Y. These shifts, though of sufficient magni- 
tude to be noticeable with accurate observations, 
are small compared with the effect discussed in 
§2, and will be in general altogether insignificant 
for heavier molecules. 

The effect of these corrections on w, and B, 
may be of sufficient magnitude to change ap- 
preciably the empirical ratios of these constants 
for the different hydrogen isotopes and should 
therefore be taken into account, if an accurate 
test is made. This will be considered in a later 


paper. 


For a reliable calculation of the constants a), as, a: 
and a, which are all necessary for an evaluation of the cor- 
rection for w, it would be necessary to know accurately 
some of the constants which occur only as small terms in 
the empirical expression of the energy. It would be a waste 
of time to attempt this unless the best available data can 
be used. As we are still improving constantly the empirical 
data for the three hydrogen molecules, we are postponing 
a detailed discussion of all these corrections to a later date. 





TRY ae SO Si 





MAY 1, 1935 


PHYSICAL REVIEW 





VOLUME 47 


A Modified Aston-Type Mass Spectrometer and Some Preliminary Results 


D. D. TayLor, Norman Bridge Laboratory of Physics, California Institute of Technology 
(Received January 7, 1935) 


A modified Aston-type mass spectrometer is described, in which ionization is produced by 
controlled electron bombardment. The electric field used is radial, and the shape of the 
magnetic pole faces was calculated to give the necessary velocity focusing. A special type 
of filament assembly is employed to minimize the energy spread of the emitted electrons. An 
amplifier circuit is described using a Western Electric D-96475 tube which gives high sensitivity 
combined with good stability. A new device is desc-ibed for giving a continuous indication 
proportional to the square of the magnetic field. Some preliminary results are included on the 


ionization of Ns, CO, NHs, N,H,g. 


N order to obtain increased resolving power for 

the study of positive ions resulting from con- 
trolled electron bombardment, a modified Aston- 
type mass spectrometer has been constructed 
and some preliminary results will be given. The 
general arrangement is shown in Fig. 1. The 
gas to be investigated is admitted at G into a 
capillary tube. The gas is released from the 
capillary at about the center of the filament 
assembly F where it is bombarded by electrons 
accelerated radially inward. The resulting ions 
are then accelerated by a 1500-volt potential 
applied between F and the anticathode*A. A is 
grounded through a microammeter to permit 
measurement of the positive ion current inci- 
dent on it, which is used as an indication of the 
total current through the mass spectrometer. 
The collimator consists of two slits, C, 40 cm 
apart and each 2.5mm X0.1 mm. The collimated 
beam next passes through the electric deflecting 
field E whose plates are equiangular sectors of 
coaxial cylinders, and in this field the beam is 
spread out into an energy spectrum. This di- 
vergent beam next enters the magnetic field M 
whose leading and trailing edges were calculated 
to focus all particles of a given value of e/m on 
the focal slit S, and any desired value of e/m 
can be focused by suitable adjustment of the 
magnetic field strength. The leading and trailing 
edges of M are slightly concave but no attempt 
has been made to indicate their true shape for the 
drawing is too small. 

With the exception of necessary electrical 
insulation the apparatus is constructed of metal. 
The magnetic circuit is made of Rema dead soft 
electric furnace charcoal iron. The remainder of 


the apparatus is of brass, with gold slit jaws and 
gold plated deflecting plates. The plug of the 
large ground joint carrying the filament as- 
sembly is of Pyrex ground into a steel sleeve. 

The path of a charged particle in the radial 
electric field was kindly obtained by Mr. G. C. 
Munro of this laboratory. Neglecting edge 
effects the field between the two plates is an 
inverse first power field. All the rays are as- 
sumed to be parallel and to enter the field at the 
same point, on the leading edge midway between 
the plates. It is required to find the emergent 
paths and velocities for the velocity determines 
the radius of curvature in the magnetic field and 
the emergent path determines the point of 
incidence in the field. The point of emergence of 
the ray on the trailing edge of the electric field, 
the emergent velocity, and the angle between 
the emergent ray and a reference ray were each 
determined in terms of one of three convergent 
infinite series in a parameter, x= $0/2E In (p2/p:). 
¢@» is the voltage applied between the electric 
field plates of radius pz and p;, and E£ is the 
energy of the particle in electron volts. If the 
center of curvature of the plates is taken as the 
origin, a circle of radius a about this point 
passes through the point of incidence of the 
rays and the radius to the point of emergence 
of any ray is p=wua. u is described by u=1 
+(1—x«)f(x). 

The tangent of the angle between the emergent 
ray and the reference ray is (1—«)F(x). The 
relation between the emergent velocity V and the 
incident velocity v is V?=27[ 1 —«(1—x«) d(x) ]. 

A range of «x from 0.7 to 1.3 covered the entire 
usable part of the field but the calculations 
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Fic. 1. Diagram of mass spectrometer. G, gas inlet to capillary; F, filament assembly—see Fig. 2; 
A, graphite anticathode; C, collimator slits; E, electric deflecting field; M, magnetic deflecting held; 
S, focal slit; 7, Western Electric Co. D-96475 tube; V, vacuum pumping connections. 


were extended from «=0.5 to 1.5. Since the 
analysis of the electric field disclosed no analytic 
expression for the characteristics of the beam 
after traversing the electric field, the cnly 
available method of determining the proper 
shape of the magnetic field lay in assigning 
appropriate values to the parameter « and then 
determining radius of curvature and point of 
incidence on the magnetic field for a sufficient 
number of rays. It was found that the calcula- 
tion could be greatly simplified by introducing 
the further restriction on the magnetic field that 
it be symmetrical with respect to the bisector of 
the angle between the incident and emergent 
median rays. A total of fifty-five points was 
calculated on the leading edge of the magnetic 
field and then by means of a Schwartz transfor- 
mation the edge effect was investigated and the 
points were corrected accordingly. The results 
of the edge effect calculation show that to a 
reasonable approximation the line integral of the 
magnetic field along the path of a particle is the 
same as if the uniform field were taken to extend 
beyond the edge a distance equal to eight- 
sevenths of the air gap. 

An attempt was made to approximate the 
curve by means of an ellipse which could be 
generated on a milling machine, but no conic 
section gave a satisfactory fit. 

The half-width of the iron at the bottom of the 
field is 1.789 inches, the minimum half-width 
is 1.768 at 1.600 inches above the bottom, and 
the half-width is 1.808 inches at a point 3.0 
inches from the bottom, above which no rays 
enter. 


A sylphon is introduced between the mag- 
netic field and the focal slit to permit radial and 
tangential motion of the focal slit for adjustment 
of focus. 

Unfortunately the theoretical work of R. 
Herzog! was not available when this apparatus 
was designed but it has since been applied. The 
electric deflecting field has a radius of curva- 
ture for the median ray of 12 cm and is equivalent 
to a cylindrical lens of 24.5 cm focal length plus 
a prism. The magnetic field has a radius of 
curvature of 10 cm for the median ray and is 
equivalent to a cylindrical lens of 11.9 cm focal 
length plus a prism. In both fields the inter- 
section of the principal planes coincides with the 
intersection of the median rays. The first collima- 
tor slit has an inverted image 33.3 cm in front 
of the focal slit S, and the image of the second 
collimator slit is 23.8 cm in front of S. The 
width of image observed at S is consistent with 
that calculated from the positions and dimen- 
sions of their optically determined images. 

Vacuum pump connections are provided at 
the three points marked V. 

The filament assembly is shown in more 
detail in Fig. 2, and has proven quite satis- 
factory in most respects. The entire body is 
made of Acheson graphite and consists of two 
disks and a short cylinder attached to one of 
them. Eight helical tungsten filaments F are 
arranged as shown and no electrical connection 
to them is required other than a sliding fit in 
the holes drilled in the graphite. The grid G is 


'R. Herzog, Zeits. f. Physik 89, 457 (1934).« 
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Fic. 2. Filament assembly diagram. G, helix of tungsten 
wire for grid; F, helical tungsten filaments. 


also helical and is so proportioned that it with- 
stands a potential drop equal to and in the same 
direction as that in the filaments, thus largely 
eliminating the variation in electron energy due 
to the filament drop. This also results in an axial 
field tending to draw the positive ions out of the 
filament assembly where they can be affected by 
the larger field applied between F and A. 
The chief disadvantage of this arrangement is 
the penetration of the larger field into the fila- 
ment assembly so that as yet no check on ioniza- 
tion potentials is possible. 

The amplifier circuit employed, as shown in 
Fig. 3, is notable chiefly for its simplicity for it 
has not been found necessary to modify the 
circuit for individual tube differences. The re- 
sistance values given were calculated only to 
give the recommended potentials and the be- 
havior of the circuit has been most satisfactory 
from the start. The zero drift is consistently 
less than 2 cm per hour with a sensitivity of 10-"” 
amp./mm, and random fluctuation of 3-4 mm, 
which compares favorably with the various 
compensated circuits published up to this time. 
Two factors contribute to this stability. The 
first is the use of large batteries, 240 amp. hr. 
capacity, so that the voltage drift is negligible. 
The second is the careful shielding against 
thermal disturbances in addition to the neces- 
sary careful electrical shielding. All of the 
resistances in the circuit, excepting only the 
grid resistance of 10" ohms, are immersed in 
transformer oil which effectively prevents tem- 
perature fluctuations. The number of variable 
contacts is reduced by soldering all possible 
contacts, leaving only those needed for adjusting 
filament current and galvanometer zero. Some 
trouble was caused by poor insulation resistance 
in the leads from the control box to the tube 





TAYLOR 











oO 
F 
A ~ MV 
YWWVYLYV VAG 
BAMIVYV 


I 
MAA] oH 
| 3 

J 


Fic. 3. Circuit diagram. Resistance AO=1.3 x 10"; 
AB=11.1; CD=12.25; DE=8.65; EF=25,000; CG =200; 
EH=200 (dual); JJ/=2; KL=10; KN=200; LM=40. 
Battery—i2v, 240 amp. hours. Galvanometer Type R, 
sens. about 10~*.amp./mm. 
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but the substitution of lead covered telephone 
cable eliminated the leakage. (This cable can be 
obtained having a guaranteed minimum insula- 
tion resistance of 500 megohm miles.) 

The tube employed is the Western Electric 
D-96475 and it is a very definite improvement 
over the similar tubes previously available. It is 
mounted, as shown at 7 in Fig. 1, directly in the 
high vacuum of the mass spectrometer, and is 
protected by a half-inch of Rema iron from the 
stray field of the magnet M. In this construction 
the tubulation for the control grid lead is most 
convenient. The high resistance is mounted on 
the grid cap along with a small collecting elec- 
trode immediately behind the focal slit. 

An auxiliary device is employed which gives a 
continuous indication of magnetic field strength. 
This takes the form of a plate of diamagnetic 
material, in this case Acheson graphite, sus- 
pended in the upper part of the magnetic field 
from a torsion balance. The force exerted on the 
plate tending to push it out of the field is pro- 
portional to the square of the field strength, and 
thus directly proportional to the m/e value of 
the ion brought to focus by that field strength. 
A mirror is mounted on the torsion balance for 
optical reading and the molecular weights so 
obtained readily fall within one percent of 
the true values, which is amply close for iden- 
tification. 
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An attempt to utilize piano steel wire for the 
torsion fiber resulted in excessive zero shift due 
to elastic hysteresis, but the use of pure tungsten 
wire, even without annealing, has made ‘this 
effect inappreciable. 

The taking of observations involves plotting 
simultaneously the magnetic balance deflection 
and the positive ion current to the amplifier 
against an independent parameter, which for 
convenience is taken as the resistance inserted 
in the magnet circuit. This resistance is con- 
tinuously variable from its maximum value of 
twenty-two ohms to about one-tenth ohm. 

The results are presented in tabular form 
with one specimen curve from each substance 
to give a better idea of the form of the results 
obtained. To facilitate comparisons the peak 
heights are all given as percentages of the height 
of the primary peak. The column headed ya 
gives the positive ion current received on the 
anticathode, the next column shows the electron 
current received on the grid, and the volts 
column shows the voltage impressed between 
the filament and the grid. However, because of 
interpenetration of the accelerating field for 
the ion, this is in all cases higher than the true 
value of the energy possessed by the bombarding 
electrons. The next column shows the true 
height in cm of the primary peak, one cm repre- 
senting a positive ion current to the collector of 
approximately 10-" amp. The pressure in the 
ionizing chamber was of the order of 10-* mm. 


NITROGEN 


The nitrogen used was prepared by the thermal 
decomposition of sodium azide NaN; which is a 
very convenient source of pure nitrogen. Nitro- 
gen within the voltage range considered yields 
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Fic. 4. Run No. 1. Nitrogen. 118 volt electrons. Height of 
N, peak = 705 cm. 


predominantly molecular ions and shows satis- 
factory agreement with previous work on this 
substance. Although the source of the nitrogen 
precludes any possibility of its contamination 
with water vapor, nevertheless the walls of the 
apparatus apparently yield some water vapor 
and some organic vapor may be present from 
stopcock grease to account for the hydrogen 
found. The carbon peak is satisfactorily ac- 
counted for in the use of a graphite mounting for 
the filaments. No satisfactory explanation has 
been found for the two peaks appearing con- 
sistently at 39 and 41. The NH; and NH, peaks 
are certainly enhanced somewhat by OH and 
H,O but no separation of the two was possible. 
The peak corresponding to atomic hydrogen 
ion was beyond the range of the apparatus so was 
necessarily omitted. Run No. 1 is shown as 


Fig. 4. (See Table I.) 


TABLE I. Jon intensities observed for nitrogen bombarded with electrons of various energies. Intensities calculated to basis 


N,* = 100. wa = positive ion current incident on anticathode; ma = electron current incident on grid. 











True 
Run ya ma _ @=~Volts N; H; * N NH NH, NH, NH, CN N; 39 41 
1 12 3 118 705 5.5 0.31 14 0.43 0.57 0.50 1.0 0.17 1009 0.17 0.26 
2 2 100 270 8.9 56 18 56 74 .67 1.1 37 100 37 74 
3 11 3 80 110 8.6 27 29 45 1.8 1.8 3.6 -- 100 - 1.8 
4 i) 3 60 86 7.0 1.7 17 46 1.2 1.2 2.1 - 100 — 
5 5 7 40 56 13. - 2.1 1.1 9 2.0 5.4 100 — 1.1 
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TABLE II. Jon intensities observed for carbon monoxide bombarded with electrons of various energies. Intensities calculated to 
basis CO* = 100. 




















True 
Run 4a ma Volts CO H, fe CH CH, CHs O OH H,O CO 39 41 CO; 
6 9 4 120 2075 15 1.7 0.05 0.38 0.29 0.77 0.39 1.2 100 0.10 0.14 0.63 
7. 5 120 1175 xX 3.6 .04 43 38 1.3 47 1.5 100 17 17 .60 
8 25 12 120 430 xX 19 46 1.2 58 5.8 .93 2.3 100 35 58 .93 
9 50 26 120 14 xX 370 7.0 18 7.0 160 100 
10 20 12 120 1560 4.5 5.1 13 32 22 1.2 38 7 100 13 .06 
11 11 1.5 120 130 4.2 18 - .6 2.3 3.8 20 100 a7 54 
12 8 3 100 365 6.8 9.8 -— 27 .02 3.2 .6 » 100 55 44 $1 
13 11.5 18 80 29 52 21 - - 14 3.5 6.9 100 2.8 
14 9 5 80 125 4.4 19 38 56 48 6.4 1.6 100 
15 10 6 60 60 5.3 15 - — _ 3.3 2.5 100 
3.7 5 2 1.5 2.5 100 


16 5 10 40 40 


X =not covered inrun. — =not found. 


TABLE III. Jon intensities observed for ammonia bombarded with electrons of various energies. Intensities calculated to basis 











NH;,* = 100. 
True 
Run pa ma Volts NH; H, a N NH NH, NH; NH, Ne 
17 10 0.1 117 18 14 33 78 8.3 86 100 139 680 
18 20 2.5 117 7 43 36 190 21 130 100 71 400 
19 9 35 114 17.5 31 7.1 7 10 104 100 74 40 
20 19 8.5 110 21 76 16 12 12 107 100 100 79 
21 30 11 110 17 47 7.3 13 13 130 100 12 35 
22 9 = 97 15 13 33 90 11 93 100 103 640 
23 7 0.1 81 9.7 - 26 103 10.3 8 100 90 850 
24 9 3 60 5.7 70 13 17 20 160 100 26 53 
25 9 2 60 3.3 75 45 120 - 100 100 100 710 
26 8.5 3.5 39 3.1 19 . . - 105 100 32 48 
27 7.5 7.5 40 1.7 60 60 — 60 100 94 590 





CARBON MONOXIDE 


The carbon monoxide used was prepared by 
the dehydration with sulphuric acid of formic 
acid. 

An interesting and unexplained 
shown in runs 7, 8 and 9 which were made con- 
secutively with everything unchanged except 


effect is 
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Fic. 5. Run No. 7. Carbon monoxide. 120 volt electrons. 
Height of CO peak = 1175 cm. 


the 
anticathode. It 


the electron current and correspondingly 
the 
seems strange that a fourfold increase in the 


positive ion current to 
positive ion current should result in a decrease 
of the primary peak height to scarcely more than 
one percent of its initial value. It is presumably 
this abnormal depression of the CO peak which 
causes the apparent large increase in the C 
peak. Run No. 13 was made at a reduced pres- 
sure of CO, about one-fourth of the normal 
pressure of about 10-* mm. Run No. 7 is shown 
(See Table II.) 


as Fig. 5. 


AMMONIA 


The ammonia used was produced by heating 
a mixture of ammonium sulphate with a large 
Runs 21, 24 and 26 
were made on a different sample which was 
further dried with KOH pellets, and show that 
the first sample was not satisfactorily dried 
(Fig. 6). It seems likely that the presence of 


excess of calcium oxide. 
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Fic. 6. Run No. 23. Ammonia. 81 volt electrons. Height of 
NH; peak =9.7 cm. 


water vapor contributes to the high intensity 
observed for the Ne» peak, but this point will 
require further observation. It is to be noted 
that there is substantial equality between the 
NH: and NH; peaks, as has been reported 
previously by Bartlett.2 (See Table III.) 


H YDRAZINE 
The sample of hydrazine used was kindly 
supplied by Professor Beckman of the chemistry 
department. It had been distilled three times 
over KOH and then stored for some months in 


TABLE IV. Jon intensities observed for hydrazine bombarded 


NH,* = 100. 
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Fic. 7. Run No. 28. Hydrazine. 113 volt electrons. 
Height of NH, peak = 222 cm. The peak at 27 is not satis- 
factorily accounted for. 


contact with KOH so that its water content 
must have been far below that contributed by 
the walls of the apparatus. 

Because the tube from which the vapor is 
released in the filament assembly is necessarily 
heated by radiation from the filaments, the 
hydrazine seems to be almost completely 
thermally decomposed. Thus the primary ion in 
this case is NH; and not NeH,y Run No. 28 
is shown as Fig. 7. The entry t in the table 
means the peak was definitely present but of 
uncertain height. (See Table IV.) 


with electrons of various energies. Intensities calculated to basis 





NH: NH; NH« CN 27 Ne NoH NoHs NeHs NeHy 

















True 

Run ua ma Volts NH; H: c CH N NH 
28 «21 7.5 113 222 47 a 6S 68 OU 
29 20 7 114 171 79 1.5 —- 17 14 
30 8 4 80 47 49 — — 19 13 
31 30 15 60 eo Bee 6 6 5 
32 15 10 40 40 57 5 2 4 
33 9 ‘ 30 4 25 — 
34 0.1 4.7 20 . _ 


116 100 16 0.67 1 72 t t t 2.7 
123 100 65. 4 1.5 61 t t t 2.9 
105 100 7.4 t t 47 1.5 1 2 3 
102 100 141 — §4 t t t 3 
81 100 20 t ee . -- 
15 100 — — — 








Because of their great stability, nitrogen and 
carbon monoxide yield predominantly molecular 
ions. However, ammonia and hydrazine decom- 
pose readily to yield fractional ions, and can 


2 James H. Bartlett, Jr., Phys. Rev. 33, 169 (1929). 


recombine in many different ways depending on 
collision probability. 

Acknowledgment is made to Professor R. A. 
Millikan, under whose direction this work has 
been done, for his interest and helpful sugges- 
tions. 
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A theoretical investigation is given of the Zeeman effect 
in molecular states intermediate between Hund's case }’ 
and case d’, which are the limiting cases for small and large 
L-uncoupling, respectively. Zeeman patterns for the 3d 
and 4d complexes of the helium molecule are calculated 
and curves obtained for the variation of overall widths with 
K which are compared with experiment. For small K the 
energies are near the case }’ curves, but as K increases a 
breaking away occurs which results in their approaching 


the case d’ curves. When the state in the limiting case d’ 
has the electronic angular momentum quantized anti- 
parallel to the nuclear rotation, a crossing over of the 
magnetic energy curves takes place with the passage 
through a state in which the molecule is relatively insensi- 
tive to the field. This crossing over occurs between K equal 
to 4 and 8 for the A states of the helium molecule consid- 
ered (magnetic field 30,000 gauss). A crossing over which 
occurs in one of the II states is discussed. 





the 
has 


HE experimental work of Millis' on 
Zeeman effect in ortho-helium bands 
shown that a wide variety of patterns is ob- 
tainable of a shape intermediate between the 
shapes of those patterns predicted by Hund’s 
limiting cases for a molecule with electrons 
having orbital angular momentum. In Hund’s 
case b’ the angular momentum vector L is 
tightly coupled to the axis of figure connecting 
the nuclei; its component along that axis is 
quantized and can be designated by the quantum 
number A. On the other hand in case d’ the 
coupling forces between L and the axis of figure 
are negligible; A no longer has meaning as a 
quantum number. It is now ‘the mechanical 
angular momentum of the revolving nuclei which 
is quantized and to which can be assigned the 
quantum number R. In both cases, spin being 
absent, K represents the quantum number for 
the total angular momentum of the molecule. 
The phenomenon of the passage from case b’ to 
case d’ as the rotation of the molecule increases 
is closely related to A-doubling and is called L- 
uncoupling. It has been treated by Hill and Van 
Vleck, Kronig, Van Vleck, Dieke, Kronig and 
Fujioka and others.? The effect of an external 
magnetic field has been discussed for the two 
limiting cases by many authors.’ In the inter- 
mediate case the effective magnetic moment of 
1 Millis, Phys. Rev. 37, 1005 (1931); 38, 1146 (1931). 
See also Crawford, Rev. Mod. Phys. 6, 112 (1934). 

* Hill and Van Vieck, Phys. Rev. 32, 267 (1928); Kronig, 
Zeits. f. Physik 50, 347 (1928); Van Vleck, Phys. Rev. 33, 
467 (1929); Dieke, Zeits. f. Physik 57, 71 (1929); Kronig 
and Fujioka, Zeits. f. Physik 63, 168 (1930); Fujioka, 


Zeits. f. Physik 63, 175 (1930). 
*For a general discussion see the article by Crawford, 


reference 1. 


the molecule varies with K, because of the 
uncoupling, in a way that causes the patterns to 
depart from both of the limiting shapes. It is 
the purpose of this paper to present a theoretical 
investigation of this subject. 


THE PERTURBATION FUNCTION 


We shall start with case b’ as the unperturbed 
system and neglect all the effects of the spin. In 
the bands examined by Millis the spin is unre- 
solved and it is quite satisfactory to regard the 
triplet systems as singlets. This merely means 
that at quite low fields the Paschen-Back effect 
becomes total in all the triplet states and there- 
fore the selection rule AM,=0 is rigorously 
obeyed. The patterns fall on top of each other 
since the term positions are built up around 0 
and +2Arp, for all such triplets. 

Our perturbation function is the sum of two 
parts.‘ The first part, Hj, is that due to the 
rotational distortion of the molecule and serves 
to remove the degeneracy which comes from the 
fact that A can be parallel or antiparallel to the 
axis of figure. The perturbing Hamiltonian 
consists of certain terms in the Schrédinger 
equation for the rotating molecule which are 
usually neglected in separating the electronic and 
rotational parts of the wave function. We take 
the value of the matrix elements of this Hamil- 
tonian from Van Vleck’ or Kronig.® It is 


‘ For a similar computation dealing with spin-uncoupling 
in molecules in a magnetic field, see Hill, Phys. Rev. 34, 
1507 (1929). 

*Kronig, Band Spectra and Molecular Structure, Cam- 
bridge, University Press, 1930. 
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H(A, K;A+1, K)=2(BL,)(A, A+1)[(K¥#A)(K4A+1)}, (1) 
all other elements being zero. In this expression 
2 L(A, A+1)h? 
(BL,)(A, A+1)= | Parte [Pasir (2) 
6 8x°cMr* 


where P(r) is the vibrational wave function of the molecule, M is the reduced mass of the nuclei, 
L,(A, A+1) is a matrix element of one of the components of the orbital electronic angular momentum 
operator L perpendicular to the axis of figure. We can take L,(A, A+1) as equal approximately to 
4(L(L+1)—A(A+1) }', which multiplied by a kind of a mean value of the rotational constant B 
gives us the approximate value of the above integral. 

For the part due to the external magnetic field, we take the Hamiltonian function 


IT,=Av,L., L,=L, sin 6+L; cos @, Av, = He/4amc*. (3) 
Here z is the direction of the field J/, (énf) is a set of Cartesian axes fixed in the molecule with the 
origin at the center of gravity of the nuclei. The ¢ axis coincides with the internuclear line, the & 
axis lies in the xy plane, and @ is the polar angle between the internuclear line and the direction of 
the field. L is the electronic orbital angular momentum operator of quantum mechanics mentioned 


above. To obtain the matrix elements of this Hamiltonian, we utilize the integrals given by Reiche 
and Rademacher* for the symmetrical top. We have for the non-zero elements 


H(A, K ; A, K) = 282, 
H(A, K;A+1, K) = —B8{L(L4+1) —A(A+1) }*((K —A)(K+A+1) }}, 
Av, Ap (K —M+1)(K+M+1)(K—A)(K+A+1)}! 
Pre Okenares —] 
H(A, K;A+1, K+1) (4) 


Av, <r | 


=[L(L+1) —A(A+1) }———— tes scan 
2(K +1) (2K+1)(2K+3) 


H(A, K;A, K+1)= 


where B= Av,M/2K(K +1). 
The unperturbed energies for case b’ are given by the expression 
W,° =f(A)+ BLK(K +1) —A*+L(L+1) —A’], B=h/8x*cMr,?, (5) 


where f(A) is the coupling energy of an L-complex. If one makes the hypothesis of pure precession, 
f(A) becomes A A*, where A is the coupling constant. 


CALCULATION OF ENERGIES 


We shall neglect the terms which are not diagonal in K. If the rotational states involved are not 
close together, in doing this we are neglecting a second order effect in the field which is small. In 
computing overall widths of Zeeman patterns this second order effect would cancel out because it 
only gives a shift of the pattern and a distortion of its shape. Classically we are neglecting the effect 
of that component of the vector L which is perpendicular to the vector K. K precesses around the 
direction of the field and while L precesses rapidly around K, still the effect of its perpendicular com- 
ponent does not average out completely for large fields. Actually we cannot expect the rotational 
energy differences to be always large compared to the Zeeman displacements. However, for extreme 
magnetic levels (K parallel and antiparallel to the field), the effect of the perpendicular component 


® Reiche and Rademacher, Zeits. f. Physik 39, 444 (1926); 41, 453 (1927). 
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is a minimum since this component is then approximately perpendicular to the field. The truth of 
this fact is borne out on examining those elements of our perturbation function which are not 
diagonal in K. They are seen to be smallest for M=-+K. Thus even when the effect of the non- 
diagonal terms is not small, we see that it can be expected that, excepting under unusual conditions 
such as two levels very near each other, their effect is not to change greatly the energy differences of 
levels corresponding to extreme values of M, that is the overall widths, but to distort the shape of 


the energy patterns. 
For L=1 the secular equation is a cubic which can be solved exactly. The three roots are 


Ws =A/2+B(K*+K+1)+8+B[a*/4—ay+y72(2K+1)*]}, 
W=A+B(K*+K) +28, 


where a=A/B, y=1-—8/B, A=f(1). If a>0, we have a normal p-complex and W_ corresponds to 
the = state and W and W, to the two II states. For zero fields these expressions reduce to those 
given by Hill and Van Vleck for A-doubling. The correlation with case d’ is easily found by expanding 
in inverse powers of 2K +1 and is given in the papers by Hill and Van Vleck and by Dieke. 

In the case of L=2, which is of interest in the helium bands, the secular equation is a quintic 
which can be factored into the following quadratic and cubic determinants 


(6) 


(W.°+88 — W) IT,(21) 0 
\(W.°+88 — W) IT,(21) 
0, IT,(21) (W,°+28-— W) 271,(10) |=0, (7) 
IT,(21) (W,°+28— W)| 
0 IT,(10) (W,°— W) 


where /7,= /1,+H1,. The quadratic we can solve at once. For the cubic approximate solutions are 
obtainable which are good for extreme values of the parameters and which reduce to those given 
by Dieke’ for zero field if one makes the hypothesis of pure precession. Since these formulas are 
rather clumsy, we shall not give them here. However if in these expressions we retain only second 
order terms in the L-uncoupling and first order terms in the field, we obtain the simplified approximate 
solutions given below in Eqs. (8) and (9). More directly, we can use the simple perturbation theory. 
Thus for @ large, K and 8/B small, that is near case b’, we have (a >0, K=2) 


WA, =4A+B(K?+K —2)+4(K°+ K —2)B*/3A+88[1—4(K?+ K —2)/3a]= Wasp, 
WIl,=A+B(K?+K+4)+8(2K+1)*B?/3A +28[1—8(2K+1)?/3a], 


Wil, =A+B(K?+K+4) —4(K?+K —2)B?/3A +26[1+4(K2+K —2)/3a], ” 
W>, = B(K?+K+6) —12(K?+K)B?/A+248(K?+K)/a. 
For near case d’, K large, a and 8/B small, we likewise have 
Wi:=A+B(K —2)(K —1)+ B(12a— 3a’) /4(2K+1)+48(K +1) — $8a?/(2K+1 
Wb_2=A+B(K+2)(K+3)+B(12a—3a*)/4(2K +1) —48K+ }80°/(2K+1 
Wé,=5A/2+B(K —1)K+B(36a—9a?)/8(2K +1)+26(K +3) —(9/8)Ba?/(2K+1), 9 


Wié_,=5A/2+B(K+1)(K+2) —B(36a—9a*)/8(2K +1) —28(K —2)+(9/8)a?/(2K+1), 
Wi9=3A+BK(K +1) —B(12a—3a?)(2—a)/2(2K+1)?+68+68a7(3—a)/(2K+1). 


Here we are using the subscript a to indicate the higher energy component of a A-doublet in the 
normal case and } for the lower energy component. The cubic gives the higher components as well 
as the = level to which we have also given the subscript a. 6 in Eq. (9) means that it is a 6-complex 


7 See Dieke, reference 2. 
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(L=2, case d’) and the subscript indicates the quantized value of L along the direction of nuclear 
rotation. These expressions reduce to the Zeeman expressions for the two cases for limiting values 
of a. 

It should be pointed out that we can just as well begin with case d’, in which case we assume the 
interaction Hamiltonian to contain the square of the cosine of the angle between the L vector and 
the axis of figure. The matrix elements are readily obtained from the relations given in the articles 
by Van Vleck and Hill. On neglecting the elements not diagonal in K, we obtain the same secular 
equation as by the above method. 

In the case of the cubic we have the following simple way of obtaining values for the magnetic 
energy levels from the zero-field levels. In order to use the data given by Fujioka, we take the 
unperturbed (case 5’) energies as approximately 


W,°=A,+B(K+})*—8'(K+}), A=0, 1, 2. (10) 


Then if Wo, Wi, We are the zero-field roots of the cubic and if we neglect the squared and cubed 
terms in 8, we can write our secular equation in the form 


—(Wy— W)(W,— W)(W2.— W) = 28[ (Wo — W) (We — W) +4(W0 2 — W) (WY — ) 
—48B°K(K + 1l)ay?+12BK(K + 1)ayo( We — W) +4B(K2+K — 1) ail Wo — W) ). (11) 


To obtain the energies corresponding to Wo, we positive terms, we have remaining three states 
substitute W, for W throughout excepting in the for K odd and two for K even. These correspond 
first bracket on the left and then solve the result- 
ing linear equation. The a’s are constants, nearly 
unity, introduced and evaluated by Kronig and 
Fujioka in their zero-field secular determinant 
so as to get the best fit with experiment. They 













































































allow for the fact that (BP,)(A, A+1) in the 
first part of our perturbation function is not SN _|---|---|---]t2-*] Ab 
1 j i riceK nd 
accurately equal to }B[L(L+1)—A(A+1) }}. aa ——, af 
A, oe ees oes oe... “1 
> 
APPLICATION TO THE HELIUM BANDS : T 
Experimentally we find that half of the ” 
rotational states are missing. This is because Be 
He: is a homonuclear molecule the nuclei of 5} 
ae c 
which obey Bose statistics and therefore the wf -J. 

: ’ ° . = — | woe am MIS OK TTy 
allowed states are symmetrical in the nuclei. MK] ong 
a , = n 
Since the states we are interested in are also 7. M=-K 
odd (u), they are all negative under reflection eens Oe AR GR 4. 
in the origin and we must omit the positive ones. PY 2 nated wateneen eateries weet Un 
It is easy to see just which the negative ones are. 
dS Sy, ~~ : io o |4 |e le he |20 Kee 
For =,* states® the terms are alternately positive ono 
and negative beginning with a positive term for al \ b> 
K=0. The two rotational terms of a given K Pe a ee ee 
for II and A states are always one positive and Fic. 1. Magnetic levels of a d-complex for the limiting 


ee ae ; —— : a ole. cases drawn around the zero-field energy (M=0). Only 
asia en _— —_— - K be ing alter- M=0, +K levels are plotted. The solid waived are case v 
nately positive and negative. Discarding the and the broken curves case d’. The correlation between the 
states is that for a>0. The energy scale is in units of Av,. 

*See Weizel, Handbuch der Experimental Physik, Band- K isthe quantum number associated with the total angular 
spektra, p. 136. momentum of the molecule. 
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Fic. 2. Magnetic levels of the 4d complex of the helium 
molecule calculated from the present theory using Kronig 
and Fujioka’s constants. a is approximately 10. The cross- 
ing over of the levels is due to the passage from quantiza- 
tion of the orbital angular momentum of the electron along 
the internuclear axis to quantization along the axis of 
nuclear rotation. During this process the component of 
the magnetic moment along K passes through zero thus 
passing through a state in which the molecule is very 
insensitive magnetically. The circles give the observed 
widths. In this figure and in Fig. 3, Av, is 2.415 cm™ 
which corresponds to a field strength of 30,000 gauss. 


to our cubic and quadratic given above because, 
our Hamiltonian being symmetrical, states of 
unlike symmetry cannot perturb each other.*® 

In Fig. 1 is given for comparison the overall 
widths of the magnetic energy levels for the two 
limiting cases, drawn using the expression 


E=28A? (12) 


for case b’ and 


E=p(K(K+1)+Z(L+1)—R(R+1)] (13) 


for case d’. In Fig. 2 are the calculated curves 
for the 4d complex of molecular helium and in 
Fig. 3 those for the 3d complex. All the curves 
show clearly the rapid change which takes place 
in passing from case 5’ with low K to case d’ 


*In Mulliken’s notation the states A,, A», Ile, Ils, 2. 
considered in this paper are designated by A,, Ay, Hy, I., Z., 
respectively. Rev. Mod. Phys. 3, 94 (1931). 


JULIAN K. 
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Fic. 3. Magnetic levels of the 3d complex of the helium 
molecule calculated using Kronig and Fujioka's constants. 
a being about 20, the binding to the internuclear axis is 
tighter than it is for the states in Fig. 2 and the transition 
to case d’ is therefore slower. 


with high K. The change is more rapid for the 
4d complex since its angular momentum vector 
is less tightly bound to the axis of figure. The 
crossing over of the extreme levels occurs when 
the component of the magnetic moment along 
K changes its direction with respect to K. The 
circles represent experimental points obtained 
from Millis’ data by interpreting them in the 
light of the present theory. The values for the 
overall widths so obtained are not all self-con- 
sistent; we have tried to select those values 
which were the most reliable from the experi- 
mental point of view. The differences which 
exist are perhaps due to intensity distributions 
and failure to obtain well resolved patterns in 
some cases. 

The II, states, for which unfortunately there 
are no experimental data, are of special interest 
since they are the only states for which the cal- 
culated curves do not lie between the limiting 
curves. An understanding of this behavior is to 
be found in examining the solutions for small and 
large K using Eqs. (8) and (9). The equation for 








PERTURBED SERIES 


II, shows that for small K the calculated curve 
for M=K will lie below the limiting case }’ 
curve (a>0) and that the difference between the 
two will increase rapidly with K. The equation 
for 59, which is the limiting case for a II, state 
with uncoupling, shows that for a>3 the cal- 
culated curve lies below the limiting curve but 
that the difference decreases rapidly with in- 
creasing K. For M=K both limiting curves lie 
for all K above the curve for M=0. The cal- 
culated curves start out for small K below the 
first limiting curve and then approach the 
second limiting curve from below. For the 4drIl, 
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state the extreme levels cross over at about 
K=10. For a slightly greater than 20, the 
extreme levels would cross over first for small K 
and a second time for very large K. 

Curves of a similar nature could be drawn for 
the magnetic levels of the three states of a p- 
complex using Eq. (6). The II state correlated to 
the 6_; state would show a crossing over of the 
levels for suitable a and K in a manner very 
similar to that shown in the cases described above. 

The writer wishes to express his sincere thanks 
to Professor Kemble and Professor Van Vleck for 
valuable suggestions and criticisms on this work. 





MAY 1, 1935 PHYSICAL REVIEW VOLUME 47 


On Perturbed Series, Especially in C III, BI and OIV 


N. G. WaHIreLaw AND J. E. Mack, Department of Physics, University of Wisconsin 
(Received March 6, 1935) 


Consideration of possible extra-configurational perturbations in atomic spectral series 
shows that often some series are preferable to others for the determination of series limits. 
In particular: In C III, several 2snd *D terms are recalculated, but Edlén’s series limit value is 
verified in spite of perturbations; in B I, the term values are increased by 90 cm™ or 0.011 volt, 
so that B I 2s*2p*P,—B II 2s? 'So = 66,930 cm =8.257 volts; and in O IV, the old 5s *S, term 
is replaced by a new one at 85,440 cm™', ?4p*P and ?6f*F are repudiated, and ?7f*F and ?8f *F 








are verified. 


LTHOUGH the term values in most of the 

spectra of the lightest elements have re- 
cently been determined with a high degree of 
accuracy and completeness by Edlén,' there are a 
few instances in which series limits may be 
slightly altered and term assignments recon- 
siderec. after a study of the influence of extra- 
configurational perturbing terms.?:* The pur- 
pose of this note is to point out the general 
preferability of some series over others in the 
determination of series limits, and to make minor 
alterations in three of Edlén’s spectra. 

An irregular relationship between the term 
values and the corresponding quantum defects 
for a spectral series generally is an indication of 
the presence of a perturbing series, of which 


1B. Edlén, Nova Acta Reg. Soc. Scient. Upsaliensis, 
Ser. IV IX, No. 6 (1933). 

? A. G. Shenstone and H. N. Russell, Phys. Rev. 39, 415 
(1932). 

*N. G. Whitelaw, Phys. Rev. 44, 544 (1933), 


usually only one or two terms appreciably affect 
the series in question. The extra-configurational 
series members have the same Laporte parity 
and (in the LS coupling case, which is all that 
need be considered in Edlén’s spectra) the same 
multiplicity and L value as the series under 
consideration. Irregularity arises when the two 
series overlap or nearly overlap, in which case the 
members of the one series share their properties 
with neighboring members of the other with a 
consequent ambiguity as to configuration. Thus 
the members of a series may show an energy 
trend leading toward a false series limit, even 
though the term or terms mainly responsible for 
the complication may remain undiscovered. Al- 
though perturbing terms can be pointed out for 
almost all of Edlén’s spectral series, we shall list 
here only those cases in which an analysis ap- 
pears to lead to definite improvements in term 
values. 
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C III: The 2snd *D series has been used to 
determine the series limit for C III. In the case of 
series members for n=7, it seems impossible to 
assign J values and consequently the energies 
should be calculated with respect to the center 
of gravity of the 2s2p *P term. If this is done one 
obtains the following recalculated term values: 


> 


¥ n 
2s7d*D 20,575 cm™ 6.9283 
8d 15,722 7.9258 
9d 12,412 8.9202 


The resulting »—n* versus v plot exhibits quite 
a curvature near v= 0 as is to be expected because 
of the presence of the 2p4p *D term near the 
series limit. The reason for the smallness of the 
shift of the 2s3d *D term due to perturbation by 
2p3p *D is not apparent. In spite of the perturba- 
tions present, the *D series is probably the most 
reliable from which to determine the series limit 
because of its length. The 'F and *F series are 
very badly perturbed as is easily seen from the 
variation of the quantum defect with v. Also the 
separation of the components of 2s4f *F would 
be nearly zero except that because of the per- 
turbation, the term acquires part of the multiplet 
width of 2p3d *F. 

B I: The series limit of B I was determined by 
Selwyn.‘ He fitted both the *S and *D series to a 
Ritz formula and took an average value for the 
series limit. However, upon examination of the 
series, one finds that a term, 2s2p* *S lies near to 
the 5s *S series member and undoubtedly per- 
turbs that term. The 2s2p* °S term has not been 
found experimentally but by examining the 
2s—Zp energy differences in members of the 
isoelectronic sequence, one can place it within the 
range 5000 to 9000 cm™', whereas 5s *S lies in 
the neighborhood of 6694 cm~'. The nd 2D series 


*E. H. W. Selwyn, Proc. Phys. Soc. 41, 392 (1929). 
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Anew 3.°8. MACK 

is perturbed by a 2s2p* *D term which lies much 
deeper than any of the *D regular series members 
and consequently the *D series is not seriously 
perturbed. It seems more reasonable therefore to 
determine a series limit from the *D series alone, 
and if this is done, we find a value about 90 cm™ 
greater than the Selwyn value. 

O IV: The series limit of O IV is quite accu- 
rately determined by the (experimentally) un- 
perturbed *D series. However, several other series 
are badly perturbed, and we think involve several 
wrong classifications. (a) The quantum defect of 
the 5s *S term as classified is impossible. The 
perturbing term, 2s2p3p *S lies below 4s *S and 
consequently should not affect 5s 2S as much as 
4s *S. We have identified a line® of intensity 
= 1, \= 185.544A, v= 538,956 cm“, attributed by 
Edlén to O IV on experimental grounds, which 
we think gives a more probable quantum defect 
for the 5s 2S term. The new term value for 5s 2S 
becomes 85,440 cm™. (b) Edlén questions his 
classification of 4p *P and one can reason that it 
must be wrong since the relative quantum defect 
for 3p *P and 4p *P is opposite to what it should 
be, because of the presence of the perturbing 
2s2p3s *P term between them. (c) Edlén also 
questions his nf *F terms for »=6. The terms for 
n=7 and n=8 seem all right but 6f *F must be 
wrong, The nearest perturbing term, 2s2p4d ?F, 
lies just below 8f ?F, and 2s2p3d *F lies below 4/ 
?F so the abnormal quantum defect of the present 
6f *F term would seem to disqualify it. 

It is easy to determine what series are likely to 
be free from disturbance by extra members in 
any particular case, although there seems to be 
no rule for them except to study the list of terms 
yielded by each configuration. 


* B. Edlén, Zeits. f. Physik 85, 85 (1933). 
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The Normal Helium Atom 


J. H. Bartiert, Jr., J. J. Grppons, Jr. anp C. G. Dunn, Department of Physics, University of Illinois 
(Received March 6, 1935) 


A least-squares criterion for the goodness of approximate wave functions is proposed. The 
root-mean-square energy deviation for the six-term Hylleraas function is calculated to be 
3.5 volts. It is shown that no ascending power series in the variables r;, rz and ry». can be a 


forma] solution of the Schrédinger equation. 





HE problem of determining the Schrédinger 

eigenfunction and eigenvalue for the lowest 
state of helium is one of fundamental importance 
to the theory of atomic structure; if it were 
solved, we should know whether or not the pres- 
ent formulation of the nonrelativistic many-body 
problem is correct. The subject has received the 
attention of several writers, notably Slater,' 
Kellner® and Hylleraas.’: *: ° Slater, by numerical 
methods, arrived at an eigenvalue for the atom 
and a charge density for an s-electron, but did 
not represent the wave function analytically; 
Kellner and Hylleraas used the Ritz variational 
procedure, in which one guesses at the form of the 
eigenfunction y and allows certain parameters to 
vary so as to minimize the integral {yHydr. 
The latest work on this by Hylleraas® resulted in 
an approximate eigenvalue lower than the ex- 
perimental term value by an amount which could 
be attributed to relativistic effects. 

The Ritz method results in an upper bound to 
the true eigenvalue; a lower bound has been 
sought by Weinstein®: 7 and MacDonald.* Since 
their methods seem somewhat difficult to justify 
rigorously, we shall in the present paper try to 
show how a least-squares method can aid one in 
judging the relative merits of different approxi- 
mate wave functions. In addition, we propose to 
prove that no eigenfunction of the Hylleraas 
type can be even a formal solution of the wave 
equation. 


1J. C. Slater, Phys. Rev. 32, 349 (1928). 

2G. W. Kellner, Zeits. f. Physik 44, 91 (1927 

+E. A. Hylleraas, Zeits. f. Physik 48, 469 (1928). 

*E. A. Hylleraas, Zeits. f. Physik 54, 347 (1929). 

*E. A. Hvileraas, Zeits. f. Physik 65, 209 (1930). 

*D. H. Weinstein, Phys. Rev. 40, 737 (1932); Phys. 
Rev. 41, 839 (1932 

7D. H. Weinstein, Proc. Nat. Acad. Sci. 20, 529 (1934). 

*J. K. L. MacDonald, Phys. Rev. 43, 830 (1933); 
Phys. Rev. 46, 828 (1934). 





The method of least squares suggested by 
Boussinesq® consists in adjusting the constants in 
an approximating function ¥ so as to make the 
integral J= {[ /7—E)y Pdr as small as possible." 
Let us set (/7—E)y= ey, where « will in general 
be a function of x, y and s. For any approximate 
eigenvalue E, ¢(x, y, z) is an indication of how well 
the wave equation is satisfied at any point. The 
integral ]= [Yedr may then be looked upon as 
a mean square energy deviation." We shall calcu- 
late its value, using an eigenfunction due to 
Hylleraas. 

For convenience in computation, the final 
eigenfunction given by Hylleraas in reference 4 
will be used. A further refinement was made in 
reference 5, but the value of E was lowered by 
only 0.013 volt, so presumably the value of J 
would not be influenced much by such a change. 

The wave function is, then, ~= o(ks, ki, ku), 
where 


o(s, t, u)=e-**[egteutcnP testes +o }. 


The variables are s=rj;+f%, t= —ri+% and 
u=f\2; the constants are co=1, ¢c,=0.0972, 
C2= 0.0097, C3= — 0.0277, c,= 0.0025, Cs 
= —0.0024, and k= (L/2M)=0.9089645. 

Since it is purely a routine matter to apply 1/ 
to this wave function and to calculate f(//y)*dr, 
we shall give only the result. The value of this 
integral is 2.1115. If J be minimized” with re- 
spect to E, then E= fy~Hydr, and J= f(Hy)*dr 
— E*. The constants in ¥ may now be varied to 


* J. Boussinesq, Theorie de la chaleur (1), p. 316. 

1° We should like to express our appreciation to Mr. 
R. K. Cook for calling our attention to this method 

One could, alternatively, weight the square of the 
deviation ef to emphasize regions where the eigenfunction 
is large. 

2 This fixes the relation between E and the constants 
Co + -¢. Up to this point EZ can be assumed as independent 
of the constants. 
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make either E or J a minimum. (One would in 
general expect that E and J do not take on 
minimum values simultaneously, since the wave 
function used is only an approximate one.) We 
shall suppose it is £ that is made a minimum. The 
value given by Hylleraas is E= — 1.45162 Rh, so 
that J=0.0043. It follows that /'=0.065 Ra, 
corresponding to about 3.5 electron volts. This 
root-mean-square deviation would vanish if the 
wave function were actually a solution of the 
differential equation. 

One cannot, however, necessarily assert that 
the true eigenvalue lies about 3.5 electron volts 
below that calculated by Hylleraas. The r.m.s. 
deviation may be viewed simply as a measure of 
the goodness of the approximate wave function 
used. That is, if two different approximate func- 
tions should give the same value of E, then one 
would suppose that function with the smallest 
r.m.s. deviation to be the better of the two. 

Since the actual deviation for some regions of 
space will in general be larger than the r.m.s. 
deviation, it seems to us that the above value of 
3.5 electron volts is not small. It can be easily 
verified, in fact, that the deviation becomes 
infinite at any of the singularities, r,=0, r= 0, 
or r2=0. The eigenfunction is thus probably 
very badly in error in the general neighborhood 
of any of these points. This prompts an investi- 
gation of the relationship of the function to the 
differential equation. It might seem that matters 
would be improved by using a power series in- 
stead of a polynomial (times an exponential fac- 
tor). We shall now prove that no series in ascend- 
ing powers of r;, re and 7,2. can be a formal solution 
of the wave equation. 

For convenience, put x=7r);, y=fe and 2= frye. 
the wave equation is then: 


Wert (2/x)Wet yy t+ (2 y)vyt+2¥.2+(4 s)v; 
+[ (x? — y?+3*)/xs j¥..+([(y? —x°+3°)/ys Wy; 
+[(A/4)+(1/x)+(1/y) — (1/28) Wy =0. 


GIBBONS AND DUNN 


os 
If we expand, ¥= }-cm,x'y"z", substitute, and 
0 

equate coefficients of powers, we find the recur- 
sion formula, 
(+2) (14+34+m)ci.2, m,n 

+(m+2)(m4+-3+N)o1, mid, n 

+ (n+2)(2n+64+]+m)cr. m. nad 

—(+2)(m+2)cir2, m—2, na2 

— (m+2)(m+2)ci~2, m2, naa (A/4)Cimn 

+Ci41, m, aw HCr, m41, 2 — $C, m, 041 =0. 

Substitution of various sets of values of 1, m, and 


n will result in specific relationships enabling us 
to determine the c;,,,'s. We list a few as follows: 


l= —1, m=2, n= —1; cy, =0, 

[= — . m =(), n =(); 2100 +Coo0 = Q, 

j= z.. m =(), n=—1 ° 5« 101 ~~ x 109 = O. 
Assuming that Cooo= 1, we thus obtain on the one 
hand ¢,9,:= 0 and on the other ¢,9;= — (1/20). The 


fact that at least two inconsistent values of the 
same coefficient may be obtained means that the 
power series in question will not satisfy the differ- 
ential equation, even formally. 

In our opinion, the fact that a minimum value 
of E exists for functions of the Hylleraas type 
does not shut out the possibility that appreciably 
lower values of E may result from a use of func- 
tions which come closer to satisfying the differ- 
ential equation. To improve the existing situation 
one may either construct functions with lower 
values of E and J, or else look for exact solutions. 
Attempts along the latter line are now being 
made. It turns out to be rather easy to find formal 
solutions, and whatever difficulties there are 
come in with the imposition of the boundary 


conditions. 
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Atomic Wave Functions for Some Excited States of Helium 


W. S. Wicson anv R. B. Linpsay, Brown University 
(Received January 16, 1935) 


Normalized one-electron wave functions and the corresponding fields and energy parameters 
have been computed by the self-consistent field method for the following excited states of 
helium: (1s)*, (1s)(2s), (1s)(2p), (2s)*, (2p)? and (2s)(2p). The results have been compared 
with those obtained by other methods and with experimental values. 





INTRODUCTION 


HE “‘self-consistent field’’ method of calcu- 
lating atomic wave functions has been ap- 
plied successfully to a number of atoms. So far 
as the authors are aware the following constitutes 
a reasonably complete list. O** to O, Ne, Nat, 
Sit‘, Cl-, K*, Cut, Cs* have all been computed 
by D. R. Hartree and collaborators with high 
accuracy. He, Li and Li*, Be*+, Be** and Be, B, B** 
and B**, F~ and F, Ca**, C, N, Na and Hg have 
been calculated with less emphasis on extreme 
accuracy.? Rather recently Sit? and Si** have 
been worked out approximately by H. L. Donley 
at Brown University (at present unpublished). 
The present paper contains the results of 
similar calculations for certain excited states in 
helium. Wave functions for some excited states 
in helium-like atoms have already been obtained® 
by the use of variational methods. Opportunity 
is taken here to compare the results of the “‘self- 
consistent field’’ calculations both as to wave 
functions and energy parameters with those of 
other methods. 


METHOD OF THE CALCULATIONS 


The method of the “self-consistent field” has 
been described so often‘ and is now so familiar, 


1 See, for example, Hartree, Proc. Roy. Soc. Al41, 282 
(1933) for a good survey of the method and some results. 
Hartree and Black, Proc. Roy. Soc. A139, 311 (1933); 
Hartree, Proc. Roy. Soc. A143, 506 (1934); J. McDougall, 
Proc. Roy. Soc. A136, 549 (1932); A138, 550 (1932). 

? Hartree, Proc. Camb. Phil. Soc. 24, 89 (1928); J. Har- 

reaves, Proc. Camb. Phil. Soc. 25, 75 (1929); F. W. Brown, 
| H. Bartlett and C. G. Dunn, Phys. Rev. 44, 296 (1933); 

artree, Phys. Rev. 46, 738 (1934); C. C. Torrance, Phys. 
Rev. 46, 388 (1934); E. H. Kennard and E. Ramberg, 
Phys. Rev. 46, 1034 (1934). 

? See, for instance, C. Eckart, Phys. Rev. 36, 878 (1930); 
J. P. Vinti, Phys. Rev. 37, 448 (1931); F. G. Fender and 
J. P. Vinti, Phys. Rev. 46, 77 (1934). Also E. A. Hylleraas, 
see reference 12. 

*e.g., Hartree, Proc. Camb. Phil. Soc. 24, 89 (1928). 


that it is unnecessary to redescribe it in detail 
here. The limitations of the method have also 
been amply dealt with elsewhere.* We shall con- 
tent ourselves with discussing briefly the tech- 
nique used in the application of the method in the 
present calculations. 

The states of helium investigated are: (1s)* 
(already worked out by Hartree, but repeated 
for the sake of completeness and to tabulate the 
wave function using smaller intervals of (r), 
(1s)(2s), (1s)(2p), (2s)*, (2p)* and (2s)(2p). 

The only important point in which the method 
used in this paper differs from that recommended 
by Hartree is that direct integration was used 
throughout in preference to the so-called varia- 
tion method* in obtaining the join of the inward 
and outward integrations of the appropriate 
Schrédinger equation. The extra work involved 
in this procedure was probably compensated for 
by the greater assurance of the elimination of 
errors. In any case the labor involved is by no 
means so great as might at first be thought. The 
successive approximations to the self-consistent 
field were in general made in each case by using 
the final field of the previous approximation. An 
important consideration is, of course, the assump- 
tion of the initial field, i.e., the initial Z, (the 
effective nuclear charge for potential.t Several 
schemes have been investigated. The one gen- 
erally followed may be briefly summarized as 
follows. 

The wave function of He* in the (1s) state was 
first tabulated and the corresponding field calcu- 
lated. This at once provided an approximate 
initial Z, for the (2s) electron in the singly ex- 
cited (1s)(2s) state. From the distribution of 


* J. C. Slater, Phys. Rev. 35, 210 (1929); V. Fock, Zeits. 
f. Physik 61, 126 (1930). 

* Cf. Hartree, reference 4, p. 95. 

+ Cf. Hartree, reference 4, p. 115. 
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Fic. 1. Plot of the normalized wave function P. 


charge of the (2s) electron the initial Z, for the 
2s electron in the state (2s)* was computed. 
Similarly the distribution of charge for the (1s) 
electron in the (1s)(2s) state was used to com- 
pute an initial Z, for the (2p) electron in (1s)(2)), 
and the (2s)(2p) case was handled by starting 
from (2s)? and the (2p)? case from (2s)(2p). The 
number of successive approximations required to 
attain self-consistency in these various cases 
varied from two to nine. 

An alternative method of estimating the initial 
Z, holding out some promise of value consisted 
in extrapolating the wave functions for the atoms 
in the second period of the natural system as 
calculated by Brown, Bartlett and Dunn.? If the 
reciprocals of the positions of the maxima for the 
(2s) wave functions from boron to neon as given 


Tape |. Normalized wave 
1a}? le)( 2s 2s 28) (2p 2 ls)(2p 
r ls Is) 2s 2s le 2p 2p le 2p 
0.00 0.0000 0.0000 0.0000 0.0000) 606.0000 0.0000 0.0000 0.0000 0.0000 
065 2152 240 O49 076 O701 0023 0021 2556 0007 
10 3900 4597 093 1382 1265 .0O80 .OO79 4626 .0025 
15 5308) =— 8240 1349 INO 1705 0190 0169 6279 0055 
2 28 7530 1628 2228 20388 .0322 0286 7576 0003 
3 TOM 92s lye 2654 429 0656 0582 9304 0193 
4 8853 1.0110 2178 2759 2526 1057 0938 1.0157 0317 
J 9241 1.0858 2213 2620 2403 1497 1329 1.0306 0459 
iy 9206 1.0189 2141 2303 2118 .1054 1736 1.0215 .0617 
7 9123 9746 lgsy 1S59 1719 «6.24138 2145 9759 0785 
8 S708) =O 134 1778 1330 1242 2859 23845 0134 0062 
4 S374 8427 1522 0750 0720 3284 «2027 )«=—M416 SO.1145 
1.0 7TR03 7680 1235 O44 0174 3681 3287 .7659 1332 
1.2 6859 6202 0601 1065 o919 4371 394 =—.6166 1709 
4 5836 0=— 4880 0069 1184 140 «64913 «4440 48260 «6.2081 
1.6 4891 3746 O737 3155 2338 .5306 4836 6.3700 2440 
18 4054 2837 1377 3953 3592 5561 5117 2704 2776 
2.0 3331 2122 —.1972 4576 4198 5606)«=— 5208) =—.2084 3085 
2.2 2718 1572 2514 5034 4664 5728 5390 «6.1539 «=. 
24 2205 1155 2005 5345 5002 «4.5676 «66408 «= .1128 = «3809 
2.6 1780 OmM43 M13 54528 $227 5557 5365 0822 3821 
2 1431 0612 3768 ydOS 53455 S388 5272 0595 4000 
3 147) = 442 4061 5590 5400 «6CBI80) «0140 S420) S445 
3.2 0916 Osis 4206 5507 5378 4046 4978 0308 4258 
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by them are plotted and extrapolated to beryl- 
lium, a ‘‘stretch’’ factor may be obtained allow- 
ing one to compute very simply the normalized 
(2s) wave function for beryllium, which has in 
the normal state 2 (1s) electrons and 2 (2s) 
electrons. The latter move in the field of a nucleus 
of 4 positive charges and the 2 (1s) electrons. 
Hence to a first approximation the field for large 
r at least may be considered to be like that for 
the (2s) electrons in helium. This choice proves 
to be a rather satisfactory one as Fig. 1 shows. 
This presents the plot of the normalized wave 
function P corresponding to the assumption just 
indicated (curve 7) compared with that obtained 
by using the initial Z, from the method described 
in the previous paragraph (curve //) and with 
the final P function (curve J//) for the (2s 
state. In this case the alternative method yields 
a rather better first approximation. *® 


RESULTS AND DISCUSSION 


The principal results of the calculations are 
exhibited in Tables I, II, II] and IV, presenting, 


*In this connection it may be remarked that a study 
was made of the use of the Fermi statistics as applied by 
Baker (cf. E. Fermi, Zeits. f. Physik 48, 73 (1928); E. B 
Baker, Phys. Rev. 36, 630 (1930 order 
This showed that for the heavier atoms the values ob- 
tained from Baker's table provide an excellent initial field 
for Hartree calculations. In particular in the case of Si** 
the former values of Z, agree with those for the self 
consistent field to within a maximum deviation of 2 percent 


to tons of any 
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TABLE Ia. Values of P/r'** for small r. 


(28) (2p) (2p)* (la)(2p) 


r 2p) (2p) (2p) 
0.00 0.00000 0.00000 0.00000 
05 O4583 04143 01336 


10 08721 734 02547 


respectively, the normalized wave functions 
(P values) for the various excited states, the 
contributions to Z (the effective nuclear charge 
for field intensity) of the various electrons, their 
contributions to Z,,’ and finally the values of the 
one-electron parameter ¢. The figures in the first 
three tables are given to four decimal places, the 
accuracy of the computations being such that 
any error in P may be expected not to exceed 2 in 
the fourth place, and any error in Z or Z,, 1 in 
the fourth place.* All values in the tables have 
actually been found by straightforward calcula- 
tion: none has been interpolated. In fact the 
intervals used in the computation were in general 
half the size of those given in the tables. Table 
I(a) contains the values of P/r'*' for small r 
in the case of the 2p electrons: these are of value 
in starting the outward integration. Correspond- 
ing values for the (1s) and (2s) states need not be 
tabulated since they can be computed with suff- 
cient accuracy directly from P. 

It may be noted in passing that the values of 
Z and Z, for (1s)* do not agree precisely with 
those of Hartree (first paper in reference 2; 
note that the values there quoted refer to the 
two electrons and must be divided by 2 for com- 
parison with those of the present paper). It 
should be remarked that Hartree in the case of 
helium limited the consistency between initial 
and final fields to 2 in the third decimal place 
which falls considerably short of the standard of 
accuracy set in the present paper. Incidentally as 
a check on the accuracy throughout the work, 
Sodrf-°P?/rdr has been calculated for each 
electron. This gives Z, for r=0 which should be 
unity for each electron. In the present work, it 
does so in every case to within 5 in the fifth 
decimal place. 

7 Strictly speaking these are unnecessary. Their utility, 
however, dictates their inclusion here. Cf. Hartree, refer- 
ence 1. 

* This means of course that the maximum error to be 
expected in the total Z and Z, (not tabulated since they 


are readily obtainable by summing the contributions of 
the individual electrons) is 2 in the fourth decimal place. 


The values of ¢ in Table IV are given in general 
to three decimals and may be assumed to be cor- 
rect to that place: the extra place used in obtain- 
ing them has been discarded. Those values listed 
to four decimal places are expected to be in error 
by no more than 2 in the last place and the error 
for the (1s)(2s) and (1s)(2p) is probably con- 
siderably less. Table IV also includes —2* and 
values of « computed by other methods as men- 
tioned in the Introduction. Comparison with 
experimental results is rather difficult since the 
latter are usually given in terms of the ionization 
energy of the atom in a given state** while the 
present values are one-electron energy param- 
eters. These can of course be used in connection 
with the wave functions to compute the total 
atomic energy values (cf. Hartree and Black’). 
For (1s)*, (2s)? and (2p)* this calculation can be 
carried out approximately by assuming that « 
itself is a good approximation to the ionization 
energy. The computed total energy for these 
three states then appears as 5.836, 1.4606 and 
1.3861, respectively. The experimental value is 
known definitely only for the first and is 5.808. 
Some evidence with respect to the second is at 
hand in the recent experiments of Whiddington 
and Priestley® on energy losses in inelastic colli- 
sions in helium. Their work indicates a loss of 
59.25+0.12 volts in a single collision of an elec- 
tron with a helium atom. Now the difference 
between the total energy in the (1s)* state (5.808) 
and that in the (2s)* state (1.4606) is 4.347, cor- 
responding to 58.81 volts, which is very close to 
Whiddington and Priestley’s value, making it 
seem likely that the transition corresponding to 
their indicated energy loss is really (1s)? — (2s)*. If 
the present calculations are at all correct it is 
much more likely this than (1s)*—(2p)*, al- 
though the possibility is not excluded that it 
may be (1s)*— (2s)(2p). Not until the total energy 
for the (2s)(2p) is calculated can this point be 
decided. 

In the (1s)(2s) and (1s)(2p) states since the « 
value for the (1s) in each case is much greater 
numerically than that for the 2-quantum elec- 
tron, we may expect that the removal of the 


* Cf. Hartree, reference 4, p. 95. 

* e.g. Bacher and Goudsmit, Alomic Energy States. 

*R. Whiddington and H. Priestley, Proc. Roy. Soc 
A145. 462 (1934). 
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TABLE II. Contributions to Z of the various electrons in helium. 
(1s)? (1e)(2a) (2a)? (2a)(2p) (2p)? (1e)(2p) (1s)? (1a)(2a) (2a)? (28)(2p) (2p)? (1a)(2p) 
r (le) (la) (28) (28) (28) (2p) (2p) (1s) (2p) r (le) (le) (2a) (2a) (28) (2p) (2p) (la) (2p 
0.00 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000] 3.0 0058 0006 .8548 5023 6457 4011 4756 0000 .8059 
05 9902 9989 9999 9999 9999 1.0000 1.0000 9989 1.0000) 32 0036 0003 (8198 5309 5875 3496 4243 0003 _7705 
‘10 9044 9022 9906 9993 9004 1.0000 1.0000 9921 1.0000| 3.4 0023 .0002- (7812 4717 (5304 (3030 3785 0001 .7335 
‘15 ©9836 «9772 9989 9980 9983 1.0000 1.0000 9770 1.0000| 3.6 0014 0001 7399 4150 4754 2611 3325 0001 6053 
20 ©9662 9533 9978 9958 9965 1.0000 1.0000 9527 1.0000| 3.8 0009 0.0000 6966 3642 4234 2240 2022 “6564 
4.0 0006 6522 3170 3747 11913-2558 6163 
3 9181 8812 9945 9807 9914 9997 0.9008 8798 1.0000 
4 S414 786299009823 9851 9900 9002 7838 9000) 44 0002 5625 2361 2888 «1379 1987 5398 
5 .7500 6806 9853 9749 9790 9073 9079 6774 .9008| 48 0001 A758 «£1725 2184 0981-1447 4647 
6 6726 S744 9804 9688 9738 9043 9955 5706 9905| 5.2 3955 1240 1624 0690 1068 3951 
‘7 S875 ATAT «9761 «9644 «9700 9806 9918 4704 (9900/ 5.6 3236 «0878 «1191 10480 ©0780 3320 
‘8 BO71 3854 9725 9618 9679 9826 9863 3810 .9982| 6.0 2611 0615 0862 0331 0564 2758 
‘9 4832 3081 9608 9607 9669 9731 9787 3038 9971 
10 «= 3670-2432 (9679 ©9805 «9667 9610 9601 12392 9956| 6.8 1638 0202 0437 0154 ©0288 1848 
7.6 .0985 0134 0214 .0069 0143 -1193 
12 2579 1467 9661 9508 9662 0311 9429 1434 9000) 84 0573 0080 0102 0031 0069 ‘0752 
14 1773 0854 9650 9542 9619 (8876 9077 0831 9837 9.2 0823 0026 0047 0013 0083 0460 
16 1197 0484 9655 9397 9502 (8349 8644 0468 9735 | 10.0 0178 0011 0021 0006 ©0015 0276 
18 0797 0200 9632 9141 92037755 BIAT 0258 9508 | 
2.0 0525 0146 .9575 8774 (8087 7117 .7603 0140 9426 | 12 0037 0001 0008 0.0001 0002 0071 
2.2 O82 O79 8473 8310 8502 461 7030 0075 9218) 4 .0007 0.0000 0.0000 0017 
24 0221 0042 9321 7768 8123 5807 6446 0039 8974 | 16 0.0001 ‘0004 
26 0142 0022 9114 ‘7175 .7508 5173 5865 0021 8097 | 18 0.0000 0001 
28 0OV1 0011 (8856 6554 7087 4571 5298 0011-8391 | 20 
22 
TABLE III. Contributions to Z, of the various electrons in helium. 
(le)? (la)(Qa) (2s) (2a) (2p) (2p)? (1a)(2p) | (le)? (le) (2a) (2s (2a)(2p) (2p)? (1e)(2p) 
r (la) (Ie) (2s) 2s)? (28) (2p) (2p) (le) (2p) | r (1a) (la) (2a) (2a)? (2s) (2p) (2p) (le) (2p) 
0.00 1.0000 a 0000 1 0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 3.0 0006 0000, .3332 : .1638 1913 1016 . 1333 0000 3282 
05 9161 9013 9850 9788 9804 .9778 9795 9007 .9873| 3.2 0004 2996 «1372 11629 0834 “1121 2921 
‘10 (8344 «(8057 9702 9580 9611 9555 9591 8046 9745) 3.4 0002 2682 (1144 «1381 0683 0041 2634 
‘15 7564 «7154 «9556 «9875 «9422 9333 «9386 «7138 9618) 3.6 0001 2392 10050 1167 0558 0788 2369 
20 G832 «6314 «9413 9176 «©9237 9110 9182 «6205 «9481 | 3.8 0001 2126 0786 ©0982 0454 «0858 2124 
40 0000 1882 0648 0823 0369 0549 “1900 
3 5522 4849 9137 8708 8883 8666 8773 4824 9236 | 
4 A423 3667 (S875 8442 8549 8223 8366 3640 .8082/ 4.4 1462 0437 0574 0243-0379 1511 
5 3520 2740 (8624 8105 8231 (7783 7960 2713 .8727| 48 "1122 0291 0395 10158 0260 “1191 
6 2787 2027 (8383 7783 7024 (7348 7559 2002 S473/ 5.2 0852 0192 0271 10108 0177 0933 
‘7 2198 [1488 (8150 7469 7625 6018 7162 1465 .8220| 5.6 0641 0126 0184 0066 0120 0724 
‘8 1728 (1084-7922 7160 ©7330 6498 6772-1065 7968) 6.0 0478 0082 0124 0042 0081 0559 
9 1856 (0786 :7609 6854 7037 6087 6300 0770 7717 
10 1061 0566 (7478 «6548 96745 5688 «= GOIS «0553-7467 | 6.8 0260 0034 0055 0017 0036 0327 
7.6 O138 OO14 O86 OOF 0016 0187 
12 0847 0201 .7039 5937 6161 4933 5306 0282 .6973| 84 0071 0006 0011 0002 0007 0105 
14 0393 O17 6608 5330 «=—.5580 0 ADM Ab 0142 6489 9.2 0036 .0003 0005 8.0001 0003 0058 
16 0237 0074 6166 A737 5010 (3617 4041 0071 6017 | 10.0 0018 0001 0002 0001 0032 
18 0143 0037 5731 4169 4460 3064 3495 0035 5561 
2.0 0086 0018 5300 3635 3938 2579 3007 0017 5121/12 0003 0000 0000 0007 
22 0051 0009 4878 3143 3452 2150 2575 0008 .4700| 14 0001 0001 
24 00831 0004 4466 2696 3005 1799 2196 0004 4300 | 16 
2.6 0018 0002 4070 2297 2509 1492 1865 0002 3921 | 18 
28 0011 0001 (3601 1945 2236 (1233 1579 0001 .3569 | 20 
' 
TABLE IV. One-electron energy parameiers. latter electron will not materially affect the wave 
———— Sa function of the (1s) electron and hence the 
€ CALL. BY e . ° ° ° . . 
Stare ELactron —2% € CALC. omen uemnons §=6. 1 Zation energy to a fair approximation will be 
(1)? (1s) 3.374 1.836 1.80% equal to the (2s) or (2p) value of «. We can then 
(1a) (2a) (3 — yA compare with the experimental values quoted by 
, at . = ; f > a _ ? ? 
on a sii eli ai Bac her and Goudsmit (reference 8a, p. 220). 
a. one one They give 0.3506 for the (1s)(2s) *S and 0.2921 
(22)(2p) ‘x “ yr . - _ . . 
\@p) naan 0.4048 for (1s)(2s) 'S; the latter is in fair agreement with 
' . ' , 
Ge (2p) ~— anaes our computed value. The experimental value for 
) 95 5; a " ' 
(14)(2p) {{39) aa by od (1s)(2p) *P is 0.2665 and that for (1s)(2p) 'P is 
= === 0.2478; the latter again agrees fairly well with 
*B. A. Hylleraas, Zeits. f. Physik 65, 209 (1990). . 
oF G. aerees P. + ved Pee Rev. 46. 7S (1934). our 0.2522. In each case the agreement appears to 
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Fic. 3. Radial charge density for helium. 


be better with the singlet than with the triplet 
state. One might indeed expect that the self- 
consistent field value would agree better with 
some average of the singlet and triplet energies. 
It must be emphasized again, however, that the 
« values are one-electron energy parameters only 
and further discussion of this question must 
await the calculation of the éofal energies of the 
various states. These are now under way and will 
be published shortly. It may be remarked that 
Hylleraas'® has recently computed by the general 
method of his previous papers many discrete 
energy levels for singly excited helium, such as 
both the singlet and triplet terms of (1s)(2s) and 
(1s)(2p). He neglects the effect of the excitation 
on the inner electron. The agreement of his re- 
sults with experiment is in general excellent. 


” E. A. Hylleraas, Zeits. f. Physik 83, 739 (1933). 
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Fic. 4. Total radial charge density for He(is)(2s). Solid 
line, Vinti; dotted line, present paper. 
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Fic. 5. Total radial charge density for He(1s)(2p). Solid 
line, present paper; dotted line, Eckart (1s)}(2p)*p; dash 
line, Eckart (1s)(2p)'p. 


In line with the procedure of some other 
writers, curves are presented (Figs. 2 and 3) 
exhibiting the radial charge density —dZ/dr as a 
function of r for the individual electrons as well 
as the completed group in the various states. 
These are self-explanatory and need no comment. 

A further comparison of the results of the 
present calculations with those of other methods 
is suggested in connection with the radial charge 
density. Using a variational method Vinti’ has 
found an analytic wave function for the (1s)(2s) 
1S state. After the function has been suitably 
normalized to unity the radial charge density for 
the two electrons may be computed by multiply- 
ing the integral 2f,°4rr,*| ~/|*dr_ (a function of 
r,=r) by 4xr*. The result is plotted in Fig. 4 and 
compared with —dZ/dr for the corresponding 
case of the present paper. The agreement is 
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rather good considering the difference in the 
methods. Fig. 5 shows a similar comparison be- 
tween the radial charge density for (1s)(2p) 'P 
and (1s)(2p) *P as computed by Eckart*® by the 
variational method and —dZ/dr for this state 
from the self-consistent field calculations. The 
agreement is again seen to be very gratifying. It is 


AND J. Y. 





BEACH 


interesting to note that the density for the singlet 
state of Eckart agrees better with the present 
values than does that for the triplet. 

The authors wish to express their gratitude to 
Dr. F. E. White for his kind assistance in connec- 
tion with the drawing of the curves and with 
some of the calculations. 
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The van der Waals Interaction of Hydrogen Atoms 


(Received March 18, 1935) 


The van der Waals interaction energy of two hydrogen atoms at large internuclear distances 
is discussed by the use of a linear variation function. By including in the variation function, in 
addition to the unperturbed wave function, 26 terms for the dipole-dipole interaction, 17 for 
the dipole-quadrupole interaction, and 26 for the quadrupole-quadrupole interaction, the 





interaction energy is evaluated as 
: 6.49903 e? 
Ww" = — 


ao p”* 


in which p 


124.399 e? 


ao p” 


adop ” 


R/ao, with R the internuclear distance. Some properties of the functions F,,,(£, 3, ¢), 


which are orthogonal for the volume element id sin 6d@dy, are discussed, and their usefulness 


in atomic problems is pointed out. 


INTRODUCTION 


N approximate second-order perturbation 
treatment of the inverse sixth power inter- 
action energy of two hydrogen atoms a large 
distance apart (corresponding to the so-called 
dipole-dipole van der Waals attraction) was given 
in 1930 by Eisenschitz and London.' This treat- 
ment led to the result W’’ = —e® A/aop", with 
p=14p/ dQ (faz being the internuclear distance for 
the two atoms), A being evaluated as 6.47. Ap- 
plications of the variation method by Hassé* and 
by Slater and Kirkwood* verified this result es- 
sentially, the constant A being shown to be 
equal to or greater than 6.4976. 
As early as 1927 this problem had been at- 
tacked by Wang,‘ using the method developed 
by Epstein® for the treatment of the Stark effect. 
1R. Eisenschitz and F. London, Zeits. f. Physik 60, 491 
(1930). 
*H. R. Hassé, Proc. Camb. Phil. Soc. 27, 66 (1931). 
<- C. Slater and J. G. Kirkwood, Phys. Rev. 37, 682 


(1931). 
*S. C. Wang, Physik. Zeits. 28, 663 (1927). 
+P. S. Epstein, Phys. Rev. 28, 695 (1926). 


Wang claimed to have obtained an exact solution; 
it was, however, pointed out by Eisenschitz and 
London that Wang’s result is necessarily in error. 
It seemed to us possible that Warty’s work might 
have contained only a numerical error, and that 
the method might actually be capable of giving 
an exact solution. Because of the usefulness 
which a method of exact solution of problems of 
this sort would have, we thought it worth while 
to study the problem thoroughly. We have found 
that the method used by Wang does not give an 
exact solution,’ but that it can be extended to 
give as closely approximate a solution as is de- 
sired. The treatment 
municated in this paper. 

A rough treatment of the dipole-quadrupole 
and quadrupole-quadrupole interactions of two 
hydrogen atoms has been published by Mar- 


results of the are com- 


* Wang does not present the final steps in his calculation 
in detail, but states that he set up a sixth degree secular 
equation from which he obtained an accurate value for the 
energy. We believe that the error in his treatment occurs 
at this point. 
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genau.’ We have applied our method fo obtain 
reasonably accurate expressions for these inter- 
actions also. 


FORMULATION OF THE PROBLEM 


In the wave equation for two hydrogen atoms 
Hy=Wy let us put H= H°®+H’, H® being the part 
of the Hamiltonian corresponding to two isolated 
hydrogen atoms and //’ representing the interac- 
tion of the two hydrogen atoms. In order to do 
this we neglect the resonance phenomenon (which 
is unimportant at large distances), taking as the 
unperturbed wave function the product function 
Vi00(A1) Wwo(B2); that is, we consider electron 1 
to be attached to nucleus A and electron 2 to 
nucleus B. We also write W= W°+ W” (W’, the 
first-order perturbation energy, being equal to 
zero). We shall consider only the interaction of 
two normal hydrogen atoms, so that W° is equal 
to —e*/as. Making the substitutions &= 2r4;/a9 
and £:= 2rg2/do, the wave equation becomes 


2 a 
(vit+ ss )v+ (vet - yy 
& 4 f& 4 


in which \= — W’’ao/ 2e’*. 
The interaction //’ is equal to 


—¢", rage, rate’, rapte’/Tip. 
If R(=r,z) is large, this can be expanded’ in in- 
verse powers of R to give the expression 
TT’ = (e?/ R*)(xyx2+ vive — 22,22) 


+(3/2)(e?/R* 17 1"Ze— 12") 


2+ 2(xyXe+ vive t+42i22)*} +---. (2) 


In this expression x,y,z; are Cartesian coordinates 
of electron 1 relative to nucleus A, and x2yez2 
those of electron 2 relative to nucleus B, the z 
axis for each being directed towards the other 
nucleus. The first term represents the mutual 
energy of two dipoles; this term alone is impor- 
tant for large values of R. The second and third 





7H. Margenau, Phys. Rev. 38, 747 (1931). 
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terms represent the dipole-quadrupole and 
quadrupole-quadrupole interactions, respectively. 
It can be easily shown that in the calculation of 
the second-order perturbation energy the terms 
can be considered separately, their contributions 
being additive. 

In the solution of the problem we shall make 
use of the functions Fy, (£, 3, ¢) discussed in the 
appendix. Each of these functions can be made 
identical with a hydrogen-like wave function 
Ynim (7, 8, ¢) by choosing a suitable linear rela- 
tion between r and £; the functions Fy, all con- 
tain the same exponential function, in contra- 
distinction to the functions ¥,im. We have defined 
£, and £ in such a way that Fi (:, 31, ¢:) and 
Fyo (£2, J2, 2) are identical with Yo (71, 31, ¢1) 
and yoo (72, Je, ¢2), respectively; that is, the un- 
perturbed wave function can be written as 
Fiw (é:) Fie (£2). We now apply the variation 
method in treating the perturbed wave equation, 
using as the variation function a linear combina- 
tion of the product functions F,..,,, (&, 1, ¢1) 
F,.d.0, (€2, 82, ¢2), with arbitrary coefficients. It 
can be seen that the second-order perturbation 
energy for the perturbation function 


IT’ = — (2e*/ao) aki &s Cos 3; COS De 


+ BE, ts" cos 3;(3 cos* dg—1) 


> 


+vé,7ko"(3 cos? 3; —1)(3 cos* d2—1) 
A ae 
in which 
a=(6)'ay?/8R®, 8 =(30)!at/32R4, 
and y=(70)'a,'/128R®, (4) 


is identical with that for the function of Eq. (2), 
and, moreover, that to obtain the first-order 
perturbed wave function and the second-order 
perturbation energy, the variation function used 
need contain, in addition to the unperturbed 
part Figo (£1) Fioo (€2), only the terms F,,1 (£1, 01) 
F,.w (£2, 82) (for the dipole-dipole term in a), 

: , I I 
F, 1 (&, 31) F,,20 (2, 82) (for the term in 8), and 
F,,2 (£1, 34). F,,20 (£3, 32) (for the term in y).* 

* On application of the ordinary methods of perturbation 
theory, it is seen that the first-order perturbed wave func- 
tion for a normal hydrogen atom with perturbation func- 
tion f(r)T(d, ¢), where T is a tesseral harmonic, has the 
form ¥,00(7) +7) T(d, ¢), the perturbed part involving 
the same tesseral harmonic as the perturbation function. 


The statements in the text can be verihed by an extension 
of this argument. 








_—— 
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For ‘a linear variation function }-,a,4, (where and 
#, represents the product functions Prd it , &.*O.dridre 
(&:, Hi, 1) Fyayu, (€2, 82, ¢2)) the secular equation i j ’ 
corresponding to the wave equation (1) is the 


determinantal equation® dr, being equal to &,* sin @déidd\dq, (with dr 


differing only in the subscripts) and the integrals 


HjP + Kix’ —Ajrdr| =O, (5) extending over the configuration space of the 

in which system. It is to be noted that the volume element 
drydtz is not such as to make the functions ® 

xunr= f far(or+(1 fi)—4 mutually orthogonal. The integrals can be 
. evaluated with the help of the relations given in 


+9?+(1/t)—4})didridre, the appendix. To obtain the second-order per- 

turbation energy we need introduce the term in A 

Ia! -{{ *( —(ao/2e*)H’)@id rd 72, in the row and column corresponding to the wave 
function for the unperturbed system only. 





THe Dipo_te-DieoLte INTERACTION 
The secular equation for the dipole-dipole interaction is 


») Ay mhadeee 


100100 —4 32a — l6a — 16a 8a 0 0 0 0 0 
210210 32a -8 2 2 0 0 0 0 0 0 

210310 — 16a 2 —14 0 4 v¥(10) 0 0 0 0 

310210 — l6a 2 0 —14 4 0 y¥y(10) 0 0 0 

310310 8a 0 4 4 —24 0 0 2V(10) 2¥(10) 0 

210410 0 0 v¥(10) 0 0 — 30 0 6 0 O «+. ] 

410210 0 0 0 v (10) 0 0 —20 0 6 0 ---!] =0, (6) 
310410 0 0 0 0 2¥(10 6 0 —34 0 3¥v(10 

410310 0 0 0 0 2y¥(10) 0 6 0 —34 3y(10) : 

410410 0 0 0 0 0 0 0 3v(10) 3V¥(10) —48 . | 





the rows and columns corresponding to the values of »+ - - 4: shown at the left. We obtain successive 
approximations to the solution of this equation by neglecting rows and columns beyond the nth. This 
process has been carried out for n= 2, 5, 10, 17 and 26. Some simplification is achieved by combining 
rows (and columns) with » and v; interchanged, the corresponding functions having the same 
coefficient. The results of the calculation are given in Table I, in terms of the constant A in the ex- 





TABLE I. The dipole-dipole TABLE II. The dipole-quadrupole TABLE III. The quadrupole-quad- 
interaction constant A. interaction constant B. rupole interaction constant C. 
Degreeof Terms Degreeof Terms Degreeof Terms 
approx. included A approx. included B approx. _includec Cc 
2 2, v2 6 2 yy 2, » 3 115.7 2 y, 3, v3 1063.1 
3 3 6.4822 5 3 4 124.10 5 4 4 1132.6 
10 + 4 6.4984 10 4 5 124.386 10 3 5 1134.35 
17 5 5 6.49899 17 5 6 124.399 17 6 6 1135.12 
26 6 6 6.49903 26 7 7 1135.21 











pression W” = — Ae*/agp*, with p= R/ay (A being equal to 3X/ 16a”). 
It is seen that the convergence is rapid, the final value of A, 6.49903, 
one unit in the last decimal place. 


being trustworthy to within 


*See, for example, L. Pauling and E. B. Wilson, Jr., Jntroduction to Quantum Mechanics, with Applications to 


Chemistry, McGraw-Hill Book Co., 1935, Chap. VII. 
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THE D1IPOLE-QUADRUPOLE INTERACTION 
The secular equation for the dipole-quadrupole interaction is 
My AP aAgue 
100100 —4 —192V(3)8 96V(3)8 96v28 —48v278 0 0 0 0 0 
210320 —192¥(3)8 —14 4 v (6) 0 0 0 0 0 0 
310320 967 (3)8 4 —24 0 2¥(6) 2¥(10) 0 0 0 0 
210420 9628 v(6) 0 —2Ww 6 0 y(14) 0 0 0 
310420 — 4828 0 2v(6) 6 —34 0 0 3¥(10) 2¥(14) 0 
410320 0 0 2yv(10) 0 0 -34 0 3v(6) 0 0 
210520 0 0 0 y14) 0 0 —26 0 8 0 =0. (7) 
410420 0 0 0 0 3v(10) 3y¥(6) 0 —48 0 3y(14) 
310520 0 0 0 0 2¥(14) 0 8 0 —44 4y(10) 
410520 0 0 0 0 0 0 0 3Vv(14) 4v¥(10) —62 





We have solved this equation approximately, the results being given in Table II in terms of the 
constant B in the energy expression — Be*/ao¢’. 
The error in the final value of B we estimate to be less than one unit in the last figure quoted. 


THE QUADRUPOLE-QUADRUPOLE INTERACTION 


The secular equation for the quadrupole-quadrupole interaction is 











vy Amira 

100100 —4r 34567 —576V(6)y —S76v(6)y 576y 0 0 0 0 0 
320320 34567 —24 2v(6) 2V(6) 0 0 0 0 0 0 
320420 S76V(6)y7 2v(6 —34 0 3v¥(6) 2y(14) 0 0 0 0 
420320 —576V(6)y 2v(6 0 34 3v(6) 0 2y(14) 0 0 0 
420420 5764 0 3v(6 3¥v(6) —48 0 0 3¥V(14) 3v(14) 0 
320520 0 0 2v(14 0 0 —44 0 4y(6) 0 0 
520320 0 0 0 2v(14 0 0 —44 0 4y(6) 0 =(). (8) 
420520 0 0 0 0 3vV(14) 4¥v(6) 0 —62 0 4y(14) 
520420 0 0 0 0 3v(14 0 4y¥(6) 0 -62 4y(14) 
520520 0 0 0 0 0 0 0 4v(14) 4y¥(14) —80 





The results of the approximate solution of this equation, in terms of the constant C in the energy 
expression — Ce*/ap'*, are given in Table III. 
The final value, C= 1135.21, is reliable except for the last figure. 


DISCUSSION OF RESULTS 


We have thus found for the interaction energy of two normal hydrogen atoms at the large distance 
R= pd the expression 
6.49903 e? 124.399e? 1135.21 2 
2? oe cm namningdnnnnns cs emminliagiintae chincimmapeanin ie ¢ 6 


Aop® 10 


(9) 


It is interesting to note that the value A = 6.4976 found by Hassé for the dipole-dipole coefficient by 
the use of a variation function yYyoo (7:1) vo (72) [1+H'(A+Brire+ Crfr?+Drr?)} is very close to 
our value for »,=4, »:=4, which is based on a variation function involving all terms (unsymmetric 
as well as symmetric) out to 7,*r,*. This indicates that the unsymmetric terms are of minor importance. 


© We use here shortened symbols, such as F,, F,, for Fy, r.5:(, O:, 1) Peed-ael Es, Oa, @2), etc. 
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The approximate second-order perturbation energy Wo'’= (7/')?/ W,°, with Wo°= —e?/ao, leads to 
the values A=6, B=135, and C= 1417.5, the last two values being given by Margenau.’ It is seen 
that the value of A is too low, and those of B and C are too high. This means that the dipole-dipole 
interaction is due more to excited states with negative energy (less than e®/a) above the normal state 
of the system of two hydrogen atoms) than to excited states with positive energy, whereas the dipole- 
quadrupole and quadrupole-quadrupole interactions are due more to the latter than to the former 


states. 
As has been pointed out by earlier authors, the van der Waals forces are more important than ex- 


change forces for values of R greater than about 7a». At this distance the dipole-quadrupole force is 
about one-half as large as the dipole-dipole force, and the quadrupole-quadrupole force is about 
one-eighth as large, the dipole-dipole attraction becoming relatively still more important at larger 
distances. This van der Waals calculation, based on product wave functions, is not significant for 
values of R much less than 7a, because of neglect of the resonance phenomenon and because of failure 


of the expansion given in Eq. (2). 

While our treatment has not led to an exact solution of our problem, the use of the functions 
Fy, (&, 3, ¢) has permitted the reasonably accurate approximate solution to be made with considerable 
ease, and we feel that these functions may be found useful in the treatment of other problems of atomic 


and molecular structure." 


APPENDIX 
The functions F,,, (£, 3, ¢) 


We define 


F,,,(é, 8, ¢) =Ayn(£)P,,(8)®,(¢), 
with 
Aya(é) = {(v—A—1) !/[(e +A) EF 
L being an associated Laguerre polynomial, 


@,,(8) = {[(2A+-1)(a —p)!]/2(A+p) !} 4P,!4!(cos 3), 


P being an associated Legendre function, and 


®.(¢)=(2r ie 


The functions are orthogonal and normalized for the weight function £ sin J, satisfying the equa- 
tion 


<* 


aco eo o-F 
F’ x(k, 8, ¢) Fa, (€, 8, e)é sin Od gdddt = 6,-,6)-,4,°, 
6 a « to» 


Fy,(é, 8, ¢) satishes the differential equation 


(V24+1/&—1)F,.,= —[(»—1)/8]Froap. 


The following relations involving the A's can be easily derived from the properties of the associated 
Laguerre polynomials: 
€A,. = — | (v—A) (wo +A41)}4A,1 »+2vA,,— | (¥+A)(v—-A-1 ;7A, L- 3 


> 


EFA = — | (h—A)( oP +A41)(P4+A42)(P4+A43)} Aye vc) 
+2(2»—A)}(r+A+1 (y+A+2 i ?ALoy .21— Or, y+A+1 y—A-—1 ; 7A, hal 


+2(2>»+A)}(»—A—1)(w—A—2)} PAL, agi — 1 (HP +A) (e—A—1)(v—A—2)(v—A—3)} PAL, a5: 


“\ Added in proof: Professor ]. H. Van Vieck has pointed out to us that the functions were used in 
dispersion by hydrogen-like atoms by B. Podolsky, Proc. Nat. Acad. Sci. 14, 253 (1928) 
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7A = — | (v—A) (vo +A4-1)(O+A42)(0+A43) (0 +A4+4)(9+A45)} Aras, ace 
+2(3»—2d){(v+A+1)(v+A42) (9 +A43)(¥+A4+4)] Arie, a2 
—5(3v—X) | (v—A—1)(v+A4+1)(¥+A42)(V4+A4+3)} Arar. vce 
+20v{(v+A4+1)(v+A4+2)(v—A—1)(v—A—2)} FAY, vue 
—5(3v+ A) {(v+A4+1)(v—A—1)(v—A—2)(v—A—3)} PAL nae 
+2(3v+2Xd)}(v—A—1)(v—A—2)(v—A—3)(v—A—A4) J PAL e vce 
— {(v+A)(v—A—1)(v—A—2)(v—A—3)(v—A—4)(W—A—5S)} MALes, vge. 


From these we obtain similar relations in the F’s. The following special cases are needed in evaluating 
the matrix elements in Eqs. (6), (7) and (8): 
EF i00 = V2 F200 +2 F yo, 
£ cos J Fi99 = — 2V2 F319 +42 Faro, 


£3/ 3 cos *8— 1) Fro = —24 F gag + 24 (6) Feo. 


The evaluation of the matrix elements 


In order to illustrate the method of construction of the matrices in Eqs. (6), (7) and (8), we shall 


evaluate some of the integrals. Let us first consider the integral'® 


x, nna { {Pe Fe, {0 +1 


By using the differential equation for Fy», (£, J, ¢), the equation becomes 


i—-434V2°4 1 E> ANF, F dr ydry. 


wre 


KK», v2 "ie > 7s F*, a = . hee 1 aie | Es; F, F, £; E> sin Uv; sin Jed Edd id gid Exdd ed 2 
= -oe F*,.°F*, ‘i —( Py | fo—(ve- | £14 Fy Fskiks sin GO; sin Ged Edd 1d gd odd ad G2. 
For Eq. 6) (A, =A2=1, uw: = 42=0) this becomes 
K | F* F*,.'T (ve—1 vi —1)(», +2) PF, Fi.4+(»,-1 vo —1)(ve4+2)} Foy Fett 


—2 2vys Vin Po F, F., T (Ve 1); vit 1 vy,—2 \tF, - Fos 


+-(v,—1 vot1)(v2—2)}* Fe, Fe, wits sin J, sin Jedi\dd de didddd 
and hence we obtain 

HK ’ ve—1 b 1)(y,+2 4 

C°,;, » 1, y =(7,—1 vo—1)(v2+2 

x’, oo = (ve—1)} (9, +1)(%,—2)] 4, 

K no, a =(¥—1 ve+1)(v2—2)}', 
and 

HK 2 23 P= vi = Ve , 


all others being equal to zero. 
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In Eq. (7) we have 


AND 
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Xe owe. 1 -6- . | [ Pot Pts cos ,(3 cos? v> = 1) Fy Fy é,7&s? sin dv; sin Pod Edd 1d 1d fed dod ge 


-sf- . | Fy,* Frg* | 48V2 F310 F 420 — 96V (3) F 3,9 F 320 — 96V2 F2i0F s20 


+192y (3) FeoF s20} E:é2 sin 3; sin Podé dd 1d gid Eedd ad go, 


and hence 
KH’ o3,11 = 192 y (3)8, 


and %’s4, 1. =48v28, all others being equal to zero. 


K’ 33,11 = — 96 (3)8, 


KH’ 24,11 = — 96V28, 


To illustrate the evaluation of the integral A let us consider Ay y: 


An, u=f- . fi F yoo( £1) Froo( Ea) | *E12E2? sin Jy sin Ped E,dd dg id Eedded g2 


-{- ° + {Fao &s) Fron) {2 Faoo( &1) Feoo( $2) — 2V2 Feool 1) Fro0(€2) — 2V2 Fioo( €1) Feoo( E2) 
+4 Fyo0( £1) Fioo( 2) } Eré2 sin 3; sin Pedi,dd dy diddade., 


or, making use of the orthogonality and normalization of the F's, Au, n»=4 
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The Magnetic Moment of the K” Nucleus 


J. J. Gippons, Jr. AND J. H. Bartvett, Jr., Department of Physics, University of Illinois 
(Received March 16, 1935 


Starting with the K* field given by Hartree, a 4s wave 
function has been found by numerical integration. The 
orthogonality correction diminishes the value of the wave 
function at the origin by about one-third. The experimental 
value of the hyperfine structure separation of the 2S normal 
state is 0.015 cm™. From this is calculated a molecular 


HERE have recently been published two 
independent determinations of the h.f.s. 
separation of the normal state of the K*® atom. 
Jackson and Kuhn! have used a spectroscopic 
method ; Millman, Fox and Rabi’? have relied on 
the method of molecular beams. The numerical 
agreement is good, and 0.015 cm! may be taken 
as the value of the separation. 
From this, Millman, Fox and Rabi estimate 
the value of the nuclear magnetic moment to be 


1D. A. Jackson and H. Kuhn, Nature 134, 25 (1934); 
Proc. Roy. Soc. A148, 335 (1935). 
? Millman, Fox and Rabi, Phys. Rev. 46, 320 (1934). 


These authors were able to determine the spin as well as 
the h.f.s. splitting. 


magnetic moment of 1.2 nuclear magnetons, as compared 
with the value of 0.38 nuclear magnetons calculated by 
Millman, Fox and Rabi on the basis of the modified 
Goudsmit’s formula. The disagreement seems to be due 
chiefly to the possibility that the single electron wave 
functions of Fermi and Segré are not mutually orthogonal 


0.38 nuclear magneton, using the semi-empirical 
formula due to Goudsmit* and Fermi and 
Segré.*: > Since previous cases of low magnetic 
moments have been associated with isotopes of 
even mass number, the seemed 
somewhat anomalous, and we therefore thought 
it of interest to ascertain whether or not the 
use of a 


above value 


“anomaly” could be removed by the 
Hartree wave function. This we now proceed to 


investigate. 
*S. Goudsmit, Phys. Rev. 43, 636 (1933 
*E. Fermi and E. Segré, Zeits. f. Physik 82, 729 (1933). 
* E. Fermi and E. Segré, Memorie, R. Accademia d'Italia 
4. 131 (1933). 











MAGNETIC MOMENT 
WAVE FUNCTIONS 


Starting with the field of K*, which has been 
published by Hartree,* we have performed a 
numerical integration to obtain the 4s wave 
function.’ The results are given in Table I. 

For simplification of the calculations, it is 
desirable to make the wave functions mutu- 
ally orthogonal. The orthogonalization integrals 
are: (1s, 2s)=0.0113; (1s, 3s)=0.00365; (1s, 4s) 
= (0.000562; (2s, 3s)=0.0212; (2s, 4s)= 0.00422; 
and (3s, 4s)=0.0506. Let us denote by P® 
the functions after orthogonalization and _ nor- 
malization (f(P°)*dr=1). Then P®(2s)=P(2s) 
—0.0113 P(is); P%3s)=P(3s)—0.0212 P(2s) 
—0.00341 P(is); and P%(4s)=1.0013[ P(4s) 
— 0.0505 P(3s) —0.00315 P(2s) — 0.000342 P(1s) }. 
The values of the nonorthogonal wave functions 
x=(P/r) at the origin are: x,,(0)= 161.9; x2,(0) 
= 45.17; x3,(0)=15.76; and x,4,(0)=3.124. Be- 
fore orthogonalization, [x4,(0) P=9.76; after- 
wards, [x4.°(0) P= 4.548. (That is, the value of 
the nuclear magnetic moment would be smaller 
by a factor of about two if the calculations were 
made with the nonorthogonal wave functions.) 


NUCLEAR MOMENT 


The h.f.s. splitting of an s-state is given® by 
A(s)=(8xr/3)[(27+1)/I ]upo¥*(0), where J is the 
nuclear spin, uo is the Bohr magneton, u is the 
nuclear magnetic moment, and ¥(0) is the value 
of the wave function (normalized so _ that 
f¥dr=1) at the nucleus. We have® (0) 
= (1/4r)'x(0), so that A(s)=(2/3)[(27+1)/J]z 
X wox?(0). For K**, 7=3/2. Substituting numer- 
ical values, we have, if g= (w/e) X 1838, 


(2/3)(8/3)(g/1838)(9.17 x 10-*")*4.548 





0.015 = - — . 
(0.528 & 10-*)* x 19.662 « 10-"” 


Solving, g=949A(s)[J/(27+1) J[1/x7(0) J=1.17. 
That is, we obtain a nuclear magnetic moment 
about three times as large as that found by 
Fermi. This disagreement has its origin, of course, 
in the fact that different eigenfunctions have 


*D. R. Hartree, Proc. Roy. Soc. A143, 506 (1934). 

7 As usual, it has been assumed that it is unnecessary 
to make the field self-consistent. 

*E. Fermi, Zeits. f. Physik 60, 322 (1930). 

* The value of the nuclear moment reported by us at the 
Pittsburgh meeting of the Am. Phys. Soc. was too small 
by this factor 4x, which had been omitted. 
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Taste I. Normalized P(4s) for potassium. 











~ « =0.265 

r P E r P r P ” 
0.000 §6(0.0000 19.000 0.40 —0.,0598 60 0.4220 0.1199 
005 142 19.374 45 — 0312 65 — 3937 .1573 
010 0257 19.788 50 — .0003 70 — 3609 .1881 
O15 0349 20.254 55 + .0307 75 — 3264 2199 
.020 0421 20.779 60 + OOO1 8.0 — .2916 .2358 
8&5 — .2878 .2547 
03 0509 22.048 7 .1105 90 — 2261 .2712 
.04 0540 23.724 8 .1465 95 — .1967 .2857 
OS 0527 26.001 x 1679 10.0 — .1700 -2988 

.06 .0480 1.0 1761 

07 .0409 1.1 1734 i! —- .1246 .3206 
08 0321 1.2 1623 i2 — 0896 .3383 
09 -0220 13 — 0634 3531 
10 0113 14 1226 14 — 0443 3656 
1.6 iS — 0305 .3763 
12 — .0108 18 0119 16 — 0209 3856 
14 —.0321 2.0 — 0470 17 — 0138 .5937 
16 —.0511 2.2 — .1040 18 — 0095 .4008 
18 — 0673 24 — .1578 19 — 0063 4071 
20 —.0801 2.6 — .2074 20 — 042 4128 
22 — .0896 2.8 — .2524 21 — .0028 4179 
24 — .0958 3.0 — .2926 22 — O18 4225 
26 — .0989 3.2 — .3279 23 — O12 4267 
28 — .0993 34 — .3583 24 — .0008 4305 
30 — 0972 3.6 — 3840 25 — 0005 4340 
38 — .4052 26 — 003 4372 
35 — 0830 40 — 4221 27 — 0002 A402 
28 — 0001 4429 
4.5 — 4476 29 — 0001 4455 
5.0 — 4530 »w — 0000 4480 
5.5 4430 31 — 0000 4502 








been employed. We use a Hartree function; 
Fermi and Segré use a one-electron function 
based on the statistical method, which function 
contains parameters to be determined from the 
experimental term values 

Since the article of Fermi and Segré® is some- 
what inaccessible, we shall indicate their pro- 
cedure.'® Let the Schrédinger radial function be 
p¥, where p is the normalization factor. Then 
u=ry satisfies the differential equation (d*u/dr*) 
+ (82*m/h*)(E—eV)u=0. The potential V is 
found by the statistical method, and the equation 
solved" for E=0 by the method of Wentzel- 
Brillouin. One sets u=&R sin 6, where R*(dé/dr 
=1. k may be determined by matching such a 
solution with the power series development 
about the origin. 

To evaluate p, it is supposed that at a point 
t a perturbation \é(r—£) is introduced, where \ 
is an arbitrary parameter and 4 is the Dirac 
improper function. The first-order perturbation 
energy is then AE=4r\p*u*(t). At the point &, 
the value of u will not be changed, but there will 


1° Certain details have been omitted by Fermi and Segré. 
These we have tried to reconstruct as accurately as 
possible. 

“ In the article by F. & S., the term “«s eigenfunctions” 
is used. We interpret this to mean s eigenfunctions with 
E=0. Such eigenfunctions, which were used earlier by 
Fermi in his calculations of the Rydberg corrections, do 
not differ much from the optical ones in the interior of 
the atom. ° 
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be a phase change A@. If this were 7, the number 
of nodes would be different by unity. Therefore, 
put AE= —(Aé@/r)(dE/dn). The discontinuity 
at £ introduced by the perturbation into the 
first derivative of the wave function is Aw’ 
= (8x*m/h*)\u(t). Now u’=kR’ sin 0+kR cos 
@x (d0/dr). Therefore 


Au’=kR sin 0A(R’/R)+(k/R sin 6)A(sin 6 cos @) 
=kR sin 0A(R’/R)+-(k/2R sin @)A(sin 26) 
=kRR’ sin 6A(1/R)+(k&/R sin @) cos 2040 


provided that R’ is not appreciably altered.” 


Since Au=0, we have (R+AR) sin (6+A8) 
=Rsin@, or AR=—Rcot @A@. Accordingly, 
A(1/R)=(1/R) cot @A@, and Au’=kR’ cos 6A6 


+(k/R sin 6) cos 2040. Now let us choose" the 
point £ so that @=27/2. Then" (48@),2= 
— (8x*m/h*)\R*. Substituting in the expression 
for AE, we finally obtain p?= (2m/h*k*)(dE/dn). 
The value of & was determined for 15 atoms by 
Amaldi and Fermi," and correction factors with 
respect to the formula ¥(0)=(1/ra*)(Z/2Rh) 
x (dE/dn) were calculated for various values of 
Z. For small Z, the factor is approximately unity. 
If, then, we have a Rydberg series, we can 
readily compute ¥7(0) from the empirical value 


2 This will probably be the case, since R is the amplitude 
of the oscillating function. 

“ This brings up the question as to whether the pro- 
cedure is actually independent of the choice of £. Fermi 
and Segré do not discuss the point. 

“This is in agreement with formula (30), p. 152, 
reference 5, for @=x/2, but we are unable to check the 
general validity of this formula. 

® Reference 30 (reference 5). 
still unpublished. 


Apparently, this work is 


AND 
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of dE/dn, so that this formula is quite useful. 

The calculations of Fermi and Segré seem to 
us somewhat incomplete, however, since they 
have not found a set of mutually orthogonal 
single-electron wave functions. This is, appar- 
ently, the principal reason for the disagreement 
between our results and theirs, as may be seen 
by a comparison of the nonorthogonal 4s 
functions calculated by the two methods. The 
statistical method gives x,,"(0)= 13.7; while the 
Hartree method results in x,,2(0)=9.76. The 
agreement is, in our opinion, as good as is to be 
expected. 

We do not at present feel able to estimate how 
accurate the results from the Hartree method 
are. Experience to date seems to indicate that 
the Hartree functions are fairly reliable, espe- 
cially when applied to scattering problems. We 
conclude that the value of the magnetic moment 
found from the Hartree method is probably the 
most accurate so far, but may be subject to 
modification when better functions are 
known. Because of this uncertainty, we have not 
calculated a value of the magnetic moment from 
the splitting of the *P 3,2 state, for even if agree- 
ment with the value from the *S,,2 state were 
obtained, it might not be overly significant." 

Finally, our value of 1.2 nuclear magnetons is 
not as abnormally low as that found by Millman, 
Fox and Rabi, so that we do not consider that 
there is any basis for regarding the structure of 


wave 


the K** nucleus as exceptional. 


The *P state is probably perturbed, according 
Fermi and Segré (reference 5, p. 141). 





i to 
hey 


nal 


ar- 
ent 
-en 


4s 


he 
he 
be 


Ww 
od 
at 
e - 
Ve 
nt 
he 
to 








MAY 1, 1935 


PHYSICAL REVIEW 


VOLUME 47 


The Dielectric Constants of Air and Hydrogen at High Pressures 


RALPH MCNABNEY; WELLS MouLTON AND W. L. BeuscHLern, Department of Chemical Engineering, 


University of Washington 
(Received January 8, 1935) 


The dielectric constants of air and hydrogen at pressures 
from 71.8 to 334.7 atmospheres at 20°C have been meas- 
ured. The values were obtained by measuring, with a 
capacity meter of the heterodyne type, the variation with 
pressure of the capacity of a fixed condenser fitted into 
the pressure vessel. The frequency was 2500 kilocycles. 
The dielectric constants obtained are slightly lower than 


those of Tang! and Broxon for pressures about 80 atmos- 
pheres. The Clausius-Mosotti function («—1)/d(e+2) was 
calculated for each pressure. The change with pressure for 
air is small and irregular, average 1.43. For hydrogen the 
value decreased uniformly from 1.16 at 71.8 atmospheres 
to 0.99 at 334.7 atmospheres. 





INTRODUCTION 


HE development of a differential manometer 

for use with hydrogen and air over the pres- 
sure range 1 to 340 atmospheres in which the de- 
flection of a diaphragm changed the capacity 
between attached condenser plates, required 
that the dielectric constant of the gases be known. 
Measurements of the dielectric constant of air 
have been made to the pressure of 170 atmos- 
pheres; while for hydrogen, some work has been 
done on pressures to 194 atmospheres. The di- 
electric constant of air at pressures up to 170 
atmospheres has been determined at low fre- 
quencies or by static methods by Tangl,':* 
Waibel® and Broxon.* Jordan, Broxon and Walz® 
have made measurements of the dielectric con- 
stant of air at pressures up to 170 atmospheres 
and at frequencies of 1 to 70 kilocycles. Oc- 
chialini,® using a frequency of about 10,000 
cycles, obtained results which were lower than 
those of the other workers mentioned here. 


DESCRIPTION OF APPARATUS 


A fixed air condenser was made to fit into a 
pressure vessel and variation of capacity of this 
condenser with pressure was measured with a 
capacity meter of the heterodyne type, operating 
at a frequency of 2500 kilocycles. 

The gas, from a Rix compressor, was cooled to 


' Karl Tangl, Ann. d. Physik 23, 559 (1907). 
? Karl Tangl, Ann. d. Physik 26, 59 (1908 
F. Waibel, Ann. d. Physik 72, 161 (1923). 
* James W. Broxon, Phys. Rev. 37, 1338 (1931). 
*A.R. Jordan, J. W. Broxon and F. C. Walz, Phys. Rev. 
46. 66 (1934 
* A. Occhialini, Physik. Zeits. 6, 669 (1905). 


about 15°C and passed through a baffle-plate 
type separator and then through a filter consist- 
ing of six inches of tightly packed felt. This ar- 
rangement removed all of the oil and lowered the 
water concentration to 0.015 percent. Moisture 
determinations were made by drying samples 
over Desic-chora. After the filter, the gas flowed 
through a coil of 1.5-millimeter diameter copper 
tubing where it was brought to room temperature, 
which was kept constant (+1°C) throughout the 
entire series of measurements. A cylindrical steel 
pressure vessel, having inside dimensions of 9 
cm diameter by 6 cm, was connected through 
appropriate valves to the tubing. A Bourdon 
type gauge which, for this work, had been cali- 
brated on a dead weight tester was in the line. 
The accuracy of the pressure measurements was 
at least one percent. 

The condenser shown in Fig. 1 was fastened to 
the bottom of the pressure vessel and the elec- 
trical lead to it carried through an isolantite 
bushing in the top of the steel shell. A spring at- 
tached to the bushing made contact with the top 
plate of the condenser. 
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Fic. 1. Condenser construction. 
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In order to make correction for the effect of 
solid insulation, the condenser was built so that 
the parallel condenser plates, upon which the 
change in capacity with pressure was deter- 
mined, were isolated from the solid insulation by 
means of shielding and could be removed from 
the frame of the condenser and replaced without 
change in capacity. This was accomplished by 
mounting the condenser plates C on end plates 
B and D which were in turn fastened to the frame 
of the condenser in such a way that the entire 
assembly of condenser and end plates could be 
removed from the frame as a unit. The top ring 
of the frame was insulated from the bottom 
plate by four isolantite posts E which were 
shielded from the condenser plates by the cylinder 
F fastened to the bottom plate and separated 
from the top plate by a 0.3-millimeter gap. 


Capacity meter 

The high frequency heterodyne type capacity 
meter, shown in Fig. 2, was used in which the 
unknown capacity acted as the tuning condenser 
of the “manometer oscillator’? which employed 
a modified Hartley circuit. The plate end of the 
oscillatory circuit was grounded and current was 
carried to the tube filament through the double- 
wound coil Ly-L;. This circuit required the use 
of only one bushing through the shell of the 
pressure vessel and allowed the base of the con- 
denser and the shielding to be fastened directly 
to the grounded shell. 

Another oscillator, which also acted as a de- 
tector, was coupled to the manometer oscillator 
by means of a low impedance link circuit. This 
consisted of single loops coupled to the oscillator 
coils connected by a twisted pair, one conductor 
of which was grounded. The standard Hartley 
circuit was used in this detector unit. An ac- 
curate vernier dial on the tuning condenser C, 
gave a maximum accuracy of one part in one 
thousand. C; is a trimming condenser used to 
adjust the range of the capacity meter. Both 
units were shielded by means of heavy copper 
boxes. 

The beat note between the oscillators could be 
heard in the telephone receivers connected in the 
plate circuit of the detector-oscillator unit. Zero 
beat, the setting where the frequencies of the 
two oscillators were exactly equal, was the point 
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Fic. 2. Circuit for the high frequency heterodyne capacity 
meter. 


at which readings were made. The use of a 
moderately high frequency, two to three mega- 
cycles, the stability of the oscillators and the 
very loose coupling between them minimized 
the effect of ‘pulling in” to such an extent that 
the zero beat method could be used without 
difficulty. Since both the oscillators were con- 
nected to the same power supply the factors 
which would cause the frequency of an oscillator 
to drift would cause both to drift the same direc- 
tion by the same amount. This would not affect 
the frequency differences and consequently 
changes in plate and filament potential and 
variation in temperature had little effect on the 
calibration. This apparatus was built rigidly 
enough so that there could not be any change of 
its calibration because of mechanical distortion 
of its component parts. 


EXPERIMENTAL PROCEDURE 


The capacity meter was calibrated against a 
General Radio Type 222 precision condenser. 
This condenser was substituted for the unknown 
capacity and readings made at frequent intervals 
on the scale. Calibrations were made after each 
set of dielectric measurements and the maximum 
difference between corresponding points on the 
various calibrations was found to be less than 
0.2 percent. As a further check on the calibration, 
a second variable condenser was placed across the 
standard condenser and the lower end of the 
standard condenser scale checked against other 
sections by means of the capacity meter. 

Air or hydrogen at 340 atmospheres was 
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slowly released into the pressure vessel which 
contained the fixed condenser. At intervals of 
about 70 atmospheres, the inflow of gas was 
stopped and the capacity of the condenser was 
measured. The system was allowed to stand at 
each pressure until the capacity remained con- 
stant, although the change that did occur at first 
was negligible. After the maximum pressure had 
been reached, the measurements were repeated 
with decreasing pressure. 

The pressure vessel was then opened and the 
plate assembly removed from the condenser. A 
blank cover plate of exactly the same dimensions 
as the top plate B Fig. 1, was put on the frame 
and the pressure vessel again closed. The gas was 
run in and capacity measurements made for each 
pressure in the manner previously described. 
These measurements gave the variation with 
pressure of the capacity of the condenser frame, 
solid insulation, contact spring, and shielding. 

When the condenser is assembled with the 
plates in, it may be considered to consist of the 
capacity between the plates, the dielectric of 
which is entirely air, in parallel with the capacity 
of the shielding, frame, insulation, etc., some of 
which are solid dielectrics. The dimensions of the 
exterior of the condenser were exactly the same 
whether the blank cover plate or the plate B 
were in, so that the capacity between the con- 
denser and pressure vessel was the same in both 
cases. Therefore, the difference between the 
capacities with the plates C in and with them 
out represented the capacity of these plates with 
each other and with the surrounding shielding. 
No solid insulation was in the field of this capac- 
ity. Consequently the variation, with change in 
pressure, of this difference Co—co where, 

Co=capacity at 1 atmosphere with plate C in, 
Co=cCapacity at 1 atmosphere with plate C out, 
was a capacity on which the solid dielectric of 

the condenser had no effect. 

One of the difficulties encountered in this cir- 
cuit was the change in reactance of some of its 
components with the change in capacity of the 
condenser inside the pressure vessel. The worst 
offender in this respect was the bushing through 
the pressure vessel shell. The value of Cy—co 
was measured with the condenser outside of the 
pressure vessel and found to be 2.0 mmf less 
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Fic. 3. Dielectric constant of air and hydrogen asa function 
of pressure. 


than the value of Cy—c» for the condenser when 
in the pressure vessel. A variable condenser was 
placed across the terminals of the capacity meter, 
the ungrounded lead of which ran through the 
bushing. The apparent capacity of the bushing 
was then measured for various settings of the 
condenser. Between the capacity settings cor- 
responding to co and Co, respectively, the capacity 
of the bushing was found to increase 2.0 mmf. 
This accounts for the discrepancy in the values 
of Cy—co. Correction was made for this variation 
in the results. The error introduced by the 
difference in resistance between the standard and 
substituted condenser when the bushing was not 
in the circuit was in the order of 0.005 percent, 
which could be neglected. The mechanical 
strength required of a bushing to withstand such 
high pressures requires a type of construction in 
which the electrical losses are high and the 
reactance variable. 

In the dielectric constant determination of 
air the accuracy of (e—1) was about one percent 
of the maximum value while for hydrogen the 
accuracy on this basis was about two percent. 


TABLE I. Dielectric constants of air and hydrogen. Frequency 
2.5 megacycles. Temperature 20°C. 


- AR — «e-1 HYDROGEN— «e-! 
DIgELEcTRIC \- 73)* DIeLecTric + 73) . 
PRESS CONSTANT (¢) i CON NT (« 1 
Amos) Runt Rux2 (3)! Rows Rum2 (3)! 
72.0 1.039 1.039 1.47 1.019 1.021 1.16 
140.0 1.073 1.071 1.40 1.032 1.036 1.07 
208.0 1.106 1.105 1.41 1.048 1.050 1.05 
274.0 1.134 1.134 1.42 1.059 1.061 1.00 
335.0 1.162 1.159 1.45 1.070 1.070 0.99 
* Values of d were calculated from L.C.T. 
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Since the dielectric constant of air at one at- 
mosphere (¢,) is 1.00059, the dielectric constant 
at pressure P is given by the equation, 


ép= [(Cp—cp)/(Co—co) }(e:). 


The values of the dielectric constant deter- 
mined for the various pressures are given in 
Table I and are shown graphically in Fig. 3. 
The measurements for air obtained with ascend- 
ing and descending pressures checked exactly. 
The data for hydrogen were obtained with de- 
scending pressures only. This showed that there 
was no hysteresis effect. Data on the compressi- 
bility of isolantite were not found in the litera- 
ture. An estimate of the error due to the cem- 
pressibility was made by using data given for 
soapstone, a similar material. The displacement 
of the plates at 350 atmospheres pressure was 
estimated to be 0.3 percent of the plate spacing, 
causing a change in capacity of less than 0.05 
percent. The values so obtained are slightly lower 
than those of Tang! and of Broxon for pressures 
above 80 atmospheres, and the slopes of the 
curves decrease slightly as the pressure increases 
for both air and hydrogen. According to the 
Clausius-Mosotti equation: 


[ (e—1)/(e+2) ](.M/d)=constant; 


e=dielectric constant; M=molecular weight; 
d=density. The Clausius-Mosotti function 
[(e—1)/(e+2)]/d, was calculated for each pres- 
sure and is shown in the tabulated results. In 
these results, the change in the Clausius-Mosotti 
function with pressure for air is small and irregu- 
lar, while for hydrogen the function decreases 
with increase in pressure. 

The curvature of the pressure-dielectric con- 
stant line for air may be due either to an effect 
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of the high frequency at high gas densities, or to 
some irregularity in the operation of the capacity 
meter. Jordan, Broxon and Walz,* using the ap- 
paratus of Broxon‘ showed that the change of 
dielectric constant of air with frequency was 
very small over the pressure range 1 to 170 
atmospheres and a frequency range of zero to 70 
kilocycles. The experimental results for air 
shown in Fig. 3 fall below those of Tangl': ? and 
of Broxon* and above those of Occhialini.* The 
frequency used in this work was 2.5 megacycles 
while Talbot’ working at a pressure of one at- 
mosphere found no appreciable change in the 
dielectric constant over the frequency range of 
0.025 to 1.100 megacycles. The frequency and 
some of the pressures used in this experimental 
work are beyond the range of the above de- 
terminations. 

Broxon‘ has compared the results of dielectric 
constant determinations found in the literature- 
and has concluded that those made at high fre- 
quencies are not consistent among themselves 
while those at low frequencies are in agreement. 
Cagniard® pointed out that many of the investi- 
gators using high frequency methods had diff- 
culty with instability of their oscillating circuits. 
No instability was encountered in the apparatus 
described in this paper; but, as has been pre- 
viously mentioned, other factors may have 
caused apparent changes in capacity. Further 
research in the measurement of capacities at high 
frequencies is desirable because the heterodyne 
capacity meter, using high frequency oscillators, 
has characteristics which make its use advisable, 
and sometimes necessary for many practical 
applications of capacity measurements. 

’F. L. Talbot, Thesis in Physics, Catholic University of 
America, Washington, D. C., 1928. 


‘LL. Cagniard, Ann. de Physique 9-10, 460 (1928). 
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Errata: Ionization Area of He and Bethe’s Theory 


D. L. Wesster, W. W. HANseN anv F. B. Duveneck, Stanford University 
(Phys. Rev., 43, 839 (1933)) 


CCORDING to Bethe’s'! theory of ioniza- 
tion by impact the area effective for 
ionizing a K shell at any voltage U times the 
ionization voltage Vx should be 
re 2b. 4U 
&( U) =—- — In — 
Vex U B 
provided U is large enough, with 6 and B pure 
numbers depending on the atomic number but 
not on U. In reporting data on Ag and looking 
for evidence on how large U must be for testing 
this theory, we compared this equation with 
Smith’s? data on He. For this purpose we 
plotted a graph using the } and B given by 
Bethe’s Eq. (54b) for the ionization of the H 
atom, namely 0.285 and 0.048, respectively. In 
writing the paper, however, we inadvertently 
recorded for b the value 1, which belongs to the 


H atom, but for the sum of the ionization and 
excitation areas, along with B= 0.638. 

The fact that our graph does not fit our 
recorded } and B has just come to our attention 
in a paper by Liska; and in looking up our notes 
we have found an arithmetical error also. To 
correct this, our graph should be lowered by 14 
percent. This change, however, does not change 
our conclusion that Bethe’s equation holds well 
for He at high voltages but not at all in our 
range U7. Likewise this change has no effect 
whatever on our data or our comparisons with 
theories for Ag. For He the best comparison with 
Bethe’s theory is that of Liska, who sets B at 
0.048 and then treats } as an empirical constant, 
finding a value 0.291, very near Bethe’s 0.285. 


'H. Bethe, Ann. d. Physik S$, 325 (1930). 
?P. T. Smith, Phys. Rev. 36, 1293 (1933). 
3]. W. Liska, Phys. Rev. 46, 169 (1934). 


Proceedings of the Metropolitan Section of the American Physical Society 


MEETING OF MarcH 22, 1935 


HE second meeting of the Metropolitan Section of the American Physical 
Society for the season 1934-1935 was held on Friday, March 22, 1935, in 
Room 301 of the Physical Laboratory of Columbia University in New York 
City. An afternoon and an evening session were held. The presiding officer was 


Professor H. C. Urey. 


The following papers were presented: 


Afternoon Session 


Prediction of Atomic Energy States, R. F. BACHER 
Measurement of Nuclear Magnetic Moments from Hyperfine Structure Data, L. P. 


GRANATH 


Nuclear Moments from Molecular Beams, |. I. Rant 
Isotope Effect in Polyatomic Molecules, W.S. Benepict 


Evening Session 


Methods of Separation of Isotopes, H. C. Urey 
In the afternoon session the following officers were elected for the season 


1935-1936: 
Chairman 
Vice Chairman 
Secretary-Treasurer 


Elected members of the 
Executive Committee 


H. C. Urey 

A. C. G. Mitchell 
W. S. Gorton 

H. D. Smyth 

M. W. Zemansky 


W. S. Gorton, Secretary-Treasurer 


699 








MAY 1, 1935 


PHYSICAL 





REVIEW VOLUME 47 


LETTERS TO THE EDITOR 


Prompt publication of brief reports of important 
discoveries in physics may be secured by addressing 
them to this department. Closing dates for this 
department are, for the first issue of the month, the 


The Hylleraas Method Applied to Lithium 


It has been shown by the work of Hylleraas' on He and 
of the authors? on H; that two-electron problems can be 
satisfactorily solved by the variation method with wave 
functions which provide for the expression of polarization 
effects (the instantaneous distortion which each electron 
produces upon the orbit of the other). A series of terms is 
prepared, starting with a leading term of approximately 
the form which would be a solution of the problem if the 
effects of electronic repulsions could really be replaced by 
those of static screening charges. Other terms are derived 
from this by multiplying it with simple combinations of 
powers of the coordinates of the individual electrons, and 
of their mutual separation. The final wave function is a 
linear combination of these terms, with coefficients deter- 
mined by the variation process. 

We have extended this method to a three-electron 
problem, the ground state of Li. Unfortunately, the pres- 
sure of other work will force us to delay the presentation of 
a full analysis and report upon the problem; we wish here 
to communicate the bare results obtained. We give below 
the total energy as computed in atomic units (e?/a9 = 27.08 
volts). Using a wave function constructed from the same 
series, after removing one electron, we have also computed 
the energy of the ion; we give here the difference between 
these two quantities, as the ionization energy. For com- 
parison, we include the results of Wilson* (the best pre- 
vious calculation), and the observed values.‘ 


Energy of Lithium Atom 


Total Ionization 

e@/ae e*/ae volts 
Experimental 7.4783 0.1981 5.364 
Wilson 7.4192 0.1965 5.320 
James and Coolidge 7.4757 0.1975 5.347 


These results were obtained with a series of 15 terms, of 
which 5 involved distances between electrons of the K-shell 
and that of the Z-shell. Although the mutual polarization 
of the K electrons is undoubtedly greater, we contented 
ourselves with two terms, taking account of this, because our 
chief interest lay in the ionization energy, and it soon be- 
came clear that this would not be affected if the wave 
function contained small errors referring exclusively to the 
inner-shell electrons. Such errors raise the energies of the 
atom and of the ion by practically the same amounts. 
There is no reason to doubt that it would be easy to get a 
considerable improvement in the result for the total energy 
by putting in two or three additional simple terms involving 
mutual polarization of the K electrons. The residual error 
in the ionization energy seems much more obstinate. A 
survey of the probable improvement obtainable from cer- 


twentieth of the preceding month; for the second 
issue, the fifth of the month. The Board of Editors 
does not hold itself responsible for the opinions 
expressed by the correspondents. 


tain untried terms whose effect can be predicted with some 
assurance does not promise more than 0.007 volt. The 
success of the method is therefore comparable to the case of 
the first excited level of parhelium,® where also there is 
involved the reaction between 1s and 2s electrons of the 
same spin. 
HuBeErT M. JAMEs 
ALBERT SPRAGUE COOLIDGE 
Research Laboratory of Physics, 
Mallinckrodt Chemical Laboratory, 
Harvard University, 
April 4, 1935. 
1 E. A. Hylleraas, Zeits. f. Physik 54, 347 (1929). 
2H. M. James and A. S. Coolidge, J. Chem. Phys 
* E. B. Wilson, Jr., J. Chem. Phys. 1, 210 (1933) 


‘A. Ericson and B. Edlén, Zeits. f. Physik 59, 656 (1930) 
* E. A. Hylleraas and B. Undheim, Zeits. f. Physik 65, 759 


1, 825 (1933) 


1930). 


Analysis of Data on Radioactivity Induced by Neutron 
Bombardment 


Though 47 elements are reported to be activated by 
neutron bombardment the literature gives the isotope 
involved and the type of transmutation for only 27. These 
27 are classified in Table I. Group I of this table contains 


Tassie I. Elements activated to radioactivity by neutron bombardment 
Classified according to type of transmutation 


I II Ill 
SmmMpLe CAPTURE PROTON EMISSION \LPHA-EMISSION 


Z MM-Z Z M M-Z Z Ww M-Z 

Na 11 23 12 Meg 12 24 12 I 9 19 10 
K 19 41 22 Al 13 27 i4 Al 13 27 14 
Vv 23 $1 28 Si 14 28 14 P 5 31 16 
Mn 25 55 30 P 15 31 16 ( 7 35 8 
Cu 29 63 34 S 16 32 16 ™ 1 45 24 
Cu 20) «65 6 Ca(?) 20 40 20 Mn 25 55 ti) 
As(?) 33 75 42 Cr 24 52 28 4 59 32 
Br 35 79 44 Fe 26 56 0 All Z’' sand M 
Br 35 81 46 Zn w 64 44 In ever ne 
I 53 127 74 Zn 1 66 36 27 cases, M-Z i 

. - even 
All Z’sand M's odd. 8 of the 10 Z's and 

M's even 
(?) The two cases indicated with this mark in this table, Table II and 


Fig. 1 are somewhat doubtful. 


those activated by simple capture of a neutron; group II 
those activated by capture of a neutron with emission of a 
proton; and group III, those activated by capture of a 
neutron with alpha-emission. Z (the atomic number), M 
(the mass number), and M-Z are also recorded for each 


parent isotope. The following facts are apparent: 
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(1) M-—Z is even for all known parent isotopes. 

(2) M and Z are odd numbers for all parent isotopes ac- 
tivated by simple neutron capture or by capture with alpha- 
emission. 

(3) M and Z are even numbers for 8 of the 10 cases with 
proton emission, the two exceptions being Al and P. 

Table II shows all the known active products produced 








(2) (3) (5) (6) 
(1) Atomic Mass FINAL RELATIVE 
RapIoO— NUMBER NUMBER (4) STABLE ABUNDANCE 
ELEMENT (Z) (M) M-Z IsoTOPE (%) 

N 7 16 9 *O 16 99.76 

Na il 24 13 * Mg 24 77.4 

Mg 12 27 15 * Al 27 100 

Al 13 28 15 * Si 28 89.6 

Si 14 31 17 > a 100 

P 15 32 17 _ 96 

K(?) 19 40 21 *'Ca 40 97 

K 19 42 23 \Ca 42 08 3rd) 
V 23 52 29 Cr 52 81.6 

Mn 25 56 31 * Fe 56 95 

Cu 29 64 35 Zn 64 50.4 

Cu 29 66 37 \Zn 66 27.2 (2nd) 
As(?) 33 76 43 Se 76 9.5 rd) 
Br 35 a0 45 Kr 80 2.45 (5th) 
Br 35 8&2 47 Kr 82 11.79 (3rd 
I 53 128 75 X 128 2.30 (7th) 


* These six elements compose over 93% of the earth's crust 
(?) Doubtful. 


with their values of Z, M and M-—Z. We see that: 

(4) M-Z is odd in every instance (a necessary condition 
since M-Z is even for each of the parent isotopes). 

5) M is even and Z odd for every active product except Mg 
27 and Si 31, the transmutation products of the Al and P 
mentioned in (3). 


Assuming that all these active products are beta-rayers 
(as stated in the literature) it is true without exception for 
both the heavy radioactive elements and for those pro- 
duced up to the present by neutron bombardment, that 
even isotopes of odd elements emit beta-rays. 

In the fifth column of Table II are listed the isotopes 
resulting from the beta-disintegrations. Comparison with a 
table of known stable isotopes’ reveals that 


6) Every one of these 16 isotopes is a known stable isotope. 
The sixth column of Table II records the relative abundance 
in nature of each of these isotopes, expressed as a percentage 
and also in order of abundance, where necessary. This 
column shows that 

7) In all known cases of Z less than 30 the final stable 
product is the most abundant isotope (neglecting of course the 
isotopes of lesser abundance when more than one occurs for 
a given element, as is the case for Ca and Zn); 

(8) The list of final stable products includes six elements 
which constitute over 93 percent of the earth's crust.2 These 
six are: O, Mg, Al, Si, Ca and Fe. 

(9) Five of these six abundant elements, namely, Mg, Al, 
Si, Ca and Fe, are activated by capturing a neutron with 
emission of a proton; and the most abundant of all, oxygen, 
has not been found to be activated. The fact that these five 
isotopes return to their original form may be a partial 
explanation of their great abundance. 
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10) Another fact which emerges is that in every case 
except K 41 the isotope activated is the most abundant one. 
However, this may not be significant, because of the greater 
difficulty of detecting the activation of isotopes present in 
minute proportions. 

Fig. 1 isa plot of M-—Z versus Z, to represent the neutron 
transmutations of elements of atomic number less than 36. 
Five others above atomic number 36 are not shown in this 
figure. Solid circles represent the parent substances; open 
circles represent the radioactive substances produced. A 
single vertical line represents a simple capture; a right- 
angle path represents capture with proton emission; a 45- 
degree path represents capture with alpha-emission. 
Though there is considerable evidence of order in the 
transformations thus far analyzed, the results are as yet 
too fragmentary to permit of complete generalization. 

In addition to the isotopes listed in Table I there are two 
others, Ir and Au, which are known to become radioactive 
by simple capture. These are not included in the list be- 
cause the isotopes responsible are not given in the literature 
However, since their atomic numbers are odd, it is to be 
expected that they have only odd isotopes. If so, the rules 
given for the elements in group I of Table I apply also to 
these two transmutations. 

Carvin N. WARFIELD 

Woman's College of the 

University of North Carolina, 
Artuur E. Rvuark 

University of North Carolina, 

April 3, 1935. 
! Bleakney, Mueller and Barton, Rev. Sci. Inst. 6, 59 (1935 


? Rutherford, Chadwick and Ellis, Radiations from Radioactive Sub 
stances, p. 526 
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The Low Frequency Double Vibrations of the 
Deutero-Ammonias 


In the vibration spectrum of NH, the low frequency 
parallel band (»;) consists of two close components, due 
to the double minimum in potential energy with respect to 
the coordinate along the symmetry axis. If the form of this 
potential function were completely known, the distance 
between the two minima could be obtained from the 
separation of any pair of vibrational components, as 
Dennison and Uhlenbeck have shown.' This, together 
with the moment of inertia indicated by the fine structure, 
would make possible a complete determination of the 
interatomic distances from measurements on one band. 
Conversely, further information from the spectrum may 
serve to define more precisely the potential constants. 
Hence the measurement of these bands for the deutero- 
ammonias is of particular interest. 

We have observed pairs of absorption maxima repre- 
senting the zero branches of »; for each of the four varieties 
of ammonia. These converge rapidly with increasing mass, 
but not as rapidly as would be predicted by the simple 
theory assuming infinitesimal amplitudes. This suggests 
that the effects of interaction are important in determining 
the separations. The two hybrid molecules are slightly 
assymetrical, of course, and would require some modifica- 
tion of the theoretical treatment. The positions of the 
observed absorption maxima in cm™ and their separations 
are as follows: 
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The molecular dimensions can best be determined by a 
comparison of the moments of inertia of NH; and ND, 
as obtained from an analysis of the rotational structure of 
their parallel type bands. These observations are now 
almost complete and will be published shortly. 

We are indebted to Professor Taylor of Princeton Uni- 
versity for a sample of very pure ND,. An earlier mixed 
sample was supplied by Mr. Wells A. Webb, through a 


grant from the Faculty Research Fund. 
E. F. BARKER, 


MarceL MiIGeorre* 
University of Michigan, 
April 17, 1935. 


! Dennison and Ublenbeck, Phys 
*C. R. B. Fellow 


Rev. 41, 313 (1932) 


Nuclear Reactions and Their Classification by Atomic and 
Isotopic Numbers 


Since 1915 the writer has used the isotopic number (J) 
and the atomic number (Z) to exhibit the relations of 
atomic species and their transformations, and these 
variables have been found to give much more simple 
diagrams than those in which the atomic mass (M) and (Z) 
are used. The relations between some of the dependent 
variables involved have been presented by Harkins and 
Madorsky.' They give the formula of any nucleus as 
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Fic. 2. Nuclear reactions arranged according to the projectiles used: neu- 
trons, protons, deuterons and alpha-particles. The transformations due to 
-rays are given in Fig. 1. A transformation of ',;B" of the type +a—p to 
*,C™, which may be radioactive (total change +a—p—e) was omitted by 
mistake. Certain uncertain transformations have been omitted until better 
evidence is obtained. An apparently partly empty section around neon is 
due to the fact that little work has been done with this element. Trans- 
formations of the type +d—n—- are seen to be very similar to those of 
+a—n—rx. With neutrons either two or three steps of the transformation 
+n—a—e occur frequently. With protons the most common change is 
+p—a. Here the positive sign indicates that the particle is added, the 
negative that it is emitted, and p represents the proton, ¢ the negative 


electron and x the positive electron. 


Thus if » represents a neutron the formula becomes 
pn)zn1 (1) 


and this is supposed to represent any nucleus, since while J 
is in general positive it was found to be negative for the 
proton. This Eq. (1) is a sufficient definition of the isotopic 
number. However, let 


M —Z=WN (total number of neutrons) 
then? 
=M-—2Z=N-—Z=2N—M. (2) 


Table I shows that the variables Z and J give a remark- 
ably simple representation of nuclear reactions. This table 
does not give all of the individual reactions, but presents 
practically all of the well-established types. The projectiles 
for which the reactions are given are the positive and nega- 
tive electron (x*, 8) the neutron (nm), the deuteron (d), 
the a-particle, and a y-ray or photon. It is seen that all 
transformations in this type of diagram are represented by 


horizontal or vertical lines, or by straight lines of slope 
—1lor —2 (see key). 

Atomic species which are represented by open circles are 
in an excited state due to the energy which has been ab- 
sorbed on account of the capture of the projectile, while 
those represented by triangles exhibit the phenomenon of 
radioactivity. 

In the preparation of this chart use has been made of the 
proof by Harkins and Gans? that collisions in which the 
projectile is not captured do not in general lead to disin- 
tegration. 

A chart similar to Fig. 1, but in which the reactions due to 
the different projectiles are arranged one above the other, 
but with common lines for (Z) shows that all of thereactions 
given by a single type of projectile are of relatively few 
simple types (see Fig. 2). A general diagram in which all of 
the individual reactions are plotted is not too complicated 
for use, especially if it is compared with the chart for the 
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individual projectiles. Section II, Fig. 1 represents two re- 
actions in which deuterons are involved, section II gives a 
behavior characteristic of atoms of higher atomic number 
in which a neutron is added and the nucleus does not 
immediately disintegrate, but exhibits radioactivity. 
WittiamM D. HARKINS 
University of Chicago, 
April 12, 1935. 


1W. D. Harkins and S. L. Madorsky, Phys. Rev. 19, 136 (1922) 
2W. D. Harkins, Phil. Mag. 43, 315 (1921) 
*W. D. Harkins and D. M. Gans, Phys. Rev. 46, 397 (1934) 


The Seasonal Variation of Ionization in Region F, of the 
Ionosphere 


In his recent letter' in this journal Dr. E. O. Hulburt 
refers to the measurements of noon maximum ionization in 
Region F, and states that there appears to be a discrepancy 
between the results obtained in different latitudes. It is my 
object in this letter to show that by means of a different 
interpretation of the published data the discrepancy can be 
satisfactorily removed. 

If we reject the contention of Kirby, Berkner and Stuart? 
that critical frequency measurements merely indicate ab- 
sorption phenomena and accept an interpretation® of their 
electron-limitation significance in terms of a pronounced 
seasonal variation of molecular temperature, we are obliged 
to conclude that noon Region F; maximum ionization in 
temperate latitudes does not follow the expected seasonal 
variation and is actually slightly higher on a winter noon 
than on a summer noon. Dr. Hulburt correctly points out 
that if my interpretation on these lines is accepted there is 
a discrepancy between this result and other evidence. He 
refers in particular to the ratio (1.5 to 1.8) for summer to 
winter noon ionization, quoted by me in a general iono- 
spheric discussion in London, and which was, in fact, 
calculated from Dr. Hulburt’s own values of short wave 
“skipped distances."’ The discrepancy can therefore be 
narrowed down to the difference between the interpretation 
of critical-frequency measurements made with medium 
wavelengths at short distances and with short wavelengths 
at large distances. It means, briefly, that my interpretation 
of the local ionospheric measurements indicates that the 
maximum Region F; noon ionization is slightly less in 
winter than in summer, whereas Dr. Hulburt's ‘skipped 
distance” data indicate, as he himself has shown,‘ that 
Region F; ionization is 1.5 to 1.8 times as great on a sum- 
mer noon as on a winter noon. 

It would be improper for me to question the accuracy of 
Dr. Hulburt’s measurements of ‘‘skipped distances” were 
it not for the fact that other measurements of similar 
character, made by C. R. Burrows® in America, yield en- 
tirely different results which are in good agreement with 
the local ionospheric measurements, in that the measured 
critical frequency is less on a summer noon than on a winter 
noon. It is true that Burrows interprets his results in terms 
of absorption-limitation, but if we regard both local and 
long-distance measurements as referring to electron- 


THE EDITOR 


limitation phenomena (as Dr. Hulburt and I agree that 
they should be) there is accord between the conclusions 
drawn from both sets of measurements. 

It may not be superfluous to add that this discussion is 
restricted to noon measurements and not to the maximum 
value of ionization density which may be experienced 
throughout the whole of a summer or winter day. 

E. V. APPLETON 

Halley-Stewart Laboratory, 

King’s College, London, 
March 23, 1935. 

' Hulburt, Phys. Rev. 47, 422 

? Kirby, Berkner and Stuart, 
(1934). 

* Appleton, Phys. Rev. 47, 89 (1935). 


‘ Hulburt, Phys. Rev. 31, 1018 (1928); 35, 240 (1930) 
* Burrows, Proc. Inst. Rad. Eng. 19, 1634 (1931 
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Nuclear Magnetic Moment of Na** 


Through the courtesy of Professor Hartree we were 
supplied in the summer of 1933 with s.c.f. functions for 
Na* and with series electron functions for 3s, 3p, 4p 
calculated by J. McDougall in collaboration with him. 
We computed then by means of these functions the ex- 
pected ratio between the hyperfine structure interval 
factor a in cm™ and the nuclear moment yu in nuclear 
magnetons as well as the absolute value of the doublet 
splitting Av for the 3p and 4) terms. 

The theoretical values of the gross doublet splittings 
disagree with experiment by roughly the same large factor 
by which the Hartree function result for the h.f.s. disagrees 
with the Landé-Fermi-Goudsmit formulas. These dis- 
crepancies were brought out in the symposium on nuclear 
moments! in June, 1934, and it was emphasized that in 
view of them one cannot be sure of the theoretical con- 
clusions about the values of nuclear moments for any but 
the simplest atoms. 

For 3s, by taking into account only the energy of the 
series electron in the central field, the effective quantum 
number m*=1.75 as compared with 1.63 experimentally. 
Correcting for interaction with inner shells as has been 
done by McDougall, we obtain a theoretical n* = 1.68. The 
difference between this and 1.63 is small and the difference 
between the Landé-Fermi-Goudsmit result and the 
theoretical one is therefore significant. For 3p the theo- 
retical n* = 2.26 and the experimental = 2.12. It is vital to 
orthogonalize the series electron functions to the core 
functions in computing » and Av because 1/r*° and 
(1/r)(dV/dr) are sensitive to the amplitude of the functions 
at small r. The factor due to this correction is approxi- 
mately 14 for 3p. 

The effects discussed should similarly be sensitive to 
perturbations by configurations involving excitation of core 
electrons. A calculation of the first-order effect of a pertur- 
bation of (2s)?(3p) by (2s)(3s)(3p) changes the theoretical 
h.f.s. value by about 8 percent and is not sufficient to 
matter materially. Similarly the effect of (2p)°(3p)? on 
(2p)*(3p) is appreciable but not sufficient to bring about 


agreement. Effects on h.f.s. due to the first power of 
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coefficients of the perturbing configuration disappear by 
using Fock's field.’ 

The usual computation of nuclear g factors is not neces- 
sarily wrong, but it is fair to say that it is uncertain on 
account of the sensitiveness of the theoretical values to 
perturbations.‘ We intended to postpone publication until 
sufficiently detailed calculations were made to ascertain the 
g factors. Professor J. H. Bartlett kindly informed us that 
he and his co-workers® have made related calculations for 
K using the Hartree field and for Na using Fock’s field, and 
it seemed of interest to report the present state of our 
calculation. 

The results of Shoupp, Bartlett and Dunn for 3s, 3p are 
easily extended to 4s, 4. We summarize in Table I their 
values, our Hartree field values, and our Fock field values 
for 4p. The ratio w;,2(0)/¥42(0)=4.75 by using Fock’'s 
functions, and 4.28 by using Z;Z¢/n**. For Hartree’, Fock’ 
it was supposed that the ratio of the h.f.s. and the multiplet 


TABLE I. 


3s 3 Paz 4pPar 

a Experiment 0.029 8.3 «10°-«7.5 X10" 1.87 x10 
Hartree 5.85 22 (20) 10.4 
Fock 2.5* 5.1* (4.6) 3.1 

“ Goudsmit-Fermi 2.0 2.8 (2.6) 2.0 
Hartree’ 3.3 (3.0) 2.4 
Fock’ 3.1* (2.9) 2.3 
Hartree 2.6 cm™ 1.2 cm 

Av Fock 10.7* 4.0 
Experiment 17.2 5.49 


frequency separations is given correctly by the respective 
functions of Hartree and Fock. The ratio of the two 
frequency differences as obtained from experiment was then 
equated to its theoretical value which involves u and hence 
u was determined. This amounts to defining an effective Z; 
by Z;= —(dV/rdr)/(1/rP) where V is the potential of the 
central field, and by using this Z; in Goudsmit's® Eq. (6). 
The Hartree and Fock values of Z; are 8.91 and 8.54. The 
first of these was communicated by us to Ellett and 
Heydenburg.’ The Z; and » thus determined are insensitive 
to progressive errors in the computation of the wave 
function and depend essentially on its shape for small r. 
They may be expected to be nearly the same for any 
central field calculation. The values of u obtained by means 
of this Z; will be correct if there exists a central field which 
suffices for the discussion of the hyperfine and the ordinary 
fine structures; they may be wrong if there are pertubations 
affecting the multiplet and h.f.s. structures unequally. 

The values marked* were kindly supplied to us by 
Shoupp, Bartlett and Dunn. 

Values of » obtained*: * from 6s, 6p, 7p, 8p of Cs check 
each other satisfactorily. Measurements of Wood and 
Fortrat on the principal series of Na disagree with the 
Landé gross doublet formula while for K, Cs the formula 
holds well. There is thus some evidence that Na is a poor 
element for testing the theory. The difference between 
uw (3p) and uw (4p) may be due, however, partly to experi- 
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mental error as is seen by comparing values in each row of 
the 3p column. 


College of the City of New York, L. A. Wits 
University of Wisconsin, G. BReIr 


April 14, 1935. 


1G. Breit, Phys. Rev. 46, 340 (1934). 

2? J. McDougall, Proc. Roy. Soc. Al38, 550 (1932). 

* Chr. Moller and M. S. Plesset, Phys. Rev. 46, 618 (1934). 

*G. Breit and I. I. Rabi, Phys. Rev. 46, 231 (1934). 

*W. E. Shoupp, J. H. Bartlett, Jr. and C. G. Dunn, Phys. Rev. 
this issue. 

*S. Goudsmit, Phys. Rev. 43, 636 (1933). 

’ A. Ellett and N. P. Heydenburg, Phys. Rev. 46, 583 (1934). 

*L. P. Granath and R. K. Stranathan, Phys. Rev. 46, 316 (1934); 
Bull. Am. Phys. Soc. April 25-27 (1935); Paper 80. 

*N. P. Heydenburg, Phys. Rev. 46, 802 (1934). 


The Magnetic Moment of the Na** Nucleus 


From the most recent experiments with sodium’ it 
has been concluded that the nuclear spin J =3/2 and that 
the total h.f.s. splitting' for the 3 *S, state is 0.0583 cm“, 
for the 3 *Ps,, state 0.0050 cm™, and for the 3 *P, state 
0.0083 cm. Ellett and Heydenburg, using formulae 
developed by Goudsmit® and Fermi and Segré,’ have 
calculated values of the nuclear magnetic moment. From 
the 3 *S, state, they find «= 2.02 nuclear magnetons; while 
from the 3 ?P states, they obtain «= 2.3—2.6 n.m. Hartree 
wave functions were used by Wills and Breit,* who obtained 
values »=5.8 and w=22.6 from the S and P states, 
respectively. Since, however, their calculations show the 
multiplet splitting for the *P state to be 2.60 cm™, as 
compared with the observed value of 17.6 cm™, it is fairly 
obvious that better wave functions than those of Hartree 
are needed for accurate calculations. It is the purpose of 
this note to investigate whether or not the functions 
recently published by Fock and Petrashen® will serve the 
purpose. 

The h.f.s. separation of a *S state is A(s)=(8#/3){(2/ 
+1)/I}uuo*(0); that of a *P a2 state is A*Py_= (8/3){ (27 
+1)/I}uo(1/r?). The gross structure separation for the *P 
state is’ 6= (3yuo*/hcag*)(1/r)(dU/dr), where U(r)=—11/r 
+ V(r), and V(r), is tabulated by Fock and Petrashen. 


We find the following results: lim (/3./r) = 2.86; (1/r*)s» 


di r— 0) 
=0.144; (1/r)(dU/dr)s, =1.23; from *S, «4 =2.5 n.m.; from 
2P, u=5.1 n.m.; and 6=10.7 cm. The values of 


(1/r)(dV/dr) are given in Table I. 


rasie I. 
1 dt —idV —idV —idV 
r - r r r - 
r dr r @ r dr ,r @ 
0.00 0.24 147 0.60 23.1 14 4.43 
02 2400 .26 123 65 19.9 1.5 2.84 
.04 1700 28 105 .70 17.2 1.6 2.37 
06 1300 30 90.0 75 15.0 1.7 1.99 
OR 987 32 78.1 80 13.2 1.8 1.69 
10 747 44 68.5 BS 11.6 19 1.44 
12 571 36 60.7 90 10.2 2.0 1.24 
14 441 38 56.2 95 9.12 2.2 J 
16 345 40 49.2 1.0 8.04 2.4 
18 273 45 39.2 1.1 6.37 | 10 
20 220 50 32.4 1.2 5.17 r 
22 179 55 27.1 1.3 4.20 
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It is seen that the introduction of the Fock functions 
serves to improve the calculated value of the *P multi- 
plet splitting very much. If one assumes the ratio 
(1/r*)s : (1/r)(dU/dr),, to be approximately the same for 
the Hartree and Fock functions, then this means that »é 
is also approximately the same. For the Hartree case, 
ud = 58.8; for the Fock case us = 54.6. If we suppose that 
this latter value is not very different from that which would 
be obtained with exact wave functions, and that the 
present theory of multiplet separations is correct, then it is 
permissible to use the experimental value of 4 and to solve, 
and one obtains «= 3.1 n.m. 

The agreement between the value «1 =2.5 n.m. from the 
25 state and the estimated value un» =3.1 n.m. from the *P 
state seems to indicate that the true value of the magnetic 
moment of the Na™ nucleus lies in this neighborhood. It 
may be remarked that our method of estimation seems 
capable of yielding good results for p electrons even with a 
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rather poor wave function (such as Hartree’s), but that 
accurate results for s electrons will probably be obtained 
only with quite accurate wave functions, as may be seen 
by comparing the values u=5.8 (Hartree function) and 
u=2.5 (Fock function). 

W. E. SHoupp 

J. H. Bartietrt, Jr. 

C. G. DuNN 
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